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ABSTRACT 


Atomic multiplets are treated by wave mechanics, without using group theory. 
In part 1 Hund’s scheme for multiplet classification is derived directly from theory. 
Part 2 is devoted to the computation of the energy distances between multiplets, and 
comparison of these distances with experiment in some typical examples. There is no 
treatment of the separations between the various terms of a multiplet, since that has 
been done elsewhere, but only between one multiplet and another. It is found that 
Hund’s rule, that terms of large L and S values lie lowest, has no general significance; 
the present theory leads to the same results as the rule when it is obeyed experi- 
mentally, but many cases which were exceptions to that rule are in agreement with 
the theory. The method of calculation of multiplet distances is described in sufficient 
detail, with the necessary tables of coefficients, etc., so that further checks with experi- 
ment could easily be made. 


HE theory of complex spectra is treated in this paper by the method 

of wave mechanics. The results of the calculation may be divided into 
two parts: first, the classification of the terms into multiplets; second, the 
energy values of these multiplets. The first part contains no new results of 
physical interest, for it leads precisely to Hund’s' scheme of classification, 
and uses almost the identical steps that Hund uses. Its value lies in the fact 
that this well-known scheme is shown in an elementary way to follow directly 
from wave mechanics. The second part, however, is almost entirely new, 
and it leads to definite formulas for the intervals between the different mul- 
tiplets in spectra, intervals which could previously be considered only 
very roughly from an empirical rule of Hund, which proves to have no gen- 
eral significance. These calculations, of course, give the intervals in terms 
of certain integrals, which we do not calculate, but merely estimate well 
enough to permit some comparison with experiment. The agreement is in 
general fairly good. A third part would be also included, dealing with the 
intervals within the multiplets, produced by the magnetic interaction, were 


1F. Hund, Linienspektren und periodisches System der Elemente, Springer. 
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it not that Goudsmit? has already answered the question by a method closely 
analogous to that used here. 

It will be noted that the objects of the present paper resemble closely 
those aimed at by Heisenberg, Wigner, Hund, Heitler, Weyl, and others,’ 
who employ the methods of the group theory. That method is not used at 
all in the present calculation,and, in contrast, no mathematics but the sim- 
plest is required, until one actually comes to the computation of the integrals. 
This, it is believed, is in itself sufficient justification for paralleling to some 
extent work already done. The simplification is achieved largely by intro- 
ducing the spin at the very beginning of the calculation, rather than later. 
Thus we need only consider antisymmetric wave functions, which can be 
treated very simply as determinants, and can avoid the other symmetry 
characters, with which the other papers have been mostly concerned. The 
process of building up wave functions of the proper symmetry by using 
determinants is not new; it is found in Dirac’s earlier papers, and has been 
used, for example, by Waller and Hartree.* The spin function, however, is 
introduced into the determinant in the present paper in a new and more 
satisfactory way. The results which we obtain, concerning the diagonal 
term of the energy with respect to the antisymmetric wave functions, are of 
the sort found by the group theory, and also by the recent method of Dirac, 
but it seems worth while to derive them in a simple fashion. 

The essentially new results of the present paper, those relating to energy, 
in the second part, come from the fact that we consider the whole degeneracy, 
that coming from orbital as well as spin angular momentum, which most 
of the other papers have failed to do. The present method is essentially 
equivalent to the others in that it gives sums of energies, rather than in- 
dividual energies. When it is applied to the case in which the degeneracy 
with respect to orbital angular momentum is not considered, it leads, as do 
the others, merely to the sum of all terms of a given multiplicity; the writer 
is indebted to Dr. Bloch for pointing out that in that case all the results of 
Heitler can be easily demonstrated by the present method. But we make 
the observation that, byconsidering the whole problem with all its degeneracy, 
the method' of energy sums can be used much more effectively, so much so 
that in most important cases we can get the actual energies of the individ- 
ual terms. In connection with this, it should be noted that in several cases 
the actual perturbation problems have been solved directly, rather than by 
the method of energy sums, obtaining results® that hold even when the mag- 
netic energy is appreciable. 


2S. Goudsmit, Phys. Rev. 31, 946 (1928). 

8 W. Heisenberg, Zeits. f. Physik 41, 239 (1927). E. Wigner, Zeits. f. Physik 40, 492, 883, 
etc. (1927). F. Hund, Zeits. f. Physik 43, 788 (1927). W. Heitler, Zeits. f. Physik 46, 47 (1928). 
H. Weyl, Gruppentheorie und Quantenmechanik, Hirzel, 1928. 

‘I. Waller and D. R. Hartree, Proc. Roy. Soc. A124, 119 (1929). The writer is indebted 
to Dr. Hartree for calling his attention to this paper. 
’ P. A. M. Dirac, Proc. Roy. Soc. A123, 714 (1929). 
* W. V. Houston, Phys. Rev. 33, 297 (1929). J. A. Gaunt, Trans. Camb. Phil. Soc. (1929). 





























THEORY OF COMPLEX SPECTRA 


PART 1. THE CLASSIFICATION OF MULTIPLETS 


1. As a preface, we remind the reader of the theory of multiplets, as 
developed before wave mechanics. This theory considered the interactions 
of the various angular momenta of the parts of an atom, and treated the 
various energy levels resulting from different relative orientations of these 
parts. In a logical development of atomic structure it would follow directly 
after Bohr’s theory of electron orbits; for Bohr went as far as was possible 
on the assumption that each electron moved in a field with spherical sym- 
metry, so that its orientation with respect to the rest of the atom was a 
matter of no concern. We start then, as Bohr finished, with an atom com- 
posed of a number of electrons, each characterised by a total quantum num- 
ber m, and an azimuthal quantum number /, the latter determining the angular 
momentum, and being often regarded as a vector (normal to the plane of 
the orbit, in Bohr’s theory, but this is of no importance). The orientation 
of this vector in space was arbitrary, on account of the spherical symmetry. 
As an additional part of the formulation, we need to consider the electronic 
spin of Goudsmit and Uhlenbeck:each electron has an angular momentum 
(with consequent magnetic moment) entirely apart from the momentum / 
which it acquires by its motion. This additional angular momentum, de- 
noted by s, (which like /, is often regarded as a vector), can likewise be 
oriented arbitrarily in space, without reference to the direction of J. Since 
the values which n, /, and s can take on, and the notation, are rather signifi- 
cant, we observe that / can equal 0, 1, 2, - - -; can be any integer greater 
than or equal to (/+1); and s is always equal to 1/2. The values of m and / 
are denoted by symbols, as 1s, 2s,-- -; 2p, 3p: - -; 3d, 4d,-- -; 4f, Sf-- >; 
referring respectively to the cases where m and / are 10, 20, - - -; 21, 31,-- -; 
32, 42, -- -; etc. Then the structure of an entire atom in a given state can 
be given, as far as Bohr’s theory is concerned, by giving the m and / of each 
electron, denoted as in the example: (1s)? (2s)? (2p)® (3s)?, meaning that 
there are two (1s) electrons, two (2s), six (2p), two (3s). Another state of 
the same atom (in this case magnesium) would be (1s)? (2s)? (2p)® (3s) (3p); 
and so on. These two states would be often abbreviated 3s and 3p respec- 
tively. The spin, having always the same value, need not be mentioned at 
this stage. 

Having described the atom to the degree of accuracy considered by Bohr, 
we next must consider that the angular momentum vectors really are coupled 
together; there are differences of energy depending on orientation, and a 
consequent splitting up of each energy level into a number of different ones. 
In the commonest case, which alone we shall consider, the coupling can be 
considered in several stages. First, the /’s of the various electrons group 
themselves into a vector sum L, which takes on only integral values, 0, 
1, 2,---, as the separate /’s do. Each such arrangement gives a different 
energy, and a different term in the spectrum. The terms with L equal 
respectively to 0, 1, 2,3,---, are called S, P, D, F,- - - terms. Second, the 
spins of the electrons group themselves into a vector sum S, taking on the 
integral or half integral values (for any case, the orientation where all s’s are 
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parallel is allowed; if it is integral, only integral S’s are allowed, and if it 
is half integral, only half integers are allowed for S). The different values of 
S lead to terms of different multiplicities, as we shall see in a moment, and 
different energies; thus S=0 gives singlets, S=1/2 doublets, S=@ triplets, 
etc. Finally Z and S are coupled together, their sum being called J, and being 
integral or half integral to correspond with S. Since J can vary from IL +S | 
to |L—S|, there are (if L>S) 2S+1 different terms so obtained. These 
terms are the various terms of a single multiplet, which is named singlet 
(2S+1=1), doublet (2S+1=2), etc., to agree with the number of terms. 
The energies involved in going from one value of Z to another, or of S to 
another, are generally large; they are (as one sees from wave mechanics) 
electrostatic energies. On the other hand, the energy of orientation of L 
with respect to S is small, coming from magnetic interaction of the parts of 
the atom, so that the different levels of a multiplet (same L and S, different J) 
lie near enough together to be grouped together. If one neglects the magnetic 
interaction as a first approximation, the levels of a multiplet all have just 
the same energy. As to notation, we speak, for example, of *P:, meaning a 
triplet P term with J =2. And if we wish to indicate the state of the electrons 
which produce the multiplet terms, we write it as (1s)? (2s)? (2p)® (3s) (3p) *Pe, 
which would often be abbreviated (3p) *P». 

The process of coupling can also be described in the language of space 
quantization. If one has an angular momentum vector, as /, acted on by no 
torques, so that its direction is arbitrary, it is legitimate to choose any axis 
in space, and consider that the component of the vector along this direction 
can take up any one of a set of values, integrally spaced between the parallel 
and antiparallel orientations. Thus, if m; is the component of / along this 
axis, m,can equal /,/—1,/—2,--- —l. These (2/+1) values are considered 
to denote (2/+1) separate stationary states, all with the same energy. Now 
we follow by this method the process of coupling, beginning with the un- 
coupled vectors, and carrying it through to the case where the /’s are coupled 
to give L, and the s’s to give S, but the magnetic interactions are not con- 
sidered. In the uncoupled case, each electron has its own m; and m,. As 
the coupling forces are introduced, torques appear which make the com- 
ponents of the individual /’sand s’s vary with time,so that the same space 
quantization is no longer possible. But the torques are internal; they cannot 
change the total angular momentum. The total sum ofall the/’s and the sum 
of all the s’s, and their components along the axes, remain constant, and 
hence quantized: }-m, and >-m, are quantized even when the coupling has 
taken its full value. On the other hand, when the electrons are coupled, the 
atom consists of vectors L and S, each separately free, so that these can be 
space quantized, giving components M, and Ms along the axis. We need 
now only identify the Mz and Ms for the coupled system with the }>m, 
and >\m, for the uncoupled one, to get the correspondence necessary to 
show what multiplets appear from any combination of electrons. The exact 
method is easily shown from an example. 
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For illustration, we take the coupling of two p electrons. Each has 
J=1; the first has principal quantum number n, the second n’. For each, 
we can have m,=1, 0, —1; m,=1/2, —1/2; that is, 2(2/+1)=6 different 

—sets (m,, m,), each representing a stationary state. With the combined atom, 
there are 36 combinations of one state of the first electron with one of the 
second. These are partly given in Table I. Two different notations are 


TABLE I. Two p electrons. 




















First notation Second notation =m =m, Notes 
(m 114) (m’11 4) {(m 11) (m’11)} { } 2 1 Exc 
(am 104) (m’11 4) {(m 10) (m’11)} { } 1 1 A 
(m.1 —1 4) (m’11 4) {(m 1 —1) (m’11)} { } 0 1 
(n 11 —4) (m’11 4) {(m’ 11)} {(m11)} 2 0 
(n 10 —4) (m’11 4) {(m’ 11)} {(# 10)} 1 0 
(mn 1 —1 —}) (m’ 11 4) {(m’11)} {(@ 1 —1)} 0 0 
(m 11 4) (m’ 10 4) {(m 11) (m’10)} { } 1 1 A 
(mn 10 4) (n’ 10 3) {(m 10) (m’10)} { } 0 1 Exc 
(m 1 —1 3) (m’ 104) {(m 1 —1)(m’10)} { } =| 1 
etc. 








used in this table for comparison. First we simply give the four quantum 
numbers (m / m, m,) for each of the electrons. For the other scheme of nota- 
tion, which is often convenient in practice, we take advantage of the fact 
that each m, is capable of only the two values 1/2, —1/2: we set up two 
brackets, one containing the symbols (n/m,) of each electron with m, 
=1/2, the second containing the symbols for those with m,= —1/2. Such 
a bracket symbol, in which the arrangement of terms within a bracket is 
of no significance, is equivalent to a state of the first sort. After each symbol 
is given the value of }>m, and of >\m, for that state. Finally, after some of 
the terms, are notes which will be explained later. 

The values of }\m, and >>m, are conveniently plotted by giving }>m, 
as abscissa, }>m, as ordinate, and putting a dot for each state. Thus in 
Fig. 1, a, we give a single p state; and in b, the points for the sum of the two 
(the numbers indicate the number of states with the same )>m,, >>m)). 
But now in Fig. 1, b, each point should also represent the Mz, Ms of one 
of the multiplets existing in the coupled atom. This can be brought about 
in but one way: by assuming the multiplets 'S'P'D*S*P*D, whose 
separate representations are shown below (in Fig. 1, d—i), and whose 
patterns, if superposed, would just give Fig. 2b. Thus it must be that two p 
electrons produce these multiplets. But that is just what our other method 
gives; for two l’s,each equal to unity, can add to give L=0, 1, 2, (S, P, D 
terms); while two s’s, each 1/2, can give S=0, 1 (singlet, triplet). In 
every case, the two methods lead to the same result, just as here. 

The scheme as described above will result, in an atom with many electrons, 
in an enormous variety of terms, comparatively few of which are realized 
in actual atoms. The principle limiting the number is the exclusion prin- 
ciple of Pauli. It cannot be stated in the language of vectors, but only in 
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terms of the space quantization just described. It is: no two electrons can 
have identical sets of numbers (n/m, m,); further, two combinations of 
electrons which differ only in the interchange of the quantum numbers of 
two electrons are to be treated as identical. Stated in terms of the bracket 
method of writing the electronic combinations, this means that no two sets 
of quantum numbers (m” / m,) within a bracket can be identical; and that 
the order of arrangement within a bracket is immaterial. By these restric- 
tions, the number of terms is often enormously decreased. Many points on 
our diagrams are removed, and just the remaining ones must be fitted into 
multiplets, resulting in a much smaller number. For example, in our Table 
I, the two electrons have the same /; hence, if they have also the same total 
quantum number (that is, if =m’), they cannot both have the same set 
of m, and m,. Thus the two terms marked Exc (and many other terms not 
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Fig. 1. (a) Single p electron; (b) Two p electrons; (c) Two equivalent 
p electrons; (d) 'S (e) 4F; (f) "D; (g) 3S; (h) *P; (i) 8D. 


included in the table) must be excluded; and the two terms marked A (and 
many more pairs) are to be considered identical. When the terms not allowed 
are removed, the new pattern in place of Fig. 1b proves to be Fig. 1c; and 
this is the superposition of *P 1S 'D, which then are the multiplets allowed 
with two equivalent p electrons (that is, two electrons, each with /=1, and 
with the same m). As a general thing, the exclusion principle is active only 
when there are equivalent electrons; its many properties, as in limiting the 
number of electrons in a closed shell, are well known and need not be elabor- 
ated. 

2. In quantum mechanics, we meet the problem of complex spectra as 
one step in the approximate solution of the wave equation for the atom. 
There are two fundamental principles which govern the structure of matter: 
quantum dynamics, and the principle of antisymmetry which shows itself 
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in the exclusion principle and the Fermi-Dirac statistics. The first states . 
that only those energy levels are possible which are characteristic numbers 
of Schrédinger’s equation for the system; the second further restricts the 
possible energies to those connected with characteristic functions antisymmet- 
ric in the electrons. Our whole object is to find such characteristic functions 
and numbers. And as before, the first step isto approximate by supposing 
that the separate electrons move in fields of force with spherical symmetry, 
so that there is no tendency for orientation of the various angular momenta. 
There are various ways of doing this; one scheme which gives good results 
is that of Hartree.’ For simplicity in description we shall imagine that scheme 
modified slightly in one detail; according to it, each electron moves in a 
field of force slightly different from the others. We shall neglect the differ- 
ence, assuming that all the electrons move in precisely the same field. And 
this field is to be so chosen as to give the best agreement with the correct 
values even without further corrections. 

An electron moving ina central field of force, according to wave mechanics, 
is characterized by the same two quantum numbers m and / that we have 
previously described. The arbitrary direction of the angular momentum is, 
as before, most conveniently described by quantizing the component m, 
in a fixed direction. Similarly the orientation of the spin is most conveniently 
given by specifying its component m, in the same fixed direction. Thus 
each electron in an atom, in the approximation in which we can neglect the 
interactions of their rotations, is specified by the four quantum numbers 
nlm,m,, just as it was before. All relations of these numbers remain 
unchanged. But now, each electron has a wave function—a function of 
its coordinates (and, as we shall describe presently, of a coordinate representing 
its spin) depending on the numbers n / m, m,, which is a solution of Schré- 
dinger’s equation for a particle in a central field. We can denote the function 
for the ith electron by u(n;/x;), where n;stands for the four numbers ;1; mj; 
m,;, and x; symbolizes the four coordinates (three of position, one of spin) of 
the ith electron. Now it is well known that the product of these functions, for 
all the electrons (1 - - -N) of the atom, gives a function which approximately 
satisfies Schrédinger’s equation. That is, u(m,/x1)u(m2/x2) -- - u(my/xy)is 
an approximate solution. But it is not antisymmetric in the electrons, so 
that it does not satisfy the exclusion principle. To build up an antisymmetric 
solution, we first note that we still have an approximate solution, connected 
with the same energy value, if we interchange any two x’s, obtaining for 
example u(m;/x2) u(me/x1)--- u(my/xy). We still have an approximation 
with the same energy if we make a linear combination of any such solutions. 
Then we can make the one possible combination which is antisymmetric, 
and it will both satisfy the exclusion principle, and will be an approximate 
solution of Schrédinger’s equation. This combination is conveniently written 
as a determinant: 


7D. R. Hartree, Proc. Camb. Phil. Soc. 24, 89 (1928). See also J. C. Slater, Phys. Rev. 
32, 339 (1928) for discussion of Hartree’s method and application of some of present results 
to it. 
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u(mi/x1) u(m/x2) -- - u(m/ xy) 


u(my/x1) u(ny/ x2) - - - u(my/ xy) 


It is obviously antisymmetric, for interchanging, say, x, and x, interchanges 
two columns of the determinant, which by a familiar property merely changes 
the sign. It can be shown that it is the only antisymmetric combination of 
these functions. And it leads at once to the familiar interpretation of the 
exclusion principle. For if two of the functions had the same quantum 
numbers (say ”, = m2, symbolizing equality of four quantum numbers), then 
the corresponding rows of the determinant would be identical (since they 
contain the functions u(m,/x:) =u(m2/x1), u(mi/x2) =u(m2/x2), etc.) and by 
another familiar rule, the determinant will vanish. Thus there is no solution 
corresponding to the case where two electrons have the same set of quantum 
numbers. Further, the determinant treats all electrons alike; hence we cannot 
count as separate two states which differ only by the interchange of the 
quantum numbers of two electrons. Our exclusion principle then coincides 
with the one previously described. 

We now have, corresponding to each set of quantum numbers, or to 
each independent bracket expression of quantum numbers (as given in Table 
I), which is allowed by Pauli’s principle, a single, antisymmetric function 
of the electrons, which is an approximate solution of Schrédinger’s equation. 
Our next task is, beginning with this, to introduce the interaction between 
the various angular momenta, and to try to improve the agreement with 
Schrédinger’s equation, without destroying the property of antisymmetry. 
We use essentially the method of perturbations. The fundamental result 
of this method is that, if we take the approximate but incorrect wave func- 
tion, and compute the matrix of the real energy with respect to this, the diag- 
onal terms of this matrix are good approximations to the actual energy 
values of the problem. The errors remaining are of the order of the square 
of the ratio of non-diagonal to diagonal terms. Since we can easily show 
in our case that the non-diagonal terms here are really small, this method 
will give a good approximation to the energy values. We are then to take 
the real energy operator (involving the interactions between electrons, 
rather than with fictitious central fields), find its matrix with respect to 
the wave functions already determined and take these diagonal terms as 
energy levels. 

There is, however, one case in which our criterion for the accuracy of 
this approximation is not valid. This is the case where a number of terms 
lie close together. Then in the first place one can no longer say that the errors 
are as small as we have assumed; in the second place, since we are generally 
interested in the energy differences between the neighboring terms, we really 
demand a much greater accuracy than usual, to give this difference correctly. 
Thus this case—that of degeneracy—demands special treatment. We see 
that, in our case, we actually meet this difficulty, for the various wave 
functions with the same values of n’s and /’s, but different m,’s and m,’s, 
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have very nearly the same energy, and yet the energy differences, giving 
the separations between multiplet levels, are just what we wish to find. 
Thus the problem of degeneracy is the essential part of the calculation for 
complex spectra. We may describe the situation this way: the non-diag- 
onal energy terms connected with transitions between states with different 
n’s and I's for the electrons are negligible; but those between states with 
the same n’s and I/’s but different m,’s although they may be no greater 
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Figs. 2-5. Figures 2, 3, 4 represent the energy matrix schematically. In Fig. 2, it is with 
respect to the unperturbed wave functions. The terms in the double shaded squares represent 
transitions between two states with the same electron quantum numbers. These terms are all 
significant; those outside the squares are negligible. In Fig. 3 we have made the linear combi- 
nations to the perturbed wave functions, reducing each square to a diagonal matrix. Fig. 4 is 
with respect to unperturbed wave functions, as is Fig. 2; but we take account of the fact that 
there are no components between states of different }-m:, }-m,. In Fig. 4, energy levels are 
arranged according to the values of dm, S-m,. Fig. 5 (see Note 3) represents the matrix of 
angular momentum, in the final perturbed wave functions. Each square represents a separate 
multiplet. In Figs. 3 and 5 the terms are arranged according to multiplets. 


numerically, are not negligible for our purpose. This situation can be de- 
scribed graphically. Suppose we make a scheme for the energy matrix, 
giving stationary states (the incorrect ones that we determine from the 
central field, denoted by the numbers ; - -: - m,y) along the two sides of a 
square array, and putting in the matrix components. Then if we arrange 
together all terms with one set of m’s and /’s, then those with another set, 
and so on, we see that the matrix components in the double shaded regions 
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of Fig. 2 are all important (even the non-diagonal ones); but those in the 
singly shaded regions are negligible. The method of treatment is now this: 
we take a particular square of the matrix, as A, coming from transitions be- 
tween different states with the same n’s and /’s. And, by linear combinations 
of the wave functions connected with these states, we reduce the matrix 
to a diagonal one. This can always be done for a finite matrix. Then we 
do the same thing with each such double-shaded square. The result after 
that is then a matrix of the form shown in Fig. 3. Here the only important 
non-diagonal terms have disappeared; so that we can take the diagonal terms 
as the approximate energy values of the real problem, with assurance that 
the errors are not large. We observe that, after this is done, there will still 
be just'as many wave functions connected with each set of values of the n’s 
and /’s as before. 

In our problem of complex spectra, there is a feature which greatly 
simplifies the calculation. For we shall prove*® that the energy has no matrix 
components connected with transitions in which }\m, or }om, change (if 
we neglect the magnetic energy, as we are doing for the present). To show 
the effect of this graphically, we can arrange the terms, first according to 
the n’s and /’s, as we have done; but under each such classification, we can 
arrange according to Ym, and dom. Then the original energy matrix really 
has the simplified form shown in Fig. 4. And the process of removing the 
important non-diagonal terms reduces to separate linear combinations be- 
tween the functions of each set having a given set of n’s, l’s, Yomi, and 
>om,. Generally the number of functions actually present in one of these 
groups is very small; thus the problem by this method becomes very simple, 
and can be readily carried out. As a result of it, we see that the number 
of states with a given set of m’s and/’s, and a given >.m, and >. m,, remains 
unchanged as we apply the interaction between angular momenta. And this 
is the essential point required in the classification of the terms. 

3. The essential features of the perturbation theory as applied to com- 
plex atoms have been described; and we can now make connections with 
Hund’s scheme for classifying the terms. Each unperturbed wave function 
can be described by its set of (” 1 m,; m, )for each electron (Pauli’s principle 
being actually satisfied by our condition of antisymmetry). Then, just as 
before, we can make a table of all the possible terms with any set of ’s 
and /’s; we can find ym, and > om, for each of the terms; and we can make a 
diagram, as before, plotting dm: against > m, obtaining just the sort of 
point diagram that was shown in Fig. 1, b or c. The small square arrays of 
matrix components, in Fig. 4, are just the components between different 
terms represented by coincident points in Fig. 1, b and c. The components 
between terms lying at different points on the }.m;—)_m, diagram are zero. 
Now we have just proved that, even when the interactions are considered, 
we still have the same diagram representing the number of states with each 
set of })m; and >\m,. To verify Hund’s method completely, we need only 
show that we can correctly break up this diagram into a number of separate 


8 See Note 1 for the proof. 
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ones, as Fig. 1, d—i, each representing a separate multiplet. Obviously we 
can separate the diagram into rectangular arrays of points (that is, sets of 
wave functions with particular values of }>m, and }_m,), but there are other 
conditions that must be satisfied to make these wave functions really repre- 
sent a multiplet. 

In looking for the conditions that a set of wave functions must satisfy to 
represent a multiplet, we first observe that, so long as we neglect mag- 
netic energies, they must all have the same energy value. Now after we 
have made our linear combinations, we shall have perfectly definite wave 
functions, each with a perfectly definite energy value (the diagonal terms in 
Fig. 3). Without proof, it is not obvious that we can pick out sets of functions, 
with different >-m, and >-m,, but all with the same energy. Even if we can 
(and we shall prove that it is possible) ,*the set of functions must satisfyfurther 
conditions. Somehow we must work in the fact that they all represent the 
same vector L, and the same S, but with different orientations in space. 
Now the wave functions corresponding to an angular momentum vector 
with different orientations have perfectly definite relations, which are ex- 
pressed by computing the matrix components of the angular momentum. If, 
then, our group of terms of the same energy really is to represent a multi- 
plet, the matrix of angular momentum must also have the required form. But 
we shall prove that this also is the case. Thus, with these proofs given, we 
see that we have a perfect right to-separate the wave functions into a group 
of multiplets, as Hund does. Since the method of doing this is unique, the 
result must give precisely the same classification of multiplets that Hund 
finds. And since at present we wish only classification, we need not actually 
make the linear combinations, and get the diagonal matrix of H, at all. 

It remains to be proved that the terms, after making the proper linear 
combinations, really can be divided up into groups of terms, having the 
following properties: every wave function of the group has the same energy 
value; and the angular momentum matrix has just the same values that it 
would if the wave functions referred to the vectors ZL and S with their various 
orientations in space. The details of the proof are given in the notes®; but we 
can state the essential features. First we show, by fundamental methods, 
that the total angular momentum matrix can have components only between 
states of the same energy. Since we can show that the angular momentum 
matrix is not diagonal, this proves that there must be at least several states 
of each energy, in order that the angular momentum can have components 
between them. This breaks up the terms associated with a given set of n’s 
and /’s into sets each of the same energy. Next we show that each of these 
sets is just the sort we need for a multiplet. We do this by detailed considera- 
tion of the angular momentum matrix; the essential point being that the 
entire matrix can be determined uniquely from certain commutation relations 
between the x, y, z components of the angular momentum vector, and that 
these relations are just the same for the whole atom as for a single angular 
momentum vector, so that we must have just the same matrix. We then 


® See Note 2. 
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have the complete proof that the wave functions of the atom really break up 
into sets each having all the properties of a multiplet; the Mz—Ms diagram 
of such a multiplet is just as we assumed in the first section; and the result is 
that the method of finding what multiplets arise from any configuration of 
electrons, is precisely the method of Hund. 


PART 2. DETERMINATION OF THE ENERGIES 


1. To find the energy values of the multiplets, the natural procedure is 
to solve the various problems of linear combination encountered above, find 
the correct wave functions, and compute the diagonal terms of the energy 
with respect to them. But we can find a simpler way, a way that even here 
will allow us to omit the actual calculation of the wave functions altogether. 
This is done by the use of the principle that the sum of a number of energy 
values is in some cases not changed by applying a perturbation, even if the 
individual values are. The principle is this: given a set of wave functions 


u;- +--+ u;; make a set of orthogonal linear combinations of them, 1’ -: - 
u;’. Consider the diagonal terms of a matrix, as the energy matrix, referred 
to the original wave functions: say E£, --- E;. Consider also the diagonal 


terms of the matrix of the same function, referred to the new wave functions, 
as E,’---E,’. Then £,\+ --- +£;=E,'+ ---+£,'. Inourcase, we can use 
this as follows. We consider a set of terms, all with the same n’s, I's, dom, 
>> m, (that is, the terms connected with a small square in Fig. 4, or with coin- 
cident points in Fig. 2b and c). Before making our linear combinations, we 
can easily calculate the diagonal terms of the energy. We can add these terms 
for all the wave functions. Then the sum is equal to the sum of the corres- 
ponding energy values after making the linear combinations; that is, it is 
the sum of certain energy values that we wish to calculate. By proper use of 
this method, we can generally get the energies of all the multiplets. 

To make the process clear, we shall illustrate by the case of two equivalent 


Gap electrons, shown in Fig. 2c. We compute the energy of the function 


m,=0, >.m,=2; that is, the energy of the configuration {(m11)} {(m11) }. 
But there is only one function connected with this point in the diagram. 
Hence the sum degenerates to one term. The process of making linear com- 
binations is not necessary here; the corresponding square matrix array in 
Fig. 4 has but one row and column. Thus the energy we have found is also 
the energy of the term after interaction is considered; and comparison with 
(f) shows that it is the energy of the 'D multiplet. Similarly the point 
dom, =1, > m,=1, gives the energy of *P. Again, dom, =1, >m,=0 or 
—1, and }>m,=—1 >>m,=1, 0, —1, should all give the energy of P; and in 
consequence of certain identities between the energy matrix components, 
these actually do give the same result as before. Also, the point >>m,=0, 
>-m:=1, has two wave functions, for the states {(m10)} {(m1 1)} and 
{(m1 1)} {(m10)}. We calculate the energies of these terms, and add. 
And the result should, and does, come out to be the sum of the 'D and *P 
energies. Finally, the point }\m,=0, }.m,=0 has three wave functions. 
We add their energies, and the result should be the sum of the energies of 
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1D, *P, and 'S. Since the first two of these energies are known, we need only 
subtract them from the sum to get the energy of 1S. Thus the energies of 
all the multiplets are computed, without finding any perturbed wave functions 
atall. Itis readily seen that extension of the same scheme will give complete 
information in all cases except where there is more than one multiplet of the 
same kind; for example, if a given configuration of electrons contains two 
2D multiplets, the method will give only the sum of the two terms. To get the 
individual terms, we must carry through the more elaborate scheme of taking 
linear combinations. But such cases are not found among the more impor- 
tant terms, so that in practical use the limitation is not important. 

For actual use, the essential step is to find the energy values referred to 
the original wave functions; then it is simply a matter of adding and sub- 
tracting to get the energies of the multiplets. This is merely a problem of 
integration; for we know the original wave functions (they are the deter- 
minants mentioned above), and we can easily find the matrix components 
of energy with respect to them. The process is somewhat complicated, 
however, in many cases, and a careful arrangement of the arithmetic simpli- 
fies it greatly. To work out the actual method properly, we shall have to 
go more into detail than we have about the whole problem. This calculation 
will be made in sections 2 and 3; finally we give a number of examples of 
calculation of energy in section 4. , 

2. Our problem in the present section is this: to compute the diagonal 
term of the energy with respect to one of the approximate wave functions, 
written as a determinant, which we discussed earlier. If H represents the 
energy operator, and 


u(m/x1) -- - u(m,/ xy) 
re ees 
u(ny/ x1) - - - u(ny/ xy) 


is the wave function, we then have as the desired diagonal term 


f u* Hudv 
f u*udv 


The functions u are determinants; that is, sums of the NV! terms formed 
from the product u(m,/x,) - - - u(mw/xn) by carrying out all permutations of 
the coordinates x; - - - xv, even permutations having the coefficient +1, odd 
ones the coefficient —1. The product u*u is then a double sum, and each- 
of the integrals would be likewise; except that, in consequence of orthogon- 
ality relations, they really reduce to single sums, of which all the terms are 
identical. 

We consider first the energy integral. The energy operator H consists, 
as we shall see, of three parts: first, a constant, which we need consider 

















1306 J. C. SLATER 


no further; second, a sum of terms each depending on the coordinates of a 
single electron, and all the same in form, which we can write H(x)+ - - - 
+H (xy); third, a sum-of terms (the Coulomb repulsions) each depending 
on the coordinates of two electrons, which we can write explicitly as Z(all 
pairs) e*/r, where r is the distance between the two. The whole integral then ) 
consists of a constant, and two triple sums, two of the summations being over 
the permutations found in u and u*, the third over the terms in the energy 
operator. Consider a single term out of this triple sum, of the sort connected 
with the coordinates of just one electron. Then, if the permutations involved 
in wu and u*are different, at least two x’s must have different quantum numbers 
in the product; for example, we may have 


u*(n'/x,)u*(n"’/x)u(n'’/ xx) u(n’’”’/ x), 


seve 


where n’ is different from n’”’, and n” different from n’’’’. When we multiply 
this by a function of one electron, and integrate, there will always be at least 
one product (as u*(n'/x,)u(n’’’/x,)), which by orthogonality will integrate 
to zero. Thus the only terms which are not zero are those where we have 
the same permutation in uw and in u*; the double sum over permutations 
reduces to a single one. We are left, then, with 


(permutations of x - - - xx) f w(m/x) - ++ u*(ny/ xn) 


(H(x,)+--- +H (xg u(me1/ 21) - ++ u(ny/xy)dx, +--+ dxy. 


A single integral from this sum, by ‘normalization, immediately reduces to 
| f u*(m1/x1)H(x1)u(m/xi)dait--- + J 0/2 (ax) u(ne/ 29) day]. 


If we let these terms be J(m,), - - - J(mw), we see that the whole is simply the 

sum of the J’s, counted N! times;so that this part of the energy integral is 

N'\(I(m)+ ---+J(nw)), where it will be observed that the sum of terms 

H(x;), each depending on the coordinates of one electron, has been changed 

into a sum of terms each depending on a single set of quantum numbers. 

Next we take the terms e*/r;,; in the energy. As before, we take a typical 

a sien of the triple summation. If it happens that the same permutation occurs 
both in uw and u*, then the integral will reduce to 


Jit: w2f u*(n/x,)u*(n'/ x,)e2/rpu(n/ x.) u(n'/ x)dxdxi, 


which we will call J(m;7’). Each permutation will yield, from terms of this 
kind, the summation over all pairs x,x,; this reduces to the summation of 
J’s over all pairs n n’; and the N! permutations finally give, from terms of 
this sort, N!>° (pairs of n’s) J(n;n’). Next, it may be that the permutation 
in u* differs from that in u by the interchange of two electron coordinates. 
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If in addition the term e?/r;,,; happens to refer to the particular coordinates 
x,.x, which have been interchanged, we have an integral which reduces to 


Kor, 1") -f u*(n/x,)u*(n'/xp)e?/reu(n/xx)u(n’/x)dxxdx,, 


which we call K(n;n’). Such terms all occur with negative sign; for, since 
the permutation in u* differs by one interchange from that in uw, one of the 
terms is always an even, the other an odd, permutation, resulting in a coef- 
ficient —1in the double sum. Then by a similar argument to that above, 
these terms yield a contribution —N!>° (pairs of n’s) K (n ; n’) to the 
energy integral. Finally, if the term of u* differs from that of u by more 
than a single interchange, there will always be at least one coordinate whose 
wave functions will integrate to zero on account of orthogonality, so that 
no other terms exist. 

Finally we need the normalization integral fu*u dv. Formally, we can 
find this as we did the first part of the energy integral, replacing H(x,)+ ...+ 
H(xy) by unity. This results in replacing }-J(n) by unity; so that the integral 
is just N/. Thus division by the normalization factor simply removes the 
factor N/ which occurs in each term of the energy integral. We have then, 
for the diagonal term of the energy which we desire, constant +), (n’s) (mn) 
+)>°(pairs of n’s) J (n; n')— > (pairs of n’s)K(n; n’). 

3. Now we must consider the exact form of the wave function, to compute 
I,J and K. In connection with this, we must answer the question which no 
doubt will be felt at this point, as to the part which the spin plays in the cal- 
culations. We begin by considering the individual functions u(m;/x,). 
Each such function is a solution of Schrédinger’s equation for an electron 
in a central field. Let the potential of such a field be U(xyz). Then Schréd- 


inger’s equation is 
h? / 2 oe? 0? 
Aywu=| — —+—+ )+U)u=en. 
' ( Sxtm\ax? dy? ds? 





Here u is a function of x, y, z, and certain coordinates representing the spin, 
depending parametrically on m, /, m;,and m,; €is a function of the quantum 
numbers alone. To describe the spin, we proceed as Pauli!® does: we use‘as a 
coordinate the component of the spin along our fixed axis, m, (which*thus 
appears both as quantum number and as coordinate). Thus we could write 
our equation 


Hyu(n 1m, m,/x y zm,) =e(n lm, m,)u(n lm, m,/x y 2 m,). 


Since H, is independent of the spin coordinate m,, (neglecting magnetic 
interactions), we can separate variables, writing u as a product of a function 
of xyz, and a function of m,: 


u(n 1 m,m,/x y2m,)=u(n 1 m,/x y z)u(m,/m,). 


10 W. Pauli Jr., Zeits. f. Physik 43, 601 (1927). 
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The function u(m,/m,) is to be interpreted as follows: m, (the quantum 
number) can have two possible values, +4, and we have a different wave 
function for each quantum number. For example, we may have u(4/m,). 
But now assuming that the quantum number is 3, we know that the wave 
function is different from zero only if the spin points along the positive axis, 
so that m, (the coordinate) is 3; for m,= —}4, there is no wave function. 
Thus «u ($/3)=1, u(3}/—3)=0. We can then write u(3/m,)=65 (4/m,). 
Similarly u(—}/m,) =5(—3/m,); and we can write the whole symbolically 
as u(m,/m,) =5(m,/m,). We now have 


u(nlm,m,/x yzm,) =u(nl mi/x y 2)5(m,/m,). 
And for the first factor, the equation is 
Hyu(n 1 mi/x y 2) =e(nl)u(n 1 mi/x yz), 


where, since the energy does not depend on m, or m,, we have left out those 
quantum numbers in describing e. 

The solution of the central field problem is well known. We separate 
variables in spherical coordinates r 6¢. Then we have 


u(n 1 m,/x y s)=R(n1/r)O(l m,/6) &(m,/¢) 








where 
| (21+ 1)(1- ee 
Ol = P,\ ml 6 
(1 m,/@) | d+] ml)! 1!™!(cos 8) 
where 
1 di miit+i(—sin? 9)! 
P,'™!(cos 6) =—— sin! ™ 16 
21! d(cos 6)! mil+! 
and 


b(m, /p) =e'™'*/(2e)!”?. 
Thus we have 
u(n 1 m,m,/x y 3 m,) = R(n 1/r)O(1 m,/0) &(mi/o)5(m./m,) . 


To proceed further, we must investigate the value of H operating on the 
product of u’s. By definition, 


h? > e2 oe? e2 
H=- + ) V, 
82m él (— dy? + 02,2 + 








where V is the potential energy, given by 
N Ze? N e 
v=) (-=)+ Lii>j—> 
i=l vi 1 Vij 


where 7; is the distance of the 7th electron from the nucleus, 7;; the distance 
between the ith and jth electrons. Thus we have 


7 
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h? [ 9?u(11/ x1) 0?u(;/ x1) 


Hu(m;/ x) +--+ u(my/ xy) = — 
aia =~ mea 





82u(m,/ x1) 


7 Yul) - ++ u(ny/ aw) >>> +Vu(my/ x1) - + - u(ny/ xy). 


But by our assumptions, 


_ h? /0?u(m1/ 21) 4 tela 4 Balon/ si) 
8x2m\ Ox," dy,’ 02,” 





)- (€(m,) — U(x1))u(m/ x), 


etc. Thus we are left with 


Hu(n,/ x1) --- u(om/ax) =| > (e(n) — U(x) +V Jens - ++ u(ny/ xy). 


i=l 
This is of the form used above: Constant +)>H(x) +) -e?/r, where constant 
=)>(n)e(n), and H(x) = — U(x) —Ze?/r. 
We are now ready to compute the integrals J, J, and K. We must first 


note that by integrating over the coordinates of one electron, we really mean 
integrating over the dx dy dz, and summing over the spins: 


1/2 
[a= D (m) f dadyas. 
—l/2 


Thus, for example, we have the normalization and orthogonality of the 
individual wave functions: 


f u*(n'/x)u(n''/ x) 


1/2 


= Do(m,)8(m,'/m,)5(m,""/m,) 


1/2 
f u*(n'l’my'/ xyz)u(n''l''my"'/ xyz)dxdydz 


=65(m,'/m,")5(n'l'm!/n"'U'mz") . 


This orthogonality is all that is needed in the proofs of the preceding para- 
graphs. 
For the integral J, we have 


1/2 


I(nlmmm,) = >> (m,)5(m,/m,)5(m,/m,) 
-1/2 


f u*(nlm,/xyz)H(xyz)u(nlm,/xyz)dxdydz. 


The summation over m, merely reduces to the factor unity, independent 
of m,. When we insert the value of H(xyz), we note that the result is also 
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independent of m,; for H is a function of r only, the functions of angles 
integrate to unity, and we are left with 


I(nlmm,) = I(nl) = — J R*(nl/r)(U(r) +Ze?/r)4ar%dr. 


Next we are to find J and K. For J, we have 


J (nlm, ; n'l'mi'm,’) 
1/2 1/2 


= Di (ms) Do(ms)s 5 (m/ (me) x)5(m.'/ (ms) 1)6(m./(m,) x)5(m,"/(mz):) 
-1/2 -—l/2 
fee ota/ ose) omy) e*/ randy] sy asaa(n' mt / 915) 
dx yd y.dz,d xd y,d2;. 


The summation again reduces to unity, independent of m, and m,’. Then 
our quantity reduces to the integral, which we can write ' 


J (nlmym,; n'l'm,'m,’) =J (nlm; n'l'm/’) . 
Similarly we have 
K(nlmmm,; n'l'm,'m,’) 


1/2 1/2 


= Do (m,)x Do(m,)1 5(m,/(m,)1)5(m,'/(m,) x)4(m./(m,) 4)5(m,'/(m,),) 


-—l/2 -1/2 
f u*(nlm,/ xryiziu*(n'l’ my / xeyecne?/ree(nlmi/ xeynz,)u(n'l’my / xryizi) 


dxyd y.dz,d xd yid2). 
The summation reduces to 5(m,/m,'); so that we are left with 
K(nlmym,; n'l'm{'m,') =5(m,/m,')K(nlmi;n'l'm/’) , 


where the integral is symbolized by the last K. We thus observe that these 
exchange integrals only exist for electrons with spins parallel to each other. 
We are now to compute 
J(nlmi;n' l' m7 ) 
and 
K(nlmi;n'l' mi ) 
We recall the expression for u previously given. Also we use the familiar 
expansion 
1 (k—|m|)! r(a)* 


= Di(k,m 


r(xx") s (k+ | m | )! r(b) +} 








P,'™'(cos 0) P;'™!(cos 6’) exp (im(o—¢’)) 


where r(xx’) is the distance between (xyz) and (x’y’z’), and where r(a) is 
the smaller, r(d) the greater, of r and r’. Forming the expression for J, 
we have integrals of the form {f"exp(imd)dd which vanish unless m=O. 
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The double sum over k and m thus reduces to a single one over k, m being 
always zero. Then we easily have 


J(nlmy;n'l'm) = >(k)a*(lmi;l'my)F *(nl; n'1’) 








where 
“eatin = | mi |)! (2-+1)(0— | m,’ | )! 
e+ | mi)! (+ | mi|)! 
® ; 6 
J (eam(cos 0)|?P..°(cos )— dé 
0 
® ; ¢’ 
f [ Py! ™"(cos 0’) |? P.°(cos = 
and . 





F*(nl; n'l’) =e?(4x)? f € R?(nl/r) R2(n'l yt 


— ry tard’. 
r(b ro 
One notes that the a’s can be computed once for all, in terms of the spherical 
harmonics; the special properties of the atom in question appear only in 
the integrals F. 
In a similar way we form K. Here the integrals over ¢ are of the form 

*rexp i(m:—m,'+m)@ do, vanishing unless |m | = |mi:—mz,' |. Thus we 

have 
K(nlmijn'l'mi’) = >> (k)b*(Imi;3l'm/)G*(nl; n'l') , 

where 


(k— | mi—mi'|)! (20-+1)(0— | mi] )! (22’+1)(0'—| mi’ |)! 
(k-+ | mi—m,'| )! (1+ | m:| )! (l’+ | mi'| )! 





b* (imi; l'm,') = 


° sin 0 3 
| pf Paim(co 6) Py '™'\(cos 6) Py! ™—™" (cos a0 | 
0 
and 


k 
r?r'2drdr’ . 
+1 





G*(nl; n'l’) = e*(41)? J J “R(nl/r) R(n''/1) R(nl/r")R(n''/P’) * 


Using the values of the associated spherical harmonics one can compute 
the various a’s and 6’s; although the writer has not succeeded in setting up 
closed formulas for them, since this would involve the integrals of products 
of three spherical harmonics, an unfamiliar form." We give a table, includ- 
ing all the coefficients involved with s, p, d electrons (that is, /, /’ <2). 

We have now obtained the diagonal term of energy which we desired: 
in section 2 we have found it in terms of certain integrals J, J, K, and in 
section3 we have evaluated those integrals. Before passing to the examples, 
we should note one fact: that in finding the energy differences between 
multiplets, one needs only the integrals J and K, which do not depend 


11 See, however, J. A. Gaunt, I. c. Whether one has formulas or not, the table of values is 
certainly most convenient for computation. 
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TABLE OF a*(lm;; l’m,’). 
(Note: in cases with two + signs, the two can be combined in any of the four possible ways). 














Electrons l m, l’ my,’ k=0 2 4 
Ss 0 0 0 0 1 0 0 
sp 0 0 1 +1 1 0 0 
0 0 1 0 1 0 0 
pp 1 +1 1 +1 1 1/25 0 
1 +1 1 0 1 —2/25 0 
1 0 1 0 1 4/25 0 
sd 0 0 2 +2 1 0 0 
0 0 2 +1 1 0 0 
0 0 2 0 1 0 0 
pd 1 +1 2 +2 1 2/35 0 
1 +1 2 +1 1 —1/35 0 
1 +1 2 0 1 —2/35 0 
1 0 2 +2 1 —4/35 0 
1 0 2 +1 1 2/35 0 
1 0 2 0 1 4/35 0 
dd 2 +2 2 +2 1 4/49 1/441 
2 +2 2 +1 1 —2/49 — 4/441 
2 +2 2 0 1 —4/49 6/441 
2 +1 2 +1 1 1/49 16/441 
2 +1 2 0 1 2/49 — 24/441 
2 0 2 0 1 4/49 36/441 











taken together). 


TABLE OF D*(lm;; l’m,’) 
(Note: in cases where there are two + signs, the two upper, or the two lower, signs must be 














Electrons l m, yl mi =0 1 2 3 4 
ss 0 0 0 0 1 0 0 0 0 
sp 0 0 1 +1 0 1/3 0 0 0 
0 0 1 0 0 1/3 0 0 0 
pp 1 +1 1 +1 1 0 1/25 0 0 
1 +1 1 0 0 0 3/25 0 0 
1 +1 1 F1 0 0 6/25 0 0 
1 0 1 0 1 0 4/25 0 0 
sd 0 0 2 +2 0 0 1/5 0 0 
0 0 2 +1 0 0 1/5 0 0 
0 0 2 0 0 0 1/5 0 0 
pd 1 +1 2 +2 0 2/5 0 3/245 0 
1 +1 2 +1 0 1/5 0 9/245 0 
1 +1 2 0 0 1/15 0 18/245 0 
1 +1 2 Fil 0 0 0 30/245 0 
1 +1 2 F2 0 0 0 9/245 0 
1 0 2 +2 0 0 0 15/245 0 
1 0 2 +1 0 1/5 0 24/245 0 
1 0 2 0 0 4/15 0 27/245 0 
dd 2 +2 2 +2 1 0 4/49 0 1/441 
2 +2 2 +1 0 0 6/49 0 5/441 
2 +2 2 0 0 0 4/49 0 15/441 
2 +2 2 Fl 0 0 0 0 35/441 
2 +2 2 F2 0 0 0 0 70/441 
2- +1 2 +1 1 0 1/49 0 16/441 
2 +1 3 0 0 0 1/49 0 30/441 
2 +1 2 F1 0 0 6/49 0 40/441 
2 0 2 0 1 0 4/49 0 36/441 
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explicitly on the central field U at all. The integral J is needed only in finding 
the center of gravity of a multiplet system. The reason is that J depends 
only on m and 7, and so is the same for all the various degenerate states 
with which we start our perturbation problem. 

4. Examples. One electron outside closed shells. We need take but one 
case: a 3p electron outside completed K and L shells. Thus the scheme of 
electrons is (1s)?(2s)?(2p)§3p. We must now consider the various antisym- 
metric wave functions which are possible. By the method of symbolization 
mentioned above, in which we group in separate brackets the quantum num- 
bers of electrons with parallel and antiparallel spins, we see that there are 
six wave functions, which we give below, together with the values of >> mi 


and >> mM,: 


=m; =m, 
{ (100)(200)(211)(210)(21—1)(311)} {(100)(200)(211)(210)(21—1) 1 4 
{ (100) (200) (211)(210)(21—1)(310)} {(100)(200)(211)(210)(21—1) 0 4 
{ (100) (200)(211)(210)(21 — 1)(31 — 1)} { (100)(200)(211)(210)(21—1)} —1 4 
{ (100)(200)(211)(210)(21—1)} {(100)(200)(211)(210)(21—1)(311)} 1 ans 
{ (100) (200) (211)(210)(21—1)} {(100)(200)(211)(210)(21—1)(310)} 0 —4 
{ (100) (200) (211)(210)(21—1)} { (100)(200)(211)(210)(21—1)(31—1)} -1 -—3 


We note that the arrangement of >> m:, >. m,, is just that for a single multiplet 
*P,. There are no cases in which more than one term has a given value of 
>>mi, >-m,; thus there is no need of applying the sum rule at all. The diagonal 
terms of the energy, computed with respect to these six wave functions, 
should give directly the energies of the six terms of the multiplet. But now 
.we come back to our general principle; these six terms must have the same 
energy. We must actually compute the energies by our rules, and see that 
they are the same in each of the six functions. As has been mentioned be- 
fore, we need only use the J and K terms. These terms, we recall, were 
>> (pairs) J —>>(pairs with par. spins)K. Now in these sums, many terms are 
the same for each of the six wave functions; all the terms, in fact, relating 
to pairs of electrons both in the closed shells. These terms can exert no 
influence on the multiplet separations, or anything of that sort. Thus we 
can leave them out, for our present purpose, as we left out the terms depend- 
ing on the integrals I. The only terms we need retain are those in which our 
3p electron is a member of the pairs. Thus for the first wave function we must 
compute the following: 


2J(311; 100) -+27(311; 200) -+2/7(311; 211) -+2/(311; 210) +27(311; 21—1) 
— K(311; 100) — K(311; 200) — K(311; 211) —K(311; 210) — K(311; 21—1). 
For the second and third, we substitute respectively 310, 31-1 in place of 311; 


the fourth, fifth, and sixth evidently give the same three results already 


given. Let us new group these terms according to the shells that the 3p’s 
interact with: 


1s:2J(311; 100) — K(311; 100) 

2s:2J(311; 200) — K(311; 200) 

2p:2F(311; 211) +2/(311; 210) +27(311; 21—1) — K(311; 211) 
— K(311; 210) —K(311; 21-1). 
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For the interaction with 1s, we have 


2J(311;100) =2 >-(k)a*(11; 00)F*(31; 10) = 2F(31; 10) 
— K(311; 100) = — 5“(k)b*(11; 00)G*(31; 10) = —4G"(31; 10). 


Before going further, let us find the corresponding terms in the interaction 
of a 310 electron with the 1s shell: 


2J(310; 100) =2 5>(k)a*(10; 00)F *(31; 10) = 2F°(31; 10) 
— K(310; 100) = — }°(k)b*(10; 00)G*(31; 10) = —3G1(31; 10). 
That is, the interaction integrals of an outer p electron with an s shell are 
the same whether the p electron has m,;=1or 0 (or—1, as one immediately 


verifies). We shall show the same result to hold for the interaction with 
the 2p shell: for the (311) electron, the terms are 


27(311; 211) +2/(311; 210) +2(311; 21—1) — K(311; 211) 
— K(311; 210) —K(311; 21—1) 
=2 S\(k)a*(11; 11)F*(31; 21) +2 Do(k)a*(11; 10)F*(31; 21) 
+2 >0(k)a*(11; 1—1)F*(31; 21) — 55(k)b*(11; 11)G*(31; 21) 
— >o(e)b*(11; 10)G*(31; 21) — 55 (k)b*(11; 1—1)G*(31; 21) 
=2 5>0(k)(a*(11; 11) +a*(11; 10) +a*(11; 1—1)F*(31; 21) 
_ >o(k)(b*(11; 11)+6*(11; 10)+5*(11; 1—1))G*(31; 21) 
= 6F (31; 21) —G°(31; 21) — 3G°(31; 21). 
Similarly for the (310) electron interacting with the 2 shell, we have 
2 >°(k)(a*(10; 11) +a*(10; 10) +a*(10; 1—1)F*(31; 21) 
— 0(k)(b*(10; 11) +5*(10; 10) +5*(10; 1—1))G*(31; 21) 
= 6F (31; 21) —G°(31; 21) —2G2(31; 21). 
This agrees with the former value, showing that the interaction of either 
a 311 or a 310 (or, by a simple extension, a 31-1) electron with a completed 
p shell, give the same result, independent of m;. Putting all these results 
together, all six levels of the *P multiplet give the same energy, by direct 
computation. Of course, this is merely a check of our general theorem that 
all the levels of any multiplet must have the same energy. The special 
properties of the a’s and b’s which lead to this result could be proved by use 
of that theorem. 

Any configuration outside closed shells. We have just seen that the inter- 
action energy of a single electron with a closed shell is independent of the 
m, of the outer electron. We have proved this by direct computation for a 
p electron interacting with s or p shells, but we could extend the result to 
the general case, from our general theorem that all the terms of any multiplet 


have the same energy. But this has an important bearing on our general 
problem. For in any case an atom consists of a certain number of electrons 
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outside closed shells, and the central closed shells. In the various unper- 
turbed wave functions which we are to use, the outer electrons have different 
m,’s. In the energies of the terms, we are to compute all J’s and K’s con- 
nected with pairs of electrons in the atom. We have already seen that all 
pairs, both of which are in closed shells, will give identical contributions 
to each of the unperturbed terms. But now we can go further: all sets of 
pairs, in which one is an outer electron, the other one of the electrons of 
a closed shell, will, when summed over the electrons of the closed shell, 
give the same result for each unperturbed state. That is, asfar as multiplet 
separations are concerned, closed shells exert no influence at all; they affect 
only the position of the whole set of multiplet terms. It is well known that 
the classification of the terms is independent of the existence of closed shells; 
this proves that the energy relations also depend only on the outer electrons.” 
For the rest of our examples, then, we shall consider only those electrons 
which are outside closed shells. 

Two electrons outside closed shells, one in s state (helium, alkaline earths). 
Suppose we have the scheme (ms), (m’p), for example. There are, in this 
case, the unperturbed wave functions symbolized by 


=m) =m, 
{(m00)(m’11)} { } 1 1 
{(m00)(n’'10)} | } 0 1 
{(n00)(m’1—1)} { } —1 1 
{(mn00)} {(n’11)} 1 0 
{(n00)} {(n’10)} 0 0 
{(n00)} {(n’ 1 —1)} —1 0 
{(m’ 11)} {(n00)} 1 0 
{(m’ 10)} {(n00)} 0 0 
{(n’ 1 —1)} {(n00){ -1 0 
{ }{(n00)(n’11)} 1 —1 
{ ${(n00)(n’10)} 0 —1 
{ }{(n00)(m’1 —1)} —1 ae | 


By our general scheme of classification, we have a 'P and *P term. The 

terms with dom, dom, equal respectively to (11), (01), (—11), (1 —1), 

(Q —1), (—1 —1), belong to the *P state. On the other hand, the remaining 

terms, as (1 0), are each degenerate. The sum of the energies of two such 

terms equals the sum of the 'P and *P energy. Thus for example we have 
3P:{(m00)(n’ 11)} 


‘P+'P:{(n00)}{(m' 11)}+{(m' 1 1)} {(n00)} =2{ (m0 0)} {(m’ 11}. 
That is, for the energies, we have 
3P:J(n00;n’ 11)—K(n 00; n'1 1) 
8P+!1P=2)(n00;n’ 11). 
Therefore 'P=J(n 00; mn’ 11)+K(n 00; n’ 1 1) 
The singlet and triplet are thus given by a definite value + K(m 0 0;n’ 1 1). 


This checks with Heisenberg’s™ calculation of this case; our integral K is 
readily seen to be the same exchange integral which he computes. 


12 See W. Heitler, Zeits. f. Physik 46, 70 (1928). 
13 W. Heisenberg, Zeits. f. Physik 39, 499 (1926). 
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Shell of equivalent p electrons. (Elements C,N,O,F, etc.) »*. This is the 
case shown in Fig. 1c. There are 15 wave functions, of which we give those 
with ).m,=0, }.m,= 0; for simplicity we omit the total quantum number from 
our descriptions, so that (1 0), for example, stands for (m 10). Then, cor- 
relating the >-m:, >.m, with the multiplets represented (from the figure), 
we have 











Multiplet Mi M, Wave function 
D 2 0 {a1} 10) 
ID+3P 0 {(10)} {(1 1)} 

1 0 {(11)}  {(10)} 
ID+%P+1$ 0 0 {(11)} {(1 —1)} 

0 0 {(10)} {(10)} 

0 0 {(1 —1)}{ (1 1)} 
- 1 i (11) (10) f 4 
- 0 1 (14a —y) t } 








Now we have for the energies 


1 
1D:J(n11;"11)=F%(n1; n +r 1;n 1) 
2 3 
8P:J(n11;n10)—K(n11;"10)=F%(n1;n - 1;n j- se 1; 1) 


We note that, for equivalent electrons, the F’s and G’s of the same indices 
are equal. Hence 


5 
*P=F%(n 13m 1)——F%(n 13m 1). 


We can check the same value from the other *P term (2m,=0, =m,=1) 
and from the two terms giving 'D+*P. Finally we have 


1D+3P+'S:J(n11;”1—1)+J(n10;"10)4+J(n 1 —1;3;”1 1) 
6 
=3F%(n1i;n +e 1;” 1). 


10 
Therefore 'S=F°(n1;n += F2(n 1; 1). 


We note that the term F°(m 1; 1) is common to all the levels; thus we 
can leave it out, as we have all the terms which do not affect the separation. 
We note from its definition that F?(m 1; 2 1) is positive. Thus we see that 
of the three multiplets, the *P lies lowest, 'D next, 'S highest, in accordance 
with Hund’s rule that the terms of largest L and S lie lowest. The separa- 
tions are in a simple ratio: *P—'!D=6/25 F*(m 1; 1), and 'D—'!S=9/25 
F*(n 1; n 1) so that the ratio is 2 to 3. 

We should find an example in the lowest levels of C, (1s)? (2s)? (2p)?; 
these levels, however, have not been observed, as far as the writer knows. 
For Si (1s)? (2s)? (2p)® (3s)? (3p)?, however, the term values are *P = 65,615, 


™ McLennan and Shaver, Roy. Soc. Canada 18, 1, (1924). A. Fowler, Proc. Roy. Soc. 
A123, 422 (1929). 
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1D =59,466, 'S=50,370, giving *P —'D = 6149, '1D—1S= 9096, the first giving 
F*(3 1; 3 1) =25,602, the second F*(3 1; 3 1) =25,267, in good agreement. 

p*. The method is so similar to that used for p* that the calculations need 
not be given in detail. The multiplets are ‘4S *D ?P, the *S lying lowest and 
*P highest, again in agreement with Hund’s rule. The separations are given 
by 4‘S—*D=9/25 F°(n1;”1) *D—*P=6/25 F°(n 1; 1). An example is 
found in the normal N spectrum, (1s)? (2s)? (2), observed by Compton and 
Boyce.™ They find *S=117,345, *D=98,143, *P=88,537,giving *S—*D 
=19,200, ?7D—?P=9,600. From the first, F*(2 1; 2 1) =53,400, and from 
the second, it is 40,000; a somewhat poorer agreement than before, probably 
on account of the tighter binding of the electrons. 

p*. Here the relations prove to be as in ~* as regards energy, as well as 
in the arrangement of terms: *P—'D=6/25 F*(n 1; 1) ‘'D—!S=9/25 
F*(n 1; m 1). An example is normal O, (1s)? (2s)? (2p)*. Here it is estimated"® 
that, counting terms up from *P as zero, '!D =25,500, and 'S=65,000, giving 
§P—1D =25,500, 'D—!S=39,500. From the first, F?(2 1; 2 1) = 106,000, and 
from the second 109,700, a very good agreement. It is interesting to note 
the increase in the integral from N to O on account of the tighter binding. 

p' and p® yield each only one multiplet, so that they need not be con- 
sidered. 

Shell of equivalent d electrons (Iron group) d*. The multiplets are *F 
3P 1G'!D 1S. When we work out the separations, however, it appears that 
they are not arranged in this order; the singlet terms are anomalous, dis- 
obeying Hund’s rule, in that the 'D lies lower than 'G. Except for this, 
however, the arrangement is as we should expect, *F lying below *P, and be- 
ing the lowest term of the combination. The separations are given by 


3F—%P = (135/441) F°(m 2; m 2) —(75/441) F*(n 2; n 2) 
§F—1G = (108/441) F?(m 2; m 2)+(10/441) F*(n 2; n 2) 
$F—1D=(45/441) F?(nm 2; m 2)+(45/441) F*(m 2; m 2) 
3F—1S= (198/441) F?(m 2; m 2)+(135/441) F*(m 2; m 2) 


The separations now depend on the two parameters F°*(n 2; 2) and 
F*(n 2; m 2), which bear no fixed relation to each other. Nevertheless we 
can estimate their relative magnitude. For by definition F* is the integral 
of a certain function of 7; and re, multiplied by r.*/r,*+! where 7 is the greater, 
r, the less, of r; and re. Thus increasing k necessarily decreases the integrand, 
and hence the functionF: F*(n 2; n 2) <F*(n2 ; n 2). Rough calculation in- 
dicates that the decrease is about to a half. If then we provisionally take F* 
to be half of F?, we have the separations in the ratio , 


3F —3P:135—37= 98 3F—'D: 45+22= 67 
*F—'G:108+ 5=113 3F —'$:198+67 = 265. 
The order of terms is thus expected to be *F 'D *P 'G'S. 


4K, T. Compton and J. C. Boyce, Phys. Rev. 33, 145 (1929). 
4 McLennan, McLeod, and Ruedy, Phil. Mag. Sept., 1928, p. 558. These values are esti- 
mated from the energy level diagram in that paper. 
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Experimentally we find an example in the normal spectrum of Ti,!’ 
(1s)? (2s)? (2p)® (3s)? (3p)® (4s)? (3d)*. Here the order of terms is in fact just 
what our calculation predicts, the exception to Hund’s rule being found 
experimentally. The observed separations are approximately 


3F—8*P= 8500 5F—'D= 7200 
5F —'G=12100 5F —'§=15100. 


If we assume that F* is really half of F?, these give F?(m 2; m 2) equal res- 
pectively to 38,000, 47,000, 41,000, 72,000, in fair agreement except for the 
last one. This last results from the 'S, which is not nearly so far above the 
rest of the terms as the theory would indicate—probably because the second 
order corrections for this term would be large, and would have the effect 
of depressing it. The agreement with observations can be somewhat, but 
not much, improved by taking a slightly different ratio of F* to F*. 

d’..d*°. We shall give only the lowest terms for the rest of the d shell. 
The terms of highest multiplicity prove to lie in general lower than the others, 
as the rule would predict. Of these, we have for d*: ‘F *P; d‘*,5D; d5, °S; 
d°, 5D; d'’, *F 4P; d’, *F*P. Evidently the only cases where there are signifi- 
cant separations in the multiplets of highest multiplicity are the d? d* d’ d® 
F-P separations. When one calculates, one discovers the fact that all these 
are given by the same formula, (135/441) F*(m 2; m 2) —(75/441) F'(n 2; 2). 
This permits an interesting comparison with experiment: we can compare the 
observed separations for Ti (d?), V(d*), Co(d’) and Ni (d*). We should 
expect these to increase regularly with the number of d electrons. Experi- 
mentally this separation is about 8000 for Ti, 9500 for V, 14000 for Co, 
15000 for Ni, indicating a fairly uniform increase of about 1000 to 1500 for 
the addition of one d electron. 

Non-equivalent p’s: Two p’s: pp’. The multiplets are *D *P *S'D'P 'S. 
The triplets are particularly interesting, and we give their separations. If 
the principal quantum numbers of the electrons are m and n’, they are 


3) —3P = —6/25F*(n 1; n’ 1) +2G%(n 1; n’ 1)—4/25G%(n 1; n' 1) 
3) —8§ =9/25F2(n 1; n’ 1) —9/25G2(n 1; n’ 1). 


The significant feature of these results is that they depend in such a com- 
plicated way on several integrals. That is, in a case like this— and it is 
the simplest set of multiplets that can be built up from non-equivalent elec- 
trons—we must not look for simple numerical relations between the separa- 
tions, or even for a definite, fixed order for the terms. We may rather expect 
that, as we go from one element to another, the relative order of terms 
can change. 

We find an example in Nt (1s)? (2s)? 2p 3p.'8 Here the experimental 
separations are approximately *D—*P=4000 *D—*S=2300, disobeying 
Hund’s rule, in that *S lies below *P. We can make from these figures a rough 
estimate of the integrals F and G. We note that G differs from F in having 


17H. N. Russell, Astroph. Jour. 66, 283 (1927). 
18 A. Fowler, Proc. Roy. Soc. A107, 31 (1925). 
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the product R(n/x:) R(n'/x,) in its integrand instead of R*(n/x:). The first 
is less than the second—very much less if the orbits m and m’ are of decidedly 
different size, for then either R(n/x,) or R(m’/x;) will be small through most 
of the range of x;. Thus each G integral is small compared to the corres- 
ponding F. It is reasonable to suppose that G° is of the same order of mag- 
nitude as F?. Further, G? will be smaller than G°, by analogy with what we 
have already seen, and as a rough assumption we may take G?=} G*®. Thus 
we have the two equations 


4000 = — (6/25)F?+ (48/25)G° 
2300 = (9/25)F2—(9/50)G°, 


giving F’(21;31)=7900, G°(21;3 1) =3070, G’(21;3 1) =1535. 
These values are reasonable, but provide no definite check for the equations. 
Five equivalent p's, one other p: p'p’. Again we have *D*P*S 'D 

1P 1S, by combination of the p’ with the ?P of p5. But now it proves that 
the triplet separations are quite different, and have a much simpler formula. 
We have 

3D) —*P =(6/25)F2(n 1; n’ 1) 

8) —38§ = —(9/25)F2(n 1; n’ 1) 


depending on only one parameter (the others all cancel out, seemingly, al- 
though perhaps not really, by chance). Thus we can predict definitely the 
order of terms: *S lies lowest, then *D, finally *P, in direct contradiction to 
Hund’s rule. We find an example, however, in Ne (1s)* (2s)? (2p)° 3p, 
Paschen’s p terms," and it definitely verifies the contradiction. Experiment- 
ally, Paschen’s 2:0, which is the *S, lies well below any others. On account 
of the wide multiplet separations it is impossible to show any good check 
of the ratio 2:3 for the 3D —*P and *D—®8S separations. If one takes centers 
of gravity of Paschen’s terms, one finds 


3D—*P=570 3D —*S = — 2023 


the first giving F?(2 1; 3 1)=2370, the second 5620. The correct figure 
is probably between these. It is interesting to note the agreement as to order 
of magnitude of this with the same one, 7900, for N*. 


NOTES 


1. We wish to prove that the energy has no matrix components (taken 
with reference to the incorrect, approximate wave functions) connected with 
transitions in which >.m, or >.m, changes. The essential part of the proof 
is the demonstration that the energy operator H commutes with the operators 
connected with }-m,or).m,. For)>_m, the proof is obvious:>_m, is a quantity 
depending only on the spins, H only on the coordinates, and operators de- 
pending on entirely independent quantities always commute. For > mz, the 
essential point is that }\m, is the operator connected with an infinitesimal 
rotation of space about the z axis, and on account of the fact that H is in- 


°F, Paschen, Ann. d. Physik 60, 405 (1920). 
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dependent of orientation, such a rotation does not affect it, and so commutes 
with it. To be more specific, one can write the three components of angular 
momentum of an electron in operator form, using the operators for linear 
momentum. Put in spherical coordinates, the z component of angular momen- 
tum is represented by the operator (h/27i)0/0¢. Operating on a single elec- 
tron wave function with the factor e**, this operator reduces to the multi- 
plication by m,h/2z, so that it has a diagonal matrix, and m; measures the 
z component of angular momentum, in units h/27. Then the operator con- 
nected with }>m; is —i)_0/0¢, where the sum is over the ¢’s of the various 
electrons in the many-electron wave function; this operator is in the many 
electron problem again a diagonal matrix. This holds, we note, even with 
the incorrect, unperturbed wave functions. Except for the factor —i, this 
operator simply represents the change in the function it operates on, if all 
¢’s are increased by the same amount; that is, if the whole electronic system 
is rotated rigidly. Then we have 


0 oH re) 
(- i Do) a =—i) —u-iH )—u. 
dp dg dg 
But on account of its spherical symmetry, 20H/d¢=0. Hence 


(- De) au u(- i r-) “u, or 
( dom) H—H( dom, =0, 


showing that the operator >.m; commutes with H. 

Knowing that both }>m; and >-m, commute with H, and that both have 
diagonal matrices, the rest of the proof is simple. We merely write down the 
commutation laws in matrix form: 


> (n") >omi(n'/n")H(n"/n'") —H(n'/n'’) domi(n''/n'") =(). 
On account of the diagonal relation, this amounts to 
( Som(n'/n’)— Somi(n'"/n'"))H(n'/n') =0. 


This cannot be satisfied unless either }>m,(n’/n’) =)m,(n'"’/n’’’) or unless 
H(n'/n’'’")=0; that is, H has matrix components different from zero only 
if }-m, has the same values in initial and final states. The same proof holds 
for >om,. 

2. Theorems regarding angular momentum. It is assumed that we have 
wave functions which, although not exact, are the good approximations 
used in this paper; that is, linear combinations have been made so that the 
energy matrix has components different from zero only between states of 
decidedly different energy; it is diagonal as far as the states are concerned 
which come from one set of electron quantum numbers. We wish to show 
first that matrix components of angular momentum have non-negligible com- 
ponents only between states of the same energy (or diagonal term of the 
energy matrix). Suppose we let the operators connected with the three com- 
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ponents of orbital angular momentum be M,*, Mz", Mx1*, where Mz’ is the 
quantity usually denoted by Mz, and is equal to }.m;. Similarly for spin 
we have Ms*, Ms’, Ms*. Now we have seen that M,;* and Ms* commute 
with H; and, since there is nothing peculiar about the z axis, we can equally 
well show that M,*, Mz", Ms*,Ms" commute with H. Thus for example 
Ms*H—H Ms*=0, or in matrix form 


>> (n")Ms7(n'/n")H(n"/n'") aan H(n'/n'’)Ms7(n"/n'" ) =(. 
We can write this, for non-diagonal terms, 
Msg?(n'/n'") [H(n'"’/n'"’) oe H(n'/n’) ] 
_——— do (n"# n’ , n'’’) [Ms?(n'/n")H(n"/n'") _ H(n'/n’’) Mg?(n"/n'"’) ] : 


The non-diagonal terms of the energy, H(n’’/n’’’), etc., are by hypothesis 
different from zero only if m’’ and n’’’ refer to states with different 
electron quantum numbers, and even so they are small of the first order. 
If the right side were precisely zero, we should have obviously the result that 
M57(n'/n''’) wasdifferent from zeroonly if H(n’’’/n’’’) was equal to H(n’/n’) 
which is what we wished to show. As it is, this result is true only to the first 
order of small quantities: Ms* can have other, small components, and 
H(n'''/n''’) can differ from H(n'/n’) by small quantities. Thus we should 
expect that our approximate wave functions for a multiplet would give 
only approximately equal energies. This is not true, however, when we actu- 
ally calculate by our method; the energies are precisely equal. The reasén 
can be easily seen. Let the states m’ and n’”’ be both connected with the 
same electron quantum number. Then the matrix components on the left 
side of the equation above are all found in terms of integrals over wave 
functions of those particular electron quantum numbers. The matrices on 
the right, however, involve other states, m’’, which must refer to different 
electron quantum numbers, so that the right side would involve different, 
independent integrals. The two sides could not be in general equal unless 
each was zero. Hence the result is: if the approximate wave functions are 
computed by our method, the angular momentum has components, be- 
tween twostatesof thesameelectron quantum numbers, only if both states have 
the same energy. Itscomponents between statesof different electron quantum 
number are small of the first order. From now on, we can neglect the latter 
terms, as we did with the energy. Thus the matrix of the angular momentum 
is of the form shown in Fig. 5. This can be compared with Fig. 3, showing 
the matrix of the energy in the final wave functions; but Fig. 3 should now 
be numbered to show that the diagonal terms of energy within one of the 
small squares of Fig. 5 (that is, in a multiplet) are all equal. 

We have shown that the wave functions connected with a set of m’s and 
l’s are divided into groups, each wave function of a group having the same 
energy value, and the angular momentum having components only between 
different wave functions of a group, not from one group to another. We 
wish now to show that the matrix components of angular momentum within 
such a group are really such as to indicate a vector of magnitude L, and 
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another S. To do this, we first note that M,*, Mz", Mz', Ms*, Ms*, Ms? 
satisfy the commutation relations 


h 
M,YM1z'-—-M1'My,= ~—M 1", otc., 
1 


2x 


h 
Ms’M3*—M3'Ms"= “3” etc., 


Tl 
M1? Ms*— Ms* M,7=0, etc. 


We show these by proving them for the separate m,’s and m,’s of the separate 
electrons, and then combining by the relations M1*=>om:, etc. For the 
separate electrons, the results for the m,’s are well known; for the m,’s, 
the results come directly from Pauli’s theory of the electron spin. We also 
note; that M,* and Ms‘ have diagonal matrices, with diagonal values Mz 
and Ms (in units of h/27). Now it can be proved®® that if we have aset of 
matrices related by the commutation rules M,"M;?—Mi? M,Y¥=— 
(h/2mri)M 17, etc.; if these have components only between a limited set of 
statesand if M,* forms a diagonal matrix; then the whole matrix is uniquely 
determined, except for a phase constant; further, M ,* must have integral, 
or half integral, characteristic values M_, ranging from a value—L to L; the 
states can be described by the values of Mz, so that can be taken as a quantum 
number. Then the matrix components prove to be given by 


~ 1/h\? 
M1*(M1/(M1i+1))M17((M1i+1)/M1) --(5) [L(L+1)—Mi(Mz+1)]--- 
M1?(M,/M1)=M_z sat 
(M17)? +(M1")?+ (M1)? |(M1/M1') = () L(L+1)8(M1/M1z’). 


Such a set of matrix components is the unique description of a vector L, 
in its (2L+1) possible orientations. Similar results apply to Ms. 

In our case, these relations lead immediately to the following results. 
In the first place, the M,?’s depend only on orbital coordinates, the Ms*’s 
only on spins. Thus the M,?’s can have components only between two 
wave functions in which the spin appears in the same way (and which thus 
have the same Ms). Similarly Ms*’s can have components only between 
states of the same M,_. Therefore for each Ms there are a number of states of 
the multiplet, with different Mz ,’s and having matrix components of the 
Mz*’s between them as given above. These different sets of states, on the 
other hand, can differ only in the spin function, so that the components of 
M*’s are the same for each set witheach Ms. Similarly, for each Mz, there 
are a number of states, with different spin functions but with the same orbital 
functions, different Ms’s, and components of Ms?’s just as above. The 
result is an array of terms which can be plotted in precisely the rectangular 
form used in the paper, and corresponding to the (2Z+1) orientations of 
the vector L, combined with the (2S+1) orientations of S. 


2° Born, Heisenberg, and Jordan, Zeits. f. Physik 35, 557 (1926). 








NOVEMBER 15, 1929 PHYSICAL REVIEW VOLUME 34 


PHENOMENA IN OXIDE COATED FILAMENTS 


By Josern A. BECKER 
BELL TELEPHONE LABORATORIES, 463 WEsT STREET, NEw YorK, N. Y. 


(Received October 4, 1929) 


ABSTRACT 


A theory of the changes in activity in oxide coated filaments is proposed. From 
a comparison of the behavior of these filaments and filaments with composite surfaces 
such as thorium on tungsten, caesiumon tungsten, and caesium on oxygen on tungsten it 
appears probable that oxide coated filaments owe their high activity to adsorbed metal- 
lic barium. The changes in emission from a coated filament produced by changes in 
plate potential and by currents sent into or drawn from it, are ascribed to electrolysis 
of the oxide. When electrons are sent into a coated filament barium is deposited on the 
surface and the activity increases until an optimum is reached beyond which the acti- 
vity decreases. When current is drawn from the oxide, oxygen is deposited on the sur- 
face. If the oxygen is beneath the adsorbed barium, it increases the activity; if it is 
above the barium, it decreases the activity. Both barium and oxygen diffuse readily 
from the surface into the oxide and vice versa. This theory is tested, confirmed, and 
extended by numerous experiments. 

An experimental technique is employed by which relative rates of evaporation of 
small amounts of electropositive and electronegative materials can be determined with 
considerable precision. The same technique might be useful in a number of similar in- 
vestigations. Metallic barium or oxygen which evaporate from a coated filament are 
allowed to deposit on one side of a flat tungsten ribbon whose thermionic activity is fol- 
lowed. When the plausible assumption is made that an optimum activity is obtained 
when the tungsten is covered with a single layer of electropositive material, the relative 
rates of evaporation can be converted to absolute rates. This technique is also em- 
ployed to determine the factors which control the evaporation of oxygen from a coated 
filament. 

Experimental results.—(1) Metallic barium evaporates from a well broken 
down coated filament provided that the temperature is not too high and that no ap- 
{ preciable currents are drawn. (2) Oxygen evaporates if an emission limited current 
is drawn and the temperature is high enough. (3) Metallic barium which is deposited 
on the surface of the oxide either by electrolysis or by evaporation from an outside 
source, produces changes in activity of the same sort as are produced by barium on 
tungsten. (4) Barium or oxygen diffuse readily to and from the surface of the oxide 
(5) Oxygen brought to the surface of the oxide by drawing current from it, increases 
the activity under some circumstances and decreases it under others. (6) Most of the cur- 
rent through the oxide is carried by electrons but a small part is carried by ions; the ratio 
depends upon the temperature, the composition of the oxide, and the plate potential. 
It appears that the conductivity of the oxide does not obey Ohm's law and that some of the 
conductivity is due to a diffusion current which may flow against an opposing field. 
(7) The complex time changes that occur in coated filaments are described ina system- 
atic manner. 


INTRODUCTION 
METALLIC filament which is covered with one or more of the oxides 
of calcium, strontium, or barium is called an oxide coated filament. 


It has long been known that such filaments when properly treated are ex- 
ceedingly good emittors of electrons. They have found considerable use in 
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thermionic vacuum tubes and in scientific investigations. The object of 
this paper is to present a scientific description of some of the processes oc- 
curring in some oxide coated filaments and to propose a comprehensive 
theory to account for these processes. 

There are two types of oxide coated filaments: those in which some or 
all of the oxides have been made to combine chemically with the metallic 
core of the filament, and those in which the oxides are uncombined. The 
latter appear pure white. In the combined filaments, the chemical compounds 
are decomposed when the tube is pumped. This leaves a finely dispersed 
metal mixed with the oxides and gives the filaments a grayish appearance. 
A method of preparation of combined filaments is described by Arnold.! 
The preparation of uncombined filaments is described by Schottky, Rothe 
and Simon, Volume 13, Part 2 of the Wien-Harms “Handbuch der Experi- 
mental-physik” pp. 215 and 305.2 The reader will find a good discussion of 
published articles on oxide coated filaments in this same volume. The au- 
thor’s experience with barium oxide+strontium oxide combined filaments 
has been more extensive than with others. 

In either type of filament, but particularly in combined filaments, it is 
found that the thermionic activity is greatly enhanced by what has been 
variously termed the activation, formation, or breakdown process. We shall 
call it the breakdown process. This process consists in glowing the oxide 
coated platinum filament at about 1500°K for several minutes; a rather high 
positive potential gradient is then applied to the filament for something like 
two to thirty minutes. During this time the emission from the filament in- 
creases greatly. When it attains rather large values such as 0.1 to 1.0 amps. 
per cm’, the applied field and the filament temperature are so adjusted that 
the current is just limited by space charge. The detailed procedure to be fol- 
lowed for any particular filament constitutes an art and is beyond the scope 
of this article. The processes that go on during this time will be discussed 
below. The behavior of the filament before and after breakdown is con- 
siderably different. 

Throughout this article we shall frequently use the term “activity” 
of the filament. As a measure of the activity we shall take the emission cur- 
rent at one particular temperature or at a series of temperatures. These 
temperatures will all be low and the plate potential will have a small constant 
value in order that the currents may remain essentially constant during the 
time required to determine the activity. 

One of the outstanding characteristics of a coated filament is its complex 
behavior. For example, coated filaments show changes in activity with time 
following upon or accompanying (1) changes in temperature, (2) changes in 
plate potential and (3) current sent into or drawn from the filament. While 
all coated filaments show time changes with these three variables, the magni- 
tude and even the direction of the change in activity depends greatly on the 
chemical and physical composition of the oxide as well as on the previous 


1H. D. Arnold, Phys. Rev. 16, 70 (1920). 
2 The title to this volume is “Gliihelektroden und Technische Elektronenréhren ” 
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treatment both immediate and remote. From this it is apparent that the 
coated filament is an extremely complex system. In our attempt to discover 
the processes going on in this complex system, we have found it convenient 
to be guided by a tentative hypothesis or unifying picture, and the reader 
may find it helpful to consider this hypothesis. 


HYPOTHESIS OR THEORY OF THE ACTION IN COATED FILAMENTS 


In constructing such a hypothesis it is desirable to consider simpler 
types of filaments and to compare them with oxide coated filaments. The 
simplest kind of filament is one having an atomically clean and homogeneous 
surface such as clean tungsten, tantalum, or molybdenum. For any of these 
filaments the emission is determined solely by the temperature and plate 
potential; it is independent of the time of treatment at any previous temp- 
eratures, plate potentials or emission currents. 

A somewhat more complex type of filament is illustrated by thoriated 
tungsten,’ caesium on tungsten,‘ or barium on tungsten. These filaments 
have composite surfaces: the tungsten surface is more or less covered by a 
partial layer, a complete layer, or more than one atom layer of a foreign 
element. The activity, or emission at a testing temperature and fixed plate 
potential, now depends on the kind of foreign element and on 8, the fraction 
of the surface covered with it. Figure 1 shows the logarithm of the emission 
for barium on tungsten at 1100°K as a function of #6. Note that the activity 
is increased by a factor of about 10" by a complete layer of barium. Since 
the amount of barium on the surface can be changed by heating the filament, 
it follows that in this case the activity can be greatly modified by a tempera- 
ture treatment. On the other hand, experiment shows that the activity and 
6 are not affected by the plate potential or by the electron current sent into 
or drawn from the filament. 

A still more complex filament is one with a doubly composite surface such 
as caesium and oxygen on tungsten. In this case experiment apparently shows 
that the activity depends upon (1) the number of caesium atoms per cm’, 
(2) the number of oxygen atoms per cm?, and (3) the relative position of the 
caesium and the oxygen. The activity increases with the amount of caesium 
until the surface is covered with a monatomic caesium layer; then it decreases 
as the second layer is formed. If we start with a monatomic layer of caesium 
and add increasing amounts of oxygen to the hot filament, the activity is 
increased even more, provided the amount of oxygen is small; as more and 
more oxygen is added, the activity comes to a maximum and then decreases. 
Temperature treatments change the concentration of caesium and of oxy- 
gen and consequently they also change the activity. Temperature may also 
affect the activity by changing the relative position of the caesium and 
oxygen. 

That the activity depends on whether the caesium is on top of the oxygen 
or the oxygen on top of the caesium, is not yet established by direct experi- 

* Langmuir, Phys. Rev. 22, 357 (1923). 


‘ Langmuir and Kingdon, Proc. Roy. Soc. 107A, 61 (1925). 
5 Becker, Phys. Rev. 28, 341 (1926). 
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ments. This view is based on the author’s interpretation of temperature 
changes in the analogous case of barium and oxygen on tungsten and on the 
theory of the cause of the increase in activity for composite surfaces.* 
On this view,** if oxygen is deposited on top of a monatomic layer of caesium 
on tungsten and stays there, the activity should decrease; but if thermal 
agitation causes the oxygen to settle down between the tungsten and the 
caesium, the activity should increase. Such surface rearrangements take 
place apparently at temperatures lower than those necessary to produce 
other surface changes such as surface migration, diffusion inward, or evapora- 
tion. Temperature treatments can thus change the activity by changing the 
amount of caesium, amount of oxygen, or their relative positions. But even 
with this complex filament, the activity is independent of plate potential 
and the amount of current sent into or taken out of the filament. 

Since coated filaments show changes in activity with temperature treat- 
ments in much the same way as do filaments with composite surfaces, it is 
quite natural to assume that coated filaments also have a composite surface 
of presumably metallic bariumf on the oxide. This view is supported by 
another striking similarity between oxide coated filaments and filaments 
which are known to have composite surfaces. In both cases the emission 
current increases quite appreciably as the plate potential is increased more 
and more, while the emission current from an atomically clean surface 
saturates quite well.*t 

Because of these two similarities, the author in private laboratory reports 
proposed the hypothesis that the enhanced activity of oxide coated filaments 
is due to metallic barium on the surface of barium oxide. Numerous ex- 
periments made in this laboratory have fully confirmed this hypothesis. 
The same or a similar hypothesis was proposed by Koller’ in 1925, Rothe® 
in 1926, Espe® in 1927, and Detels!® in 1927. 

This hypothesis takes into account the similarities in coated and com- 
posite surface filaments but it fails to take into account the dissimilarities. 
It has been emphasized above that the activity of a coated filament 
is changed by changes in plate potential and by the current sent into or 
drawn from it, while the activity of other forms of composite surface fila- 
ments is not. It follows that an additional hypothesis must be made for 


* See 15 below. 

** This view was held by Mr. Ryde of the G. E. Laboratory at Wembley, England in 
1928. He was kind enough to discuss the matter with the author on his visit to America. 

+ Throughout this article barium will be used as representative of the group of alkaline 
earths. It is to be understood that calcium and strontium act qualitatively like barium. 

* See “On Electrical Fields near Metallic Surfaces” by Becker and Mueller, Phys. Rev. 
31, 431 (1928). 

+ In oxide coated filaments poor saturation may sometimes be produced by additional 
causes such as heating effects in the oxide or localized space charges. 

7 Koller, Phys. Rev. 25, 671 (1925). 

* Rothe, Zeits. f. Physik 36, 737 (1926). 

® Espe, Wissenschaftliche Veréffentlichungen aus dem Siemens Konzern, 5, 29 and 46 
(1927). ; 
10 Detels, Jahrbuch der Drahtlosen Telegraphie und Telephonie 30, 10 and 52 (1927. 
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coated filaments. Electrolysis of the oxides naturally suggests itself since it 
has been known for many years that heated salts and in particular some oxides 
are electrolized. This was shown by the work of Bose," Horton,” and Span- 
ner.'* From what is known of the conductivity of oxides, it may safely be 
postulated that the current perpendicular to the surface of the oxide may be 
carried partly by electrons and partly by ions, presumably of barium and 
oxygen. 

This hypothesis too was proposed by the author in private laboratory 
reports. It had previously been held by Fredenhagen™ in 1913, who believed, 
however, that the emission of electrons was caused by chemical recombination 
of barium and oxygen. The electrolytic hypothesis was also held by Rothe,* 
Espe,® and Detels.!° They and others showed that an oxide coated filament 
which has been heated to a high temperature for a considerable time does not 
give off appreciable quantities of gas if no current is drawn from it, but that 
considerable gas is given off when curernt is subsequently drawn. Detels'® 
showed that this gas was oxygen by collecting it in a Geissler tube and show- 
ing that its spectrum was identical with that of asimilar tube filled with pure 
oxygen. J. E. Harris in these laboratories first showed that if large currents 
are drawn from a barium oxide+strontium oxide coated filament, barium 
and strontium are present in the core, while if a similar filament is glowed 
even more severely but without drawing currents, no barium or strontium 
are found in the core. It is thus clearly established by numerous and varied 
experiments that some electrolysis takes place in oxide coated filaments and 
that the products are the ones to be expected. It is proper, therefore, to 
attempt to account for the observed changes in activity with drawing cur- 
rent or sending current into the filament on the basis of the electrolytic 
theory. 

On the electrolytic theory, if electrons are sent into a coated filament from 
another source, barium is deposited on the outer surface of the oxide while 
oxygen is deposited on the interface between the core and the coating. 
Since thermionic emission depends on surface conditions and not on the con- 
ditions many layers beneath the surface, the barium should directly and im- 
mediately affect the activity while the oxygen should not do so. If barium 
on barium oxide produces the same effects as does barium on the tungsten, 
the activity should increase, pass through an optimum and then decrease as 
more barium is deposited on the surface. Experiments to be described below 
show that this is precisely what happens. These experiments will also be used 
to determine quantitatively the factors that determine the amount of elec- 
trolysis. 

If, on the other hand, electrons are drawn from a coated filament, barium 
is deposited on the core-oxide interface and oxygen on the outer surface. 
The barium should not affect the activity directly and immediately but may 


1 Bose, Ann. d. Physik 9, 164 (1902). 

12 Horton, Phil. Mag. 11, 505 (1906). 

18 Spanner, Ann. d. Physik 75, 609 (1924). 

4 Fredenhagen, Leipziger Berichte 65, 42 (1913). 
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do so if it subsequently diffuses to the oxide surface. The oxygen, however, 
should affect the activity directly and immediately. From analogy with 
caesium and oxygen on tungsten discussed above, the change in activity due 
to additional oxygen should depend on the amount of barium, the amount 
of oxygen, and also on the position of the oxygen relative to any surface 
layers of barium or barium-oxide. The experiments described below ful- 
fill this expectation. They show that when current is drawn from a coated 
filament, the activity sometimes decreases, while at other times it increases 
to a maximum and then decreases. Which of these behaviors is exhibited 
depends upon the degree of breakdown of the filament, the temperature, . 
and the plate potential. 

Another supplementary hypothesis is required to account for the com- 
plex behavior in coated filaments. It is to be expected that the barium and 
oxygen which are produced by electrolysis will diffuse into the oxide and 
from the oxide to the outer surface where they produce their effect on the 
activity. On account of the open or spongy structure of the oxides on these 
filaments, we should expect diffusion to take place very readily since it has 
been shown that surface migration of adsorbed atoms takes place at tempera- 
tures at which evaporation is negligible.” At temperatures near 1000°K 
there should be a continuous interchange of barium and oxygen atoms 
between the surface and the interior of the oxide. 

The various phases of the proposed hypothesis have been thoroughly 
tested, confirmed, and extended by the experiments which are described 
below. For this reason it seems proper to speak of the theory rather than the 
hypothesis of oxide coated filaments. 


EXPERIMENTAL TUBES 


Most of the experiments were performed with two tubes designated by 
12-J and 89-B. Figure 2 shows the plate and filament assembly of tube 12-J. 
The plates were made of preglowed nickel. The dimensions are given in the 
figure. Usually only the current to the center plate was measured; the end 
plates served as guard plates. F,; and F; are oxide coated filaments. The 
core consisted of a 95 percent platinum, 5 percent nickel alloy in the form of 
a ribbon 10.2X0.030X0.0054 cm. They were coated with about 20,000 
molecular layers of barium oxide and strontium oxide and were flashed at 
1000°K in air to form some barium and strontium platinates and nickelites. 
F, was broken down on the pumps, while F; was merely glowed but thoroughly 
degassed. During the course of the experiments, F; was subsequently broken 
down. F; is a tungsten ribbon 7.05 X0.0381X0.003 cm. Material from any 
filament can be vaporized onto any other filament. All the filaments were 
kept taut and straight by means of molybdenum springs not shown in the 
figure. This plate and filament assembly was inserted in one end of an appro- 
priate lead glass tube. At the other end a Buckley ionization manometer"® 
was attached. The tube was thoroughly pumped and had magnesium “get- 


1% Becker, “The Life History of Adsorbed Atoms and Ions” Proc. Am. Electrochem. Soc. 
55, 153 (1929). 
16 Buckley, Proc. Nat. Acad. Sci. 2, 683 (1916). 
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ter.” After it was sealed off and cleaned up by glowing the thoriated tungsten 
filaments in the manometer and the broken down coated filament, the pres- 
sure was ordinarily about 1X10-* mm Hg and the activity of the filaments 
was not affected by gas even over long periods of time. 

Tube 89-B was similar to 12-J. It differed from it in three significant res- 
pects: (1) The tungsten ribbon Fw was about 10 cm long instead of 7.05 cm. 


PLATES AND FILAMENTS 
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This reduced the effects of end cooling. (2) The middle collector or plate 
was 1.27 cm high instead of 2.53 cm. This too reduced end-cooling effects 
and also minimized effects due to non-uniformities along the filament. (3) 
A movable nickel shield could be placed between the unbroken down coated 
filament and the other two filaments. This permitted experiments which 
could not be performed in 12-J. 


EXPERIMENTAL RESULTS 


The experimental results will be divided into eight groups which tend to 
show that: (1) Metallic barium evaporates from a well broken down coated 
filament; (2) Oxygen evaporates if an emission limited current is drawn and 
the temperature is high enough; (3) The high activity of oxide coated filament 
results from metallic barium on the oxide surface; (4) Barium diffuses readily 
to and from the surface of the oxide; (5) Oxygen affects the activity in various 
ways; (6) Oxygen diffuses readily to and from the surface of the oxide; 
(7) Most of the current through the oxide is carried by electrons but a small 
part is carried by ions; the ratio depends upon the temperature, the composi- 
tion of the oxide, and the plate potential; (8) The complex time changes that 
occur can be interpreted in terms of these processes. 
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(1) Evaporation of metallic barium from oxide coated filaments. That 
metallic barium evaporates from a well broken down coated filament was 
shown by the following experiments. In tube 12-J, after filament F; had been 
broken down and thoroughly aged, it was glowed at 1260°K for about 10 
minutes while the tungsten filament Fw, which had previously been cleaned 
off, was cold. F; and the whole tube were then allowed to cool and the activity 
of Fw was determined at a series of testing temperatures which were as low 
as possible so that no appreciable time changes took place during the test. 
It was found that the activity of Fw had increased appreciably. F; was then 
again glowed as previously while Fw was cold. When the tube was cool, 
the activity was again determined at a series of low testing temperatures. 
This process was repeated a number of times. Whenever necessary, the emis- 
sion from Fw was extrapolated or interpolated to 1100°K. 

When the log of these emission currents reduced to a cm? were plotted 
against the total time of flashing F; at 1260°K, a curve similar to Fig. 1 
was obtained: Log 7 rises rapidly at first, then more and more slowly until it 
passes through an optimum beyond which it decreases. By appropriate 
experiments, it was shown that in equal intervals of time of glowing F;, equal 
quantities of material were deposited on Fw so that 6, the fraction of its 
surface covered, is directly proportional to the time of glowing F; at a 
fixed temperature. That @ is equal to unity for the optimum activity is based 
on (1) analogy with caesium on tungsten for which it was proven® and (2) 
the author’s theory of the emission from composite surfaces’ from which it 
follows that @ should be equal to unity at the optimum activity for all electro- 
positive materials adsorbed on strong electronegative surfaces. The exact 
shape of the curve in Fig. 1 was obtained by correcting the observed curve 
for non-uniform deposition on various portions of one side of Fy. 

When Fw was subsequently glowed at 1100°K, the activity increased 
until it reached practically the same optimum as before; then the activity 
decreased and approached a steady value. When the temperature of Fw 
was increased again, the activity decreased, at first rapidly and then more 
and more slowly. This process was continued at successively higher tempera- 
tures until the activity of clean tungsten was reached at a temperature of 
about 1700°K. 

Apparently something evaporated from F; and activated Fw. When 
Fy was subsequently heated, this something disappeared from the surface of 
Fw. If this whole behavior could be repeated even in some of its quanti- 
tative aspects, with a source known to contain barium, it would be fair to 
infer that the material that evaporated from a well broken down coated 
filament is barium. 

This inference was established by an experiment performed by J. M. 
Eglin?’ of these laboratories. He built and worked with tube 5-E which was 
just like 12-J except that the two coated filaments were replaced by two fila- 
ments consisting of an alloy core of 95 percent platinum, 5 percent nickel 
and metallic barium. The barium was obtained from barium-azide and was 


17 Eglin, Phys. Rev. 31, 1127 (1928). 
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alloyed with the filaments in a separate tube by a process which is now well 
known.? If such barium alloy filaments are thoroughly aged to free them of 
nitrogen, and are then glowed at 1200°K, chiefly barium evaporates from 
them. The changes in activity of the tungsten filament obtained in this tube 
were just like those described in the preceding paragraphs for tube 12-J. 
Even more significant is the fact that the optimum activity obtained in these 
two tubes agreed quite closely with each other. 

Figure 3 is in part intended to show how close this agreement is. It is 
a plot of the emission vs. the power supplied to Fw per cm? on a coordinate 
paper devised by C. J. Davisson of these laboratories. If the emission from a 
filament obeys Richardson’s 7'/? equation and if the watts dissipated in- 
crease with 7‘, then the “power emission” line should be straight. Line No. 
1 shows the observed optimum for tube 12-J. Line No. 2 is for tube 5-E. This 
line has been slightly corrected for the fact that the tungsten ribbons in the 
two tubes were not quite identical. This correction, which amounted to 
about a factor of 8 at the lower currents and 2 at the higher currents, was 
based on the power emission lines for clean tungsten and on the heating 
currents necessary to make the filaments barely visible in a dark room. Line 
No. 3 is for a similar optimum from a coated filament in tube 89-B. Note 
that none of the three lines is straight. This is due to the fact that at temp- 
eratures below about 1000°K, the end cooling is appreciable and becomes 
more marked at lower temperatures. Line 4 is line 1 corrected for end cooling 
by approximate methods. Since in 89-B the end-cooling correction should 
be less than in the other two tubes because the filament was 10 cm long in- 
stead of 7.05 cm line 3 should be everywhere higher than lines 1 and 2. 
The figure shows that this is so. The agreement between line 2 for barium 
on tungsten and lines 1 and 3 for two cases of material vaporized from well 
broken down coated filaments on tungsten, is almost as good as between 
two successive determinations of the same line. Line 5 in the lower right 
corner shows the emission from clean tungsten. An approximate tempera- 
ture scale is added below the power scale. The close agreement of the op- 
timum lines and the similarity of the general behavior in the three tubes make 
it very probable that barium evaporates from a well broken down coated 
filament at temperatures near 1250°K. 

If in similar experiments the coated filament were not well broken down 
or if its temperature were too high, or if current was being drawn from it 
during the vaporization process, then another element presumably oxygen 
either free or in the form of an oxide, was vaporized over with the barium. 
This oxygen showed itself in the following four ways: (1) The log 7 vs. @ curve 
of Fig. 1 instead of rising rapidly and then more slowly, rose slowly and then 
approached the trend shown in Fig. 1. If the oxygen ratio was rather high, 
the activity actually decreased during the first few minutes of flashing; 
then it started to rise. (2) If about 0.9 layer was deposited on Fw and if Fw 
was then glowed at successively higher temperatures, the activity decreased 
to a value below that of clean tungsten. At still higher temperatures, near 
1800°K the activity increased again to that of clean tungsten. This behavior 
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is analogous to that for caesium and oxygen on tungsten.® (3) If more than 
one layer was deposited, the temperature required to pass backward through 
the optimum at a reasonable rate was considerably higher than if no oxygen 
was vaporized. The oxygen apparently holds the barium more tenaciously 
than does clean tungsten. This too is analogous te its behavior for caesium 
on oxygen on tungsten. (4) The optima were not definite but depended on 
the amount of oxygen coming over with the barium. Comparatively small 
amounts of oxygen increased the optimum line to a value of about five times 
that for barium on tungsten. Large amounts of oxygen lowered the optimum 
line more and more. This also is analogous to caesium on oxygen on tungsten, 
but in this case, the optimum caesium on tungsten line can be increased 
several hundredfold by the proper amount of oxygen. 
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An interesting deduction from these experiments is that barium and oxy- 
gen do not permanently combine chemically on a hot tungsten surface. An- 
other way of stating this, is that barium oxide is decomposed on a hot 
tungsten surface. 

The peak or optimum activity phenomenon furnishes us with a convenient 
and powerful method for determining rates of evaporation of electropositive 
materials at various temperatures. All that need be done is to determine the 
time required for the activity of Fw to reach its optimum and to compute from 
the geometrical relationships the number of layers that must evaporate from 
the primary to deposit one layer on the tungsten. In this way the number 
of layers evaporating per minute at a series of temperatures was determined 
for several coated filaments. 

In Fig. 4 the log of the evaporation rate expressed in layers per minute 
has been plotted as a function of 1/7. Curves 1, 2 and 3 are for the evapora- 
tion of barium from oxide coated filaments.* Line 4 represents the data 


* 6, the fraction of the surface covered with barium may have been different from one 
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obtained by Reerink for barium oxide from barium oxide which was repor- 
ted in an article by Zwikker.'* The points through which line 5 has been 
drawn are experimental values obtained by C. J. Davisson in these labora- 
tories several years ago. He showed also that the rate of evaporation of stron- 
tium oxide from strontium oxide was negligibly small compared to barium 
oxide from barium oxide. The filaments which furnished the data for line 
1, 2 and 3 consisted of a mixture of barium oxide and strontium oxide. Pre- 
sumably the evaporation rate of barium oxide from this mixture is somewhat 
less than it would have been from a pure barium oxide surface. From this 
and from a comparison of the slopes of lines 1 to 3, and lines 4 and 5, it fol- 
lows that at low temperatures chiefly barium should evaporate while at suffi- 
ciently high temperatures chiefly oxide should evaporate. f 

It should be pointed out that this method of determining evaporation 
rates at various temperatures might well be used for a study of the precise 
way in which the evaporation depends upon the temperature. 

That mostly metallic barium rather than barium oxide evaporates from 
well broken down coated filaments follows from experiments performed by 
C. J. Davisson in these laboratories several years ago. It was found that the 
material that evaporated produced rather large resistance leaks across glass 
surfaces between metallic terminals. When air was let into such tubes, 
the resistance increased greatly and tapidly. Commercial hydrogen, pre- 
sumably slightly impure, caused no rapid changes in resistance. Metallic 
barium would be expected to produce these results while barium oxide would 
not. 

These experiments, it seems to us, prove conclusively that metallic barium 
evaporates from a well broken down oxide coated filament. If the tempera- 
ture is not too high and if no electrons are drawn from the filament, it ap- 
pears that practically nothing else evaporates. At sufficiently high tempera- 
tures barium oxide must also evaporate. 

(2) Evaporation of oxygen. When current is drawn from a coated 
filament at a sufficiently high temperature and if this current is not limited 
by space charge, oxygen evaporates. Some of this oxygen is deposited on the 
tungsten ribbon in Fig. 2 and decreases its activity. Thus, when filament 
F; in tube 12-J was at a temperature of about 800°K and the plates were at 
a positive potential of about 100 volts, it was found that the activity of the 
tungsten ribbon decreased progressively. In one particular case, the current 
to the center plate was about 1.0 milliampere and the activity of the tungsten 
ribbon was decreased by a factor of two for every 2 hours. By drawing large 
currents we have been able to reduce the activity of the tungsten filament by 
a factor of 200.* If the tungsten ribbon was subsequently raised to tempera- 





filament to the next but for a given filament it was probably quite constant since the activity 
did not change appreciably with changes in glowing temperature. 

18 Zwikker, Physica 8, 241, (1928). 

t At about 1250°K the evaporation rate of barium from metallic barium is about 10* 
times as great as the rate of evaporation of barium from a well broken down coated filament. 

* From the fact that the activity of the tungsten ribbon can be reduced by more than a 
factor of two, we may conclude that the oxygen which strikes one side of the ribbon must 
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tures higher than about 1600°K, the activity increased for a time and then 
became almost constant. If the temperature was raised again, the activity 
increased rapidly at first and then more and more slowly.** At 1800°K, the 
activity reached that of clean tungsten in something like ten minutes. 

Practically all of this oxygen must leave the coated filament in an un- 
charged condition, for the rate at which it collects on the tungsten ribbon is 
not affected by the potential of the tungsten ribbon. If an appreciable amount 
of the oxygen were negatively charged and if the tungsten ribbon were made 
negative with respect to the coated filament, the rate of deactivation of the 
tungsten ribbon should be appreciably decreased. We have made several 
experiments in several tubes to test this point and have never found any 
evidence that the oxygen was negatively charged. This is rather surprising 
in view of Barton’s result!® which is frequently quoted in the literature. 
It may be that the condition of his filament was considerably different from 
that of ours and from that in ordinary tubes. 

By determining the rate of deactivation of the tungsten ribbon as a re- 
sult of various treatments imposed upon the coated filaments, it is a rather 
simple matter to determine the factors upon which the rate of evaporation 
of oxygen depends. Such experiments have shown that the rate of oxygen 
evaporation increases with (1) the temperature of the filament, (2) the amount 
of current drawn from it, (3) the plate potential and (4) the composition or 
state of the oxide. Thus far, the exact manner in which it depends on these 
factors has not been determined accurately. The first two factors were to 
be expected. That the oxygen evaporation increases rapidly with the positive 
plate potential is somewhat surprising since one might expect the potential 
difference across the oxide at a given temperature to depend only on the cur- 
rent through it. Experiment, however, shows very definitely that if the emis- 
sion is limited by space charge, no oxygen or very little oxygen is liberated 
compared to the amount liberated when the plate voltage is increased such 
that the current is limited by emission. When the plate is at a negative po- 
tential, no oxygen evaporates. It should be pointed out that the temperatures 
at which it is possible to obtain oxygen at an appreciable rate are considerably 
lower than those required to obtain barium at an appreciable rate. It is 
therefore possible to obtain oxygen without barium by one procedure, and 
barium without oxygen by another. 

An important observation was this: it was comparatively easy to obtain 
oxygen from filament F;, the unbroken down coated filament in tube 12-J, 
soon after the tube was sealed off the pumps. However, after a considerable 
quantity of oxygen had been removed and the filament had become more 
active, it required large currents, 10 to 20 mils, and high plate potentials, 
150 to 200 volts, to obtain appreciable quantities of oxygen. We take this 





migrate and reach the other side of the ribbon at the testing temperatures of about 1400°K. 
Subsidiary experiments especially with the shielded tube 89-B showed that no oxygen could 
reach the back side of the tungsten ribbon by being reflected from the walls of the tube. 

** The reason for this behavior is that the rate of evaporation at a given temperature 
decreases rapidly as @ decreases. 

19 Barton, Phys. Rev. 26, 360 (1925). 
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to mean that originally an appreciable part of the current through the oxide 
was Carried electrolytically by oxygen and barium ions; after the filament was 
broken down and barium was dispersed through the oxide, a much smaller 
part of the current was carried by ions. 

From the experiments described in this section, we conclude that when 
current is drawn from a coated filament at a sufficiently high temperature 
and this current is not limited by space charge, oxygen evaporates. Most 
of this oxygen must be uncharged. The rate at which it appears increases 
rapidly with the plate potential, and it decreases as more barium is dispersed 
through the oxide. 

(3 Activation due to barium on barium oxide. There are two ways by 
which barium can be caused to appear at the surface of the oxide: one is by 
sending electrons into the oxide from an outside source and thereby elec- 
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trolizing the oxide; the second is by vaporizing barium from another fila- 
ment which, under proper conditions, may be another oxide coated filament. 
The barium which is thus brought to the surface of the oxide should change 
the activity of the filament. If barium on barium oxide is analogous to bar- 
ium on tungsten the log of the activity should increase rapidly at first, then 
more and more slowly until it passes through an optimum beyond which it 
should decrease. When the oxide coated filament is subsequently glowed and 
the barium is thereby caused to diffuse or vaporize from the surface, the 
activity should again increase to an optimum and then decrease to its normal 
value. The following experiments and figures show that this is exactly what 
happens. 

Figure 5 shows a logarithmic plot of the activity of a strontium oxide un- 
combined filament as a function of the amount of electricity, expressed in 
microampere minutes, sent into it from an outside source through a P. D. 
of 140 volts. The activity is here defined as the emission current from the 
entire 10 cm length of the filament at approximately 600°K and 70 volts 
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plate potential. For some points this activity was determined by extra- 
polation from higher or lower temperatures. Note that the activity increased 
to an optimum value which was about 30,000 times its initial value after 
about 1200 microampere minutes had been sent into the filament. Subsi- 
diary tests showed that all parts of the filament did not reach their optimum 
at the same time. This causes the peak to be considerably flatter and some- 
what lower than it would otherwise be. Such non-uniformities are to be 
expected and are extremely difficult to avoid. If they had been avoided, we 
estimate that the optimum would have come at 300 mic-minutes. Subsi- 
diary tests also showed that the initial test point in Fig. 5 did not correspond 
to clean strontium oxide, for subsequent more severe glowing reduced the 
activity below this value. If the current through the oxide were entirely 
electrolytic and if the microscopic surface area were the same as the measured 
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macroscopic area, about 2.0 microampere minutes per cm? should suffice to 
cover the surface with a monatomic layer of metallic strontium. From these 
facts we may conclude that very roughly 1/200 of the current through this 
particular oxide was carried by strontium ions. 

The first part of Fig. 6 shows a similar curve for test 17, tube 89-B, fila- 
ment F;. In this case 21 microamperes were sent into the coated filament for 
a total of 22 minutes. At the end of this time F, was glowed at 800°K for 
23 minutes. The activity was tested under standard conditions at the times 
for which points are plotted. The activity passed through an optimum which 
is about 25 percent lower than the optimum for the preceding curve. This 
small discrepancy in the optima is nogreater than is usual in this type of curve 
and is possibly due to slight differences in the degree of uniformity of cover- 


ing. 
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If we compare Figs. 5 and 6 we note that the peak is much sharper in 6 
than in 5. This is undoubtedly due to the fact that in Fig. 6 the degree of 
non-uniformity at the peak is much less than it is in Fig. 5 since the activity 
in Fig. 6 refers to 1.27 cm of the center portion of the filament while that in 
Fig. 5 refers to the entire 10 cm length. 

A still more striking difference between these two figures is that in Fig. 6 
the optimum activity is only about 13 times the initial activity while in 
Fig. 5 it is about 30,000 times the initial activity. This is undoubtedly due 
to the fact that in Fig. 6 the filament surface was nearly completely covered 
with barium at the initial activity. 

The proof of this statement is based on test 61A, tube 89-B. In this test 
the broken down coated filament F,; was flashed at 1450°K in order to va- 
porize barium onto F; and also onto Fy. (See Fig. 2 for the relative positions 
of the oxide coated and tungsten filaments.) The activity of F, changed in 
the manner practically identical with that given in the first part of Fig. 7 for 
test 65 except that the curve was extended for 20 minutes instead of 5; from 
5 to 20 minutes it decreased at a slower and slower rate. During the same 
time the activity of Fy followed a course given by Fig. 1, the peak occurring 
after 14 minutes of flashing F,. Note that F, reached its optimum activity 
before two minutes of flashing while Fyreached its optimum after 14 minutes. 
Since F, was about twice as far from F, as was Fy from F,, the amount of 
material reaching F, in a given time is about 1/2 as great as that reaching 
Fy. Consequently it required only 1/14 or about 7 percent as much material 
to bring F, toits optimum as it did tocover Fy with a monatomic layer. On 
the hypothesis that the optimum activity occurs at a monatomic layer, we 
may conclude that F, was 93 percent covered at the initial value in Fig. 7. 
This initial value could be repeatedly obtained after glowing at about 1300°K 
for several minutes, and was roughly the same as in Fig. 6. 

Fig. 7 shows the change in activity of F, while F; was glowed at 1460°K 
for 5 minutes. Note that the activity passes through an optimum after bar- 
ium is deposited on F,. Subsequent glowing of F2, at the temperatures and 
times indicated in Fig. 7, caused the activity to pass through practically the 
same optimum as before. 

It is very instructive to compare Figs. 6 and 7. In one case barium was 
brought to the surface by electrolysis while in the second case it was deposited 
by vaporization. In both cases the activity of F, passed through an optimum. 
In both cases subsequent glowing caused the activity to go backward through 
an optimum practically equal to that of the optimum which preceded it. 
The slight differences between the two cases are easily accounted for by 
the treatment given the filament between tests 17 and 65. Unfortunately we 
do not yet have two such curves for the filament in identical conditions and 
especially in an initially rather inactive state. 

That the amount of barium brought to the surface by electrolysis is 
directly proportional to the quantity of electricity sent into the oxide coated 
filament may be inferred from Fig. 8. This is a plot of the logarithm of the 
emission from the central inch of filament F; at a temperature of 475°K and 
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a plate voltage of 31.5 volts vs. the mic-minutes of current sent into F; while 
F; was at 475°K. Each set of points shown in this figure represents a com- 
plete curve by itself. For one set of points, the current sent into the filament 
was 6.0 microamperes; for others, 15, 60, and 380 microamperes respectively. 
From the fact that the different sets of points all fall on the same curve, we 
may conclude that the number of ions brought to the surface depends only 
on the product of the current and the time and not upon either variable alone. 
This is readily understandable if the amount of barium brought to the surface 
is directly proportional to the quantity of electricity sent into the filament. 

In Fig. 8 the initial activity is very close to the optimum activity. While 
we have no proof in this particular instance that this means that the surface 
is practically completely covered, we have no hesitation in inferring this from 
analogy with the data discussed above for tube 89-B. Similar data in tube 
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12-J could not be obtained since the filament F; was burned out rather early 
in the history of the tube. Before the curve shown in Fig. 8 was taken, fila- 
ment F; had been well activated by breaking it down. Its activity could not 
be decreased appreciably by glowing it even at fairly high temperatures such 
as 1300°K for a time sufficient to vaporize hundreds of layers of barium from 
it. Possibly by glowing it at temperatures near the melting point of the core, 
we might have decreased the activity slightly. 

While we have no curve showing the rise in activity of F; while barium 
was being vaporized onto it from F,, we do have a curve which shows what 
happened to the activity of F; after about three layers of barium had been 
deposited on it. This is shown in Fig. 9 in the form of power emission curves. 
The line or lines on the extreme right show the initial activity in test 1 soon 
after the tube was pumped. After about three layers of barium had been 
vaporized onto F;, the emission was that given by points 1 and 2 in the lower 
left-hand corner of the figure. Even at as low a temperature as 500°K the 
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activity increased due to barium disappearing from the surface. At any 
temperature such time changes were usually rapid at first and then slowed 
down until they were fairly steady.* The temperature was then changed to 
a higher value at which the activity increased again. This process continued 
up to points 7 and 8. For still higher temperatures the activity decreased at 
each temperature. After point 12 the temperature was decreased. At the 
lower temperature the activity was now quite steady. In this condition the 
activity was still appreciably higher than it was in the initial state, pre- 
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sumably because some barium was still left on the surface even after glowing 
at approximately 1000°K. Note that the maximum line is about 100,000 
times as active as the initial line. By subsequent treatment of the filament, 
the activity was raised almost to that of the optimum line given in the upper 
left corner. This line is about 10’ times as active as the initial line. 

The evidence presented in connection with Figs. 5, 6, 7, 8, and 9, proves 
conclusively that the activity of an oxide coated filament can be greatly 
modified by depositing barium on its surface. We could give many more 
curves similar to the ones which have been given, but feel that this is un- 
necessary. 

The amount of barium on the surface of the oxide is not the only factor 
which determines the activity, for otherwise the activity of all oxide coated 
filaments should be the same for an optimum amount of barium on the sur- 
face. Experiment shows that this is not so. The optimum line for filament 


* The reason for this behavior is that the rate of evaporation at a given temperature 
decreases rapidly as @ decreases. 
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F; in tube 12-J is shown in Fig. 9. The optimum line for F,; in tube 89-B 
was approximately parallel to this line but was a factor of 25 below it. The 
optimum line for filament F; tube 89-B, was again approximately parallel 
to the optimum line for F; tube 12-J, but was a factor of 500 below it. That 
the optimum lines for the same filament at various stages ofits breakdown 
process may differ is suggested by the position of the “maximum line” and 
“optimum line” in Fig. 9. There is some question as to the part that oxygen 
played in the determination of the “maximum line,” but for this line as well 
as for the “optimum line,” the surface was quite likely covered with an opti- 
mum amount of barium. 

(4) Diffusion of barium in the oxide. That the barium which is produced 
by electrolysis diffuses into and through the oxide can be deduced from a 
number of experimental facts. One of these is the activation resulting from 
the breakdown process. The barium is deposited on the oxide-core interface 
and must diffuse through the oxide to the surface before it can increase the 
activity. 

Another argument for diffusion is this; hundreds of layers of barium can 
be vaporized from a well broken down filament without appreciably altering 
its appearance or its activity. This barium could not possibly have been on 
the surface of the oxide at one and the same time, for otherwise, (1) the sur- 
face would have had a metallic appearance, (2) the activity could not have 
been affected by adding more barium, and (3) the rate of vaporization would 
have been much higher than was found to be the case. Hence the barium 
must have diffused through the oxide to the surface. 

A semi-quantitative proof of diffusion is based on the life of these fila- 
ments. It is a well-established fact that these filaments retain an appreciable 
percentage of the original oxide even after furnishing a current of 20 milliam- 
peres from a cm? for at least two years or approximately 15,000 hours. We 
shall show below that more than one part in one thousand of the current 
through the oxide is carried by ions. Hence at 15,000 hours, 300 milliampere 
hours or 18,000 milliampere minutes of ion current must have reached the 
surface. If each ion carries a charge equal to that on the electron and ifa 
single layer of barium corresponds to about 7 X 10" atoms per cm*, it requires 
about two microampere minutes of ion current per cm? to produce one layer. 
Hence the number of layers of barium that should be deposited at the current 
and time specified above is more than 9,000,000. Actually there are only 
about 20,000 layers of oxide on a filament and only a portion of these dis- 
appear during the treatment. Hence it follows that by far the largest part, 
about 0.998, of the barium which is produced by electrolysis must diffuse 
into the oxide, and must recombine with oxygen to form barium oxide. 

The rate at which the barium diffuses from the surface into the oxide may 
be determined from curves such as that shown in the second part of Fig. 6. 
In this case the surface concentration of barium was greater than that cor- 
responding to equilibrium conditions and hence the barium disappeared 
from the surface rapidly at first and then more slowly. This barium either 
evaporated or diffused into the oxide. That only a small part evaporated 
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appears from the following calculation. In Fig. 6 we estimate that at 800°K 
about 0.05 layers left the surface between times 30 and 45 minutes. This 
corresponds to a rate of one layer in 300 minutes or about 3X10~* layers 
per minute for an average @ of about 0.98. From curve 3 in Fig. 4 for which 
6 was about 0.90 we can obtain that at 800°K, the rate of evaporation was 
about 3X10-" layers per minute or only 1/10,000 of the rate of disappear- 
ance. It seems improbable that the evaporation rate increases as much as 
one thousand-fold when @ changes from 0.90 to 0.98. Hence most of the bar- 
ium must have left the surface by diffusion and the diffusion rate must be ap- 
proximately the same as the total rate of disappearance. 

Whenever the barium surface concentration exceeds the equilibrium con- 
centration, the rate at which the surface concentration approaches the equili- 
brium concentration is rapid at first and then decreases more and more. 
This is readily deduced from the shapes of the glowing curves like those in 
Figs. 6 and 7. One would hardly expect that the diffusion inward should cease 
when the equilibrium concentration is reached. It appears very likely that 
there is a constant interchange of atoms between the surface and the interior 
of the oxide. We should expect that the rate of diffusion to the surface should 
be directly proportional to the volume concentration of barium in the oxide 
while the rate of diffusion from the surface should increase exponentially 
with the absolute value of the surface concentration. The latter expectation 
is based on the picture that it takes a definite amount of work to take an 
atom from the surface into the interior and that this amount of work de- 
creases as the surface concentration increases; in other words, that the dif- 
fusion inward process is analogous to the evaporation process.* When a 
steady state is reached at any glowing temperature just as many atoms come 
to the surface as leave the surface. From this picture we should expect that 
the surface concentration should not differ greatly when the volume concen- 
tration is changed by even as great a factor as 10. This may explain why 
the activity is not very greatly dependent on the amount of barium stored 
in the core and oxide coating. 

(5) The effect of oxygen on the activity. It is a well-known fact that 
when current is drawn from a coated filament the activity of this filament 
changes rapidly at first and then more and more slowly. The magnitude and 
even the direction of the change in activity depend greatly upon the composi- 
tion of the oxide, the plate potential, the degree of breakdown and the tem- 
perature of the filament. 

For a well broken down filament we have invariably found that the ac- 
tivity decreases rapidly and then more and more slowly. This is illustrated in 
the first part of Fig. 10. Curve 1 shows how the emission current in milli- 
amperes decreased at a temperature of 715°K and 60 volts plate potential. 
At intervals the activity was tested at a lower temperature, 475°K. This 
activity is plotted in microamperes as curve 2. 

It was shown above that if the temperature and the plate potential are 
high, oxygen is formed by electrolysis and evaporates. Presumably, there- 


* See reference 15. The reason for this behavior is, that the rate of evaporation at a given 
temperature decreases rapidly as @ decreases. 
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fore, the decay in activity is caused by oxygen which has come to the surface 
and some of which stays there because the temperature does not remove it 
rapidly. Because of its electron affinity some of the oxygen should be nega- 
tively charged and if the oxygen stays on top of the surface layer of barium 
it should produce electrical fields which hinder the escape of electrons. Con- 
sequently the activity should decrease. If more and more oxygen is brought 
to the surface the rate at which it disappears from the surface either by 
evaporation or by diffusion should increase rapidly until the rate of disap- 
pearance is equal to the rate of production. We thus have a satisfactory pic- 
ture of the cause of the decay in activity when current is drawn. 

This view is supported by experiments in which oxygen was vaporized on- 
to the broken down filament from another oxide filament. The activity of the 
broken down filament was reduced by more than a factor of ten. The slope 
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of the log 7 vs. 8 curve was practically constant while the activity decreased 
by a factor of five; then the slope gradually became less. 

Of course, one could postulate that the oxygen which comes to the surface 
combines with the barium and decreases the activity by decreasing the sur- 
face barium concentration. This, however, is unlikely since it was shown above 
that barium and oxygen do not permanently combine chemically when they 
are present on a hot tungsten surface. Even more convincing is the experi- 
mental fact (test 34, tube 89-B) that when the barium concentration was 
more than the optimum and oxygen was then brought to the surface, the 
activity decreased rather than increased. If the oxygen combines with the 
barium, the barium concentration should decrease and thus come nearer to 
the optimum concentration which should result in an increase in the activity. 

If the filament is merely glowed at the same temperature at which cur- 
rent has been previously drawn the activity increases rapidly at first and then 
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more and more slowly. This is illustrated by curve 3 in Fig. 10. In general 
the activity either never or only extremely slowly reaches the initial activity. 
If, however, the filament is glowed at a high temperature such as 1200°K, 
the initial activity can be readily restored. 

The explanation for this activation on glowing is fairly apparent. Some 
of the oxygen which was deposited on the surface by electrolysis is caused 
to leave the surface by thermal agitation. As the oxygen surface concentra- 
tion gets nearer and nearer to its equilibrium value, the rate of disappear- 
ance of the oxygen decreases more and more. 

Thus far we have described the effects of oxygen on the activity in a 
well broken down filament. We shall now proceed with the effects of oxygen 
in filaments which have not been broken down. The best data on this sub- 
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ject were obtained with an uncombined strontium oxide filament. The 
results were confirmed with an uncombined barium oxide filament and to 
some extent in a barium oxide plus strontium oxide combined filament. They 
are thus quite general and well established even though similar results 
have never been reported in the literature. 

When current was drawn from the unbroken down strontium oxide fila- 
ment in tube 105-J at a temperature less than 950°K and a plate potential of 
178 volts, the log of the emission increased slowly, then more rapidly, then 
more slowly again until it reached a peak or maximum value beyond which 
it decreased and asymptotically approached a steady value. A typical curve 
is shown in the first part of Fig. 11. The rate of activation and the height of 
the peak increased considerably with the plate potential. For sufficiently 
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small plate potentials such as 30 volts the activity never passed through a 
peak value but merely increased to its steady value. 

When the filament was glowed at the same temperature at which current 
had been drawn, the activity decreased rapidly at first, then more and more 
slowly. This activity was determined under testing conditions which did 
not alter the activity: either both the temperature and the plate potential 
were decreased, or else the plate potential alone was decreased to a low 
value. The solid line in the second part of Fig. 11 shows the decrease in the 
activity when the plate potential was reduced to 5.0 volts. The dashed line 
above this solid line shows the instantaneous emission that would have oc- 
curred if the plate potential had been changed back to 178 volts at any par- 
ticular instant. This was actually done at time 65 minutes with the result 
shown in the third part of the figure. A peak was again obtained, but this 
peak was lower and occurred later than the first peak. In order to duplicate 
the first curve it was found necessary to glow the filament at a much higher 
temperature such as 1200°K for 5 or 10 minutes. 

If the data of the first curve in Fig. 11 are used to plot a curve of the log of 
the activity vs. mic-minutes of current drawn from the filament, a straight 
line is obtained for the initial portion of the curve. The slope of this line is 
a measure of the rate at which some activating material is coming to the sur- 
face. A series of experiments has shown that the slope of this line (1) in- 
creased rapidly with the plate potential, (2), decreased as the temperature 
increased, (3) increased when the surface was slightly covered with an active 
material vaporized from the plates, (4) decreased if the filament was not thor- 
oughly glowed after the preceding test, (5) decreased markedly as the degree 
of breakdown increased. The height of the peak emission divided by the 
initial emission varied in this same manner with these five factors. The 
mic-minutes required to reach the peak (1) decreased with the plate potential 
(2) decreased as T increased, (3) and (4) changed slightly with the third and 
fourth factor given above, (5) decreased from about 20 to about 2 mic- 
minutes per cm? when the filament was broken down. In the case of a barium 
oxide uncombined filament in tube 123-J the peak came at about 60 mic- 
minutes per cm?. 

At first sight it might appear that this peak was due to barium as in the case 
of Figs. 5, 6 and 7. The following reasons show that this cannot be so: (1) 
the barium is deposited on the core and the temperature is so low that it 
seems highly improbable for the barium to diffuse rapidly enough to the sur- 
face of the oxide to account for the rapid increase in activity; (2) the glowing 
curve in Fig. 11 decreases rather than increases as it did in Figs. 6 and 7; 
(3) the activity at the peak is very variable and depends on the conditions 
prevailing during the treatment whilein the bariumcase the optimum activity 
does not depend on these same conditions; (4) the optimum activity due to 
barium is about 50,000 times the initial activity while for this peak pheno- 
menon the peak activity is less than 1000 times the initial activity. 

The hypothesis which we tentatively propose for this peak on drawing 
current is an extension of what is already known. Due to electrolysis, oxygen 
is brought to the surface. As will beshown below the factors that determine 
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the amount of electrolysis are of the same sort as those that determine this 
peak phenomenon namely, plate potential, temperature and composition of 
the oxide. Hence it is natural to ascribe this peak phenomenon to oxygen. 
From analogy with caesium on tungsten it is to be expected that a little 
oxygen increases the activity while more oxygen decreases it. Whether the 
activity is increased or decreased depends upon just where the oxygen is 
located with reference to the surface barium atoms. If the oxygen is beneath 
the barium layer, the activity is raised because the negative oxygen induces 
more positive barium ions and the moment of the additional barium ions 
and their images more than counterbalances the moment of the oxygen ions 
and images. If on the other hand, the oxygen is on top of the surface layer of 
barium, the activity is decreased because now the moment due to the oxygen 
is greater than the moment due to the barium. 

By combining the hypothesis with the experimental facts we may draw 
the following conclusions: oxygen is deposited by electrolysis just beneath the 
partial surface layer of metallic barium or strontium. If the temperature is 
high enough some or all of this oxygen rather rapidly arrives at the top of 
the barium layer where it forms an outer oxygen layer which is a more stable 
layer than the inner oxygen layer. If the temperature is low enough the 
first oxygen that is deposited stays in the inner layer. As the oxygen concen- 
tration increases, the field tending to pull electrons and negative oxygen ions 
toward the plane of the barium ions increases and some of the oxygen reaches 
the top of the barium layer. At the peak the factors tending to increase the 
activity are balanced by those tending to decrease it. Beyond the peak the 
net result of electrolysis is the increase in concentration of the outer layer of 
oxygen. We believe that some of the oxygen also diffuses inward into the 
oxide. If the filament is glowed either before or after the peak, oxygen leaves 
the inner layer; it reaches the outer layer or diffuses into the oxide. 

To be consistent we need to apply this hypothesis to the case of well bro- 
ken down filaments in which emission is accompanied by a decrease in ac- 
tivity. For well broken down filaments the field which tends to pull elec- 
trons and also negative ions into and beyond the plane of the positive Ba 
ions is very great so that even at very low temperatures the oxygen is rapidly 
pulled into the outer layer. It may also be that in well broken down filaments 
the inner oxygen layer is already near its maximum concentration, and any 
additional oxygen is forced to go to the outer layer. This latter view would 
fit in well with the observed fact that the maximum activity which can be 
obtained depends on one or more factors besides the optimum amount of 
barium. The study of the effect of oxygen on the activity is still far from com- 
plete. 

(6) Diffusion of oxygen. That oxygen diffuses from the surface into the 
oxide follows from the same sort of evidence as that presented above for 
the diffusion of barium. The argument which was presented there and which 
was based on the life of the filament can be applied word for word for oxygen. 

Another line of evidence for the diffusion of oxygen is this: if current is 
drawn from a broken down coated filament even at as low a temperature as 
550°K the activity decreases. After the decrease amounts to a factor of 
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about three, the filament can be glowed at this same temperature. It is then 
found that the activity increases due to the glowing. This increase in activity 
is presumably due to the disappearance of oxygen from the surface. We should 
hardly expect oxygen to evaporate under these conditions: at any rate, we 
have never been able to detect it. Consequently we must conclude that the 
oxygen diffuses into the interior of the oxide. 

As more and more oxygen diffuses into the interior the concentration of 
oxygen in the oxide increases. Some of this oxygen presumably recombines 
with metallic barium; some of the oxygen must diffuse back to the surface 
at a rate which is presumably proportional to the concentration of oxygen 
in the oxide. If current is drawn from the filament for a long time at a con- 
stant temperature and plate potential, the net diffusion inward decreases 
and may eventually approach zero especially if the temperature is high 
enough so that oxygen evaporates. 

(7) Factors which determine the amount of electrolysis. Only a rela- 
tively small part of the current through the oxide can be carried by positive 
barium ions and negative oxygen ions. The ratio of ion current to total 
current for barium is not necessarily the same as for oxygen. For barium 
this ratio depends upon the temperature and the composition of the oxide; 
it is independent of the speed with which electrons from the outside strike 
the filament. For oxygen this ratio depends upon the composition of the oxide 
and the field strength at the surface of the filament when current is drawn 
from it; presumably it also varies with the temperature. 

The ratio barium ion current to total current can be determined from the 
data for curves like those in Figs. 5, 6 and 8. The calculations based on Fig. 
5 are given above: it was found that roughly one part in 200 of the current 
was Carried by strontium ions. A similar computation for Fig. 6 leads to the 
conclusion that in this case 1/240 of the current is carried by barium ions. 
For this calculation we have used the following data: at the initial activity 
the fraction of the surface covered with barium was 0.90 while at the optimum 
it was 1.00; this 0.1 layer was deposited in 36 mic-minutes or 48 mic-minutes 
per cm?; 2 mic-minutes per cm? of ion current were taken to correspond to one 
layer. A similar computation for Fig. 8 gave the result that about 1/300 of 
the current was carried by barium ions. 

A series of tests showed that when the temperature of the coated fila- 
ment was purposely raised above that prevailing in the tube while electrons 
were sent into it, a curve was obtained which was similar to that in Fig. 8 
except that the mic-minute axis was increased by a constant factor. From 
this we conclude that the ratio of barium ion current to total current de- 
creases by about a factor of 2.0 when the temperature is increased from about 
300°K to about 600°K. 

In another series of tests it was shown that the curve similar to that in 
Fig. 8 did not depend upon the potential of the oxide coated filament with 
respect to the source of the electrons. To obtain this result, it was however, 
necessary that the distribution of the current from one filament to the other 
was not altered when the potential between the filaments was changed. 
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Hence, the amount of electrolysis is not dependent upon the speed with which 
the electrons strike the oxide. 

Our knowledge of the factors that determine the ratio of oxygen ion cur- 
rent to total current is based partly on the factors that determine the rate 
of evaporation of oxygen which were discussed in a previous section. It was 
there shown that the evaporation rate depended greatly on the composition 
or state of the oxide and also on the plate potential. Presumably the rate of 
oxygen production varies in the same manner as the rate of evaporation. 

The rate of oxygen production was also studied by determining the de- 
crease in log i as a function of the mic-minutes drawn from the filament. 
The initial slope of a plot of these variables is a measure of the oxygen ion 
current to the total current. From a large number of tests it follows that this 
ratio decreases with the degree of breakdown of the filament. Presumably 
this means that as more and more metallic barium is dispersed through the 
oxide a larger and larger proportion of the current is carried by electrons. 

Another series of tests on about 10 different filaments has shown that the 
amount of oxygen produced per milliampere minute increases rapidly with 
the positive plate potential and that if the current is limited by space charge 
comparatively little oxygen is produced. This conclusion is again based on 
the initial slopes of log 7 vs. milliampere minute curves. 

It has been tacitly assumed that the rate of oxygen production varies 
directly with the electron current. This tacit assumption is difficult to prove 
since every factor which changes the emission current seems at the same time 
to change the rate of oxygen production. Since, however, it is shown above 
that the rate of barium production is directly proportional to the electron 
current, we can quite safely assume that the rate of oxygen production is 
also directly proportional to the current. 

The fact that the oxygen yield depends so markedly on the plate po- 
tential presumably means that part of the potential difference between the 
filament and plate is used up in the oxide but that the P.D. across the oxide 
is not directly proportional to the electron current through the oxide. In 
other words, the conductivity of the oxide does not obey Ohm’s law. From 
the fact that little or no electrolysis takes place when the current is limited 
by space charge, we are inclined to believe that there is little or no potential 
difference across the oxide in this case. The electrons get througk the oxide 
not because of a field in the oxide but because of a diffusion pressure similar 
to that which has been shown to exist in an arc. 

We have no evidence at present that the oxygen yield varies with the 
temperature but from analogy with the barium yield, it seems probable that 
the oxygen yield per milliampere minute decreases somewhat as the tempera- 
ture increases. 

(8) Time changes in coated filaments. A rather complete study has been 
made of the changes in emission that occur in time as a function of the 
temperature, the plate potential, the immediate previous treatment and the 
remote but violent previous treatment. This has been done for various types 
of filament. We shall first describe the time changes that occur in a barium 
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oxide+strontium oxide combined broken down filament. Then we shall 
discuss the changes in an uncombined filament, first, before breakdown and 
later after some breakdown has occurred. 

These time changes can best be described in terms of power emission 
curves. It was first found that the effect of any immediate previous history 
of the filament could be wiped out by a one to ten minute glowing at about 
1050°K. The activity atalow testing temperature could repeatedly be restored 
to a definite value by thistreatment. In the case of filament F; tube 12-J, this 
activity is represented by the uppermost dashed line A in Fig. 12. Since the 
emission decreased when current was drawn, it was in some cases necessary 
to extrapolate the current to zero time. At low temperatures the time changes 
are so slow that this need not bedone. A typical decay curve has been shown 
in curve 1 of Fig. 10 for a temperature of 715°K. If current is drawn for a 
sufficiently long time, these decay curves become practically flat. The emis- 
sion for this steady value on drawing current is represented by a point on 
line C Fig. 12 at the treating temperature. By glowing at 1050°K and draw- 
ing current at other deactivating temperatures until the emission became 
practically steady, other points on the C or steady value drawing current 
line were obtained. 

If the filament is glowed at the treating temperature after a point on the 
C line has been obtained, the emission increases rapidly at first then more and 
more slowly until it becomes practically flat. This is illustrated by curve 3 
in Fig.10. The emission for this flat porton of the curve is then plotted against 
the treating temperature, and aseries of such points at various temperatures 
constitute the D line in Fig. 12. 

The rate at which these time changes take place increases rapidly as the 
temperature is increased. Thus at a temperature of 600°K it might take some- 
thing like 10 hours to go from a point on the A line to a corresponding point 
on the C line, while at 900°K it would take only several minutes for the same 
transition. 

If at any stage in a deactivation curve (A-C) ora subsequent activation 
curve (C-D) the temperature is suddenly changed to a higher or to a lower 
value and if the instantaneous emission is obtained for this new tempera- 
ture, it will be found to lie on a line such as one of the dashed lines in Fig. 
12. These lines are called instantaneous lines. They are roughly parallel to 
one another. These instantaneous lines are the only ones which truly repre- 
sent the activity. 

The B line is obtained by glowing at 1050°K, reducing the temperature 
to a lower value for a considerable time, applying the plate voltage and ob- 
serving the emission at the instant the plate voltage is applied. If the plate 
voltage is kept on, the current will again end on the C line. 

With such a diagram of steady value and instantaneous lines it is possible 
to predict the time changes that will occur in a coated filament under any 
set of circumstances provided it is known where the state of the filament can 
be represented by a point on this diagram. For example, suppose we know 
that the state of a filament is such as to be represented by a point d; on the 
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D line at a low temperature 7;. Suppose we now wish to know what hap- 
pens if the temperature is increased to 7; and (a) no current is drawn or 
(b) current is drawn. In the first case the activity will be given by 7 in Fig. 
12 and will increase to d,. In the second case the current will decrease from 
12 tO Co. 

All the lines represented in Fig. 12 are for a single plate potential of 60 
volts. For a different plate potential an entirely new set of curves would 
hold. If the filament is given a severe treatment a complete new set of lines 
results. 

We shall next consider uncombined filaments. Fig. 13 shows two sets of 
power emission lines which show the type of time changes that occurred in 
an uncombined strontium oxide filament before it was broken down (lines 
with subscript 1) and after it was partially broken down (lines with subscript 
2). The dashed line A; represents the instantaneous emission current ob- 
tained at various temperatures after the filament had been thoroughly 
glowed at 1200°K. At any temperature below about 950°K the emission 
current increased with time, came to a peak value which is represented by 
a point on the P, line at that temperature, and then decreased to a steady 
value represented by a point on the B, line. Above about 950°K the current 
no longer passed through a peak but attained a steady value which again 
fell on the B, line. 

After large currents had been drawn from the filament for considerable 
periods of time so that the filament was partially broken down, the analo- 
gous instantaneous line Ae, the peak line P:, and the steady value line B, 
were obtained. At a given temperature the instantaneous activity had in- 
creased one hundred to five hundred-fold. The angles between A: and Pz, 
and A; and B, were now considerably smaller than the analogous angles in 
the preceding case. The point at which the peak phenomena disappeared 
was shifted to lower currents and to lower temperatures. These and other 
sets of such lines for different degrees of breakdown show that the peak phen- 
omena become less pronounced and are shifted to lower temperatures as the 
breakdown process continues. The angle between the instantaneous line 
and the steady value drawing current line decreases as the breakdown pro- 
cess increases. These deductions fit in well with Fig. 12 for which no peak 
phenomena were observed and for which the angle between the instantane- 
ous line A and the steady value glowing line C is very small. 

From the discussions which have been given in preceding sections it 
follows that the only lines whose slopes may properly be used to obtain work 
functions are the instantaneous lines. If the data represented by the A; 
and A; lines in Fig. 13 and the A line in Fig 12 are replotted in the customary 
way, i.e., log 7-2 log T vs. 1/T, straight lines are obtained whose slopes yield 
work functions of 1.8, 1.3, and 0.53 volts respectively. Since it is question- 
able that the work function of a composite surface is independent of the 
temperature, no great significance is to be put on any work functions deter- 
mined in this manner. As far as they go, however, the values for the work 
functions are quite satisfactory, since the value 0.53 volts is for a filament 
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which was in an unusually active condition. It might be well to point out 
that if work functions had been based on the steady value drawing current 
lines there might have been an error by more than a factor of two. 

From the discussion of time changes given in preceding sections, the 
reader can appreciate the significance of these various lines. Their complete 
interpretation must be left until an even more exhaustive study has been 
made. 


SUMMARY 


The activity of a coated filament depends upon the surface concentration 
of metallic barium, the surface concentration of oxygen, on the relative posi- 
tions of the barium and oxygen, and on the composition of the underlying 
oxide. Barium can be deposited on the surface either by evaporation from 
an outside source, which may be another coated filament, or else by electroly- 
sis due to sending electrons into the filament from an outside source. In 
either case the barium on the barium oxide acts similar to barium on tungs- 
ten: as the amount of barium increases, the log of the emission at a given 
temperature increases rapidly at first and then more and more slowly until 
it passes through an optimum value beyond which it decreases. Oxygen may 
also be deposited on the surface either by evaporation or by electrolysis due 
to drawing current from the filament itself. The effect of this oxygen on the 
activity is quite complex: in general it is similar to the effect of oxygen in 
the case of barium on oxygen on tungsten; if the oxygen is between the sur- 
face layer of barium and the oxide, it increases the activity while if the oxygen 
is above the surface layer of barium it decreases the activity. Both barium 
and oxygen diffuse quite readily from the surface into the oxide and vice 
versa. The rate of diffusion increases rapidly with the temperature and with 
the concentration gradient near the surface. When current is sent into the 
oxide or drawn from the oxide most of-it is carried by electrons, but a small 
portion is carried by oxygen and barium ions. The ratio of ion current to 
total current depends upon the temperature, the composition of the oxide, 
and in the case of drawing current on the plate potential. The complex 
changes in activity that occur in time when the filament is treated in various 
ways can best be described in terms of various “power-emission” or equiva- 
lent Richardson lines, namely, instantaneous lines,steady value glowing lines, 
steady value drawing current lines, and peak lines. These complex time 
changes can in a general way be interpreted in terms of the proposed theory. 
Some of the features require further study. 
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ABSTRACT 


Method for eliminating secondary electrons in the study of the voltage-intensity 
relations in spectra—All previous measurements of the relation between the intensity 
of spectrum lines and the velocity of the exciting electrons have been subject to uncer- 
tainties due to the presence of secondary electrons. To eliminate these effects a beam 
of electrons is so concentrated by a magnetic field as to prevent their striking any metal 
surfaces until they reach a collector which traps all secondaries. This arrange- 
ment insures that the light emitted from the path of such a beam is produced solely by 
electrons which have a specified velocity (to within the Maxwellian distribution). 

Voltage-intensity cha:acteristics for Hg 2537 up to 7.4 volts —The above method 
was applied to a study of the variation with electron velocity of the intensity of \2537 
in the Hg spectrum. The results obtained agree in a general way with those of others. 
A maximum intensity is found at about 6.0 volts which agrees within 0.2 volts of the 
values given by White and by Hanle. The rate at which the intensity drops off beyond 
the maximum is, however, considerably greater than that given by previous workers. 
This may be attributed to the elimination, in the present work, of secondary elec- 
trons. 


INTRODUCTION 


OR a more complete understanding of the dynamics of the collision 

process between electrons and atoms, a knowledge of the probability of 
excitation of the atom to a given energy state as a function of electron energy 
is necessary. There are two experimental methods of attacking this problem: 
(1) from a study of the number of electrons which, at any given speed, lose 
the energy corresponding to the given energy state of the atom; and (2) from 
a study of the intensity of the spectrum lines which are emitted as a result 
of the excitation process. The latter method is practically limited to the 
study of the efficiency of excitation to the lowest excited state cf the atom 
and must be carried out under conditions which do not permit of anything 
happening to the atom between the excitation and the radiation processes. 

Unfortunately an examination of the existing literature shows a rather 
wide discrepancy, not only in the results as obtained by the two methods, 
but also in the results of different observers who have used the same method. 

One source of error in practically all spectroscopic studies of excitation 
efficiencies has been the presence, in the primary beam of electrons, of large 
numbers of secondary electrons of undetermined speed. In principle, the 
apparatus almost invariably used! in such studies is illustrated in Fig. 1. 

1D. R. White, Phys. Rev. 28, 1125 (1926); P. Bricout, J. de Physique et le Radium 
9, 88 (1928); W. Hanle, Zeits. f. Physik 54, 848 (1929); W. Hanle, Naturwissenschaften 
15, 832 (1927); W. D. Crozier, Phys. Rev. 31, 800 (1928); J. Valasek, Phys. Rev. 29, 817 
(1927); J. H. Lees, H. W. B. Skinner, Nature 123, 836 (1929); M. G Peteri, W. Elenbaas, 
Zeits. f. Physik 54, 92 (1929). 
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It consists of a source of electrons F, a grid G to which is applied the variable 
accelerating potential, and a plate P which serves to collect the electrons 
which pass through the grid. The region of spectroscopic observation is the 
field-free space between G and P, and it is usually assumed that the electrons 
in this region have a velocity corresponding to the potential applied to G 
(with proper corrections for contact e.m.f., etc.). No one who has studied 
the velocity distribution of the electrons in such an arrangement can fail 
to have been impressed by the large number of electrons which have speeds 
less than that corresponding to the applied potential. These slower electrons 


G P 


F 
J 





Fig. 1. Arrangement of electrodes customarily used in a study of the excitation 
function of spectrum lines. 


are presumably secondary electrons which come from P and G by reflection 
or otherwise. Their number varies greatly with the nature and geometry 
of the electrodes. The importance of correcting for their presence has been 
realized by others but, in view of the uncertainties involved, it is desirable 
to devise a method which eliminates them entirely from the region of obSer- 
vation. 

In what follows such a method will be described and its use will be illus- 
trated by some preliminary measurements on the voltage-intensity relations 
for the 2537A line of the mercury spectrum. 


APPARATUS AND METHOD 


The apparatus is sketched in Fig. 2. The boxes A, B, C, are made of 
molybdenum. To support them in rigid alignment, short tungsten rods 
fastened to the boxes are sealed intoa glass frame. The entire assembly is 
next inserted in a Pyrex tube which is then closed and attached to the vacuum 
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Fig. 2 Diagram of apparatus used in the present study. 


system. Box A is maintained at the potential of the mid-point of the fila- 
ment. A magnetic field of 200 gauss from two Helmholtz coils mounted 
with their axes parallel to the axis of the boxes prevents any lateral spread 
of the electron beam. The diaphragms in boxes B and C are of such size that 
no electrons from the filament will strike them. The box C contains a plate 
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parallel to the electron beam but removed from it by 5 mm. This plate is 
maintained at a positive potential of 1200 volts with respect to box C and 
serves to remove electrons from the beam and prevent the return of secon- 
daries to B. The electrons are accelerated from the filament F to the box B. 
The space within B is field-free. The mean free path of the electrons should 
be about ten times the distance from the filament to the collector. The 
window W in box B is one of two openings opposite each other. One of the 
openings faces a horn in the glass tube which is painted on the outside with 
black enamel and serves as a light well. The light to be studied passes out 
through the opening W and through a quartz window to whatever photo- 
metric device may be chosen for the particular purpose in view. 


VOLTAGE-INTENSITY RELATION FOR 2537A IN MERCURY VAPOR 


To illustrate the applicability of the method, it was applied to the study 
of \2537 of Hg. For this purpose the pressure of Hg was limited by the 
presence of ice on the trap. The voltage was varied by steps of 0.2 volts 
from 4.6 to 7.0 volts, except near the maximum, where a closer study was 
made. In some of the work the filament was placed farther in front of box A 
than indicated in the diagram. This necessitated a field between box A and 
the mid-point of the filament. The results obtained with both filament ar- 
rangements gave essentially the same curve. 

A quartz lens focused the light on a photographic plate. Eastman 40 
plates were used. A special plate holder made possible an entire run on one 
plate. On the same plate graded exposures to the light were made by inter- 
posing in the light path metal screens which had been calibrated with a 
thermopile to determine the percentage radiant energy transmitted. These 
exposures were taken with the applied voltage which gave the maximum 
intensity, and this enabled us to determine the relation between the blacken- 
ing on the photographic plate as read on the microphotometer and the 
energy delivered to the plate. Since all pictures were taken with 90 seconds 
exposure, there was no need to assume the validity of the reciprocity law for 
time and intensity of radiation. The plates yielded very nearly a straight- 
line relation between the logarithm of energy and the deflection on the 
microphotometer. Photographs taken at 4.1 applied volts showed no stray 
light. The current to the plate was measured and was about 27 micro- 
amperes. Unfortunately because of imperfect alignment there was a slight 
current to the front face of box B. 

The results of the measurements are plotted in Fig. 3. Their reproduci- 
bility can be estimated from a comparison of the positions of the circles and 
dots which represent two independent sets of observations. In the figure the 
voltage scale is arbitrarily shifted so that the curve, extrapolated to zero 
intensity, will cut the axis at 4.86 volts. 

The chief points of interest in the curve are the steep initial rise and the 
abrupt decline beyond the maximum, together with the indication that the 
curve is a composite of two curves, one with a maximum at about 5.5 volts 
and the other with a maximum at about 5.95 volts. Since the whole radiation 
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from the beam was focused on the photographic plate it is not possible to 
decide whether these two maxima are characteristic of \2537 alone or whether 
at about 5.7 volts a new radiation, possibly a molecular band, is added to 
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Fig. 3. Variation with electron energy of the intensity of Hg 2537. 


2537. In this connection it may be significant that in critical potential 
measurements under similar conditions a critical potential at 5.7 volts is 
always found. 
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Fig. 4. Comparison of present curve (O) with the results of others. B, Bricout; C, 
Crozier; H, Hanle; W, White (see reference 1). 


In Fig. 4 the present curve is plotted together with all curves previously 
published on \2537. The most prominent distinction of the present curve is 
the much more rapid decline from the maximum. The reason for this is 
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undoubtedly the elimination of secondary electrons which in this region 
would be particularly effective in increasing the amount of radiation. 

A detailed comparison of the present curve with the others is impossible 
because, in addition to the elimination of secondary electrons, the vapor 
pressure was here very much lower than in the other cases. The experi- 
mental conditions in Hanle’s curve are probably most nearly like those of 
the present work and the differences in the curves are probably to be as- 
cribed largely to the elimination of secondary electrons. In White’s curve 
which most nearly agrees with the present one in position of maximum, 
the vapor pressure was so high (saturated vapor pressure at about 100°C) 
that most of the fine structure is attributed by White to secondary effects. 
It is difficult to reconcile the present curve with those of Bricout or of Crozier. 


SOURCES OF ERROR IN THE PRESENT WORK 

Space charge made it impossible to take a velocity distribution. The 
voltage scale is uncertain, but is arbitrarily extrapolated to give zero inten- 
sity at 4.86 volts. Space charge might tend to compress or extend the applied 
voltage scale. A rough calculation has shown this to be 0.1 volt or less for 
the range of voltages studied. Space charge is known not to make the velocity 
distribution any broader than usual. This is also born out by the sharpness 
of the curve. Space charge probably causes a velocity gradient in the beam. 
However, in this work it was possible to study a portion of the beam 0.15 cm 
in length. Velocity gradient could not cause much loss of homogeneity in 
this length of beam. 

Electrons spiral about the magnetic lines of force and make the path in 
mercury vapor uncertain. This is due only to the component of velocity 
perpendicular to the path of the beam. This should not be more than a few 
tenths of a volt. The only effect on the shape of the curve would be the 
difference in the perpendicular component for different velocity electrons. 
This is a second order correction. In addition to this, the radius of the 
spiral is proportional to the square root of the voltage. This divides by two 
any percent of error due to this source. Furthermore the ratio of the per- 
pendicular component of velocity to the total should be practically the same 
over this range of velocities. 

It might be supposed that the radiation given out by atoms excited by 
electrons of 6.0 or 7.0 volts would have a slightly different fine-structure 
from that of the light given out by atoms excited by electrons of 5.0 volts. 
This might cause the absorption coefficient of the mercury vapor to change 
considerably, and would mean that the observed curve would not represent 
the true intensity. Hanle tested this by putting additional mercury vapor 
in the path of the light and taking another test. He concludes that the shape 
of the curve is, within experimental error, independent of the length of path 
in mercury vapor. 

In conclusion, I wish to thank Professor John T. Tate who suggested the 
problem and its method of attack and under whose guidance the work was 
done. Thanks are also due to the entire staff, especially Mr. Wm. Haliday, 
Mr. Charles Johnson, and Mr. E. Greinke for help with the apparatus and 
the glass work. 











NOVEMBER 15, 1929 PHYSICAL REVIEW VOLUME 34 


THE MAGNETIC MOMENT OF THE OXYGEN ATOM 


By O. E. Kurt! wits T. E. Puipps 
LABORATORY OF PHYSICAL CHEMISTRY, UNIVERSITY OF ILLINOIS 


(Received October 2, 1929). 


ABSTRACT 


The magnetic moment of the oxygen atom has been investigated by the atomic 
ray method of Stern and Gerlach. The atomic oxygen was formed by means of an elec- 
trodeless ring discharge and the ray was detected by its oxidizing action on freshly sub- 
limed litharge. When the ray was passed through an inhomogeneous magnetic field the 
resulting image consisted of an undeviated central line and two broad lines displaced 
symmetrically, one on each side. Hydrogen was substituted for oxygen and the pro- 
cedure was repeated. The separations obtained in the two cases provided the data 
necessary for the calculation of the effective magnetic moment of the oxygen atom. 
The experimental magneton values were zero and +1.67. The magneton numbers ex- 
pected from the spectroscopic ground term of oxygen are 0, +3/2, and +3. The 
effective magneton numbers which (due to overlapping of the images) would result from 
these values are zero and +1.71. The agreement is within the limits of experimental 
error. 


INTRODUCTION 


HE experimental method employed in this investigation is that develop- 
ed by Stern and Gerlach.? Rays of atoms possessing a magnetic moment 
are deflected into separate rays when passed through an inhomogeneous mag- 
netic field at right angles to the direction of the field. The deflection s, ex- 
perienced by those atoms in the ray having the most probable velocity is 


Sq=1/4-M/RT-dH/ds:(Ly2+2L,L2) (1) 


where M is the magnetic moment in gauss-cm per mol; L; is the length of 
the pole pieces of the magnet, or the distance the atoms must travel through 
the magnetic field; Lz is the distance from the end of the field to the target; 
dH/ds is the value of the inhomogeneity of the field in the path traversed by 
the atomic ray. In this equation 0H/ds is assumed constant throughout the 
path of the ray in the field. If the magnetic moment for one element is known, 
then the magnetic moment for another element can be calculated from the 
relative deflections obtained under identical conditions in the two cases. 
This method obviates the tedious mapping of the field necessaryinthe evalua- 
tion of dH/ds. The equation employed is 


M2= My: 5q,T2/Sa,71 (2) 


where M, and M; may be expressed either in gauss-cm or in Bohr magnetons. 
In this investigation the hydrogen atom, assumed to have a magnetic mo- 


! This communication is an abstract of a thesis submitted in partial fulfilment of the re- 
quirements for the degree of Doctor of Philosophy in Chemistry at the University of Illinois. 
? Stern, Zeits. f. Physik 7, 249 (1921); Gerlach and Stern, Zeits. f. Physik 8, 110 (1922); 
9, 349, 353 (1922), Ann. d. Physik 74, 673 (1924); 76, 163 (1925). 
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ment of one Bohr magneton, is used as the reference substance in Eq. (2). 
Magneton values for the oxygen atom have been determined and compared 
with theoretical predictions. 


EXPERIMENTAL 


The electrodeless ring discharge. A long discharge tube of the R.W.Wood 
type with aluminum electrodes was used in the early attempts to prepare 
atomic oxygen. This method was found unsatisfactory since under optimum 
conditions for the discharge the aluminum electrodes were oxidized with a 
considerable evolution of heat, and a consequent decrease in the pressure 
due to the “clean-up” of the active species. 


a 
L ) ra 






































Fig. 1. Diagrammatic sketch of apparatus. 


The electrodeless ring discharge was found to be a more satisfactory 
method. The electrical circuit is shown in Fig. 1. The secondary of a 1 k.w. 
25000 volt transformer was connected with stranded leads of No. 8 B. & S. 
insulated copper wire through a condenser C of about 0.015 uf capacity, to 
six turns of copper tubing, about 0.3 cm in outer diameter, spaced equally 
around a Pyrex bulb of 200 cc capacity. The condenser consisted of 15 
plates, each having an area of 500 cm*. Two thicknesses of double-strength 
window glass between plates served as a dielectric, and the entire condenser 
was immersed in transformer oil. The spark gap SG consisted of two zinc 
electrodes, each having an area of 2 cm?, placed about 1.5 cm apart. 

In the preliminary experiments precautions were taken to insure the pur- 
ity and dryness of the oxygen admitted to the discharge bulb. However, 
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ordinary commercial tank oxygen was found satisfactory and was used in 
all later experiments. The oxygen was admitted through a glass capillary, 
adjusted to give a suitable pressure in the discharge bulb. 

The discharge was found to pass at any pressure between 0.007 and 0.45 
mm. After experience had shown that the greatest chemical activity was 
found at a pressure of about 0.13 mm, this pressure was adopted in all later 
experiments. An intense line spectrum against a black background was ob- 
served under these conditions by means of a Zeiss pocket spectroscope. As 
seen by the eye in daylight, the discharge consisted of a general glow filling 
the bulb, with streamers of a different color following the directions of the 
coil. At low pressures the glow was a blue-gray while at high pressures a 
distinct pink color appeared in the blue. The streamers in both cases had 
a greenish blue appearance. Impurities such as mercury vapor and nitrogen 
were detected by their characteristic lines in the spectrum and by the changed 
appearance of the glow. Heating the walls of the bulb or running the dis- 
charge continuously for several hours had the effect of removing the impuri- 
ties and intensifying the oxygen spectrum. The effect of water vapor as an 
impurity was unique and will be discussed at length in a later paragraph. 

The detection of atomic oxygen. Freshly sublimed litharge was found to be 
an excellent detector at room temperature of atomic oxygen. The presence 
of water vapor, even in small traces, intensified the oxidizing conditions in 
the discharge. An extended investigation into the effect of varying amounts 
of water vapor upon the discharge led to the results summarized below. 

Mixtures of oxygen and water vapor containing as much as 40 mol per- 
cent water vapor oxidized litharge in a few seconds to brown lead dioxide, 
when in direct contact with the discharge. Lead dioxide could be distin- 
guished by chemical test from litharge and metallic lead, since it is very 
slightly soluble in 25 percent nitric acid, and gives a green color with an 
acetic acid solution of diphenylamine. If water vapor were in excess of about 
40 mol percent the action was distinctly reducing, as in a hydrogen discharge; 
lead dioxide was reduced to litharge and litharge to metallic lead under these 
conditions. It was found that lead thus formed would be reoxidized directly 
to lead dioxide without the intermediate formation of litharge, when the con- 
ditions in the bulb became oxidizing. 

The effect with dry hydrogen, with hydrogen saturated with water vapor 
at room temperature, or with water vapor alone in the discharge bulb was 
the stepwise reduction of lead dioxide to litharge, and of litharge to metallic 
lead, when these were in direct contact with the discharge. However, ata 
distance from the discharge, as under molecular ray conditions, lead dioxide 
was reduced to litharge but litharge was not reduced to metallic lead. This 
peculiar result has been noticed by another investigator.* Not only did the 
chemical effects indicate the presence of atomic hydrogen in all of these cases, 
but the Balmer series was observed to be very intense, and small particles 
on the walls of the discharge bulb glowed brightly. Both of these phenomena 
are characteristic of hydrogen atoms, and were never observed when dry 


? DeHemptinne, Ann. Soc. Sci. Bruxelles 47B, 143 (1927). 
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oxygen was admitted toa discharge bulb whose walls had been well outgassed. 
A very dry condition however, was not as conducive to intense oxidizing 
effects as was the presence of even a trace of water vapor. Hence in subse- 
quent experiments the oxygen was saturated with water vapor at room tem- 
peratures before being introduced into the discharge bulb. 

The detection of atomic oxygen with blue molybdenum oxide was also 
investigated. Since it was necessary to heat the oxide to a temperature of 
200°C before oxidation by oxygen atoms occurred, it was considered less 
satisfactory than litharge. 

The slit system. To produce a ray of atoms, an all Pyrex slit system was 
sealed to the discharge bulb as shown in Fig. 1. This consisted of two slits, 
0.12 cm long and 0.005 cm wide, aligned with each other at a distance of 
6.6 cm.‘ The space between the slits was connected with tubing 3.0 cm in 
diameter through a liquid air trap to a high speed mercury vapor pump. 
The second slit was connected to the target chamber by means of the small 
tube D. The target chamber was connected through a liquid air trap to a 
second high speed mercury vapor pump. 

The target. The target consisted of a microscope cover glass held in an 
all Pyrex clamp so designed that the cover glass could be easily removed and 
mounted, for convenience in handling and photographing, on a microscope 
slide. No sealing wax or cement was used in holding the target. Repeated 
failures to obtain images when sealing wax or de Khotinsky cement were 
present in the vacuum led finally to the conclusion that these substances had 
a desensitizing effect upon the litharge surface. (This effect was not noticed 
with molybdenum oxide surfaces exposed to rays of hydrogen atoms.) The 
litharge was active only when freshly prepared and kept in a vacuum. 
Even in a vacuum it became less sensitive if the walls of the target chamber 
were heated, for the purpose of outgassing. 

The litharge was deposited upon the microscope cover glass by sublima- 
tion. An oxyhydrogen flame was directed against a surface of molten lith- 
arge, thus producing a dense cloud of very fine particles. In order to obtain 
a fine-grained surface the cover glass was held close to the molten litharge. 
Under the best conditions surfaces were produced which showed no grain 
under high magnification. A transparent pale yellow layer about 4 10-* cm 
in thickness was produced in this manner. 

The magnetic field. A large DuBois magnet, an exact replica of the mag- 
net used in a previous investigation,’ was used in this research. The pole 
pieces were similar to those employed by Leu® and in accordance with his 
observations the inhomogeneity of the field may be assumed to be uniform 
throughout the region traversed by the ray. A cross-sectional view of the 
pole pieces is shown in Fig. 2. These were 6 cm long and were placed 0.3 cm 
apart. The slot was 0.6 cm wide. The edge of the 60° wedge-shaped piece 


‘ The method of constructing these slit systems was devised by one of the authors, and 
will be published later. 
5 Phipps and Taylor, Phys. Rev. 29, 308 (1927). 
* Leu, Zeits. f. Physik 41, 551 (1927). 
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was 0.1 cm wide. The position of the small Pyrex tube D and the path of 
the ray between the pole pieces are shown in Fig. 2. 

Procedure. The apparatus assembled for a run is shown in Fig. 1. The 
discharge bulb was placed above the magnet, since in this position the dis- 
charge suffered less disturbance from the magnetic field than in any other 
position. Even in this position the current in the magnet coils could not 
exceed 5 amperes without causing serious interference with the discharge. 
The leak LZ (Fig. 1) was adjusted to give a pressure of 0.13 mm with oxygen 
in the discharge bulb. The pressure was measured by means of a McLeod 
gauge attached through a liquid air trap to the system at the point A. 
Before passing the discharge it was necessary to seal off the McLeod gauge 
to avoid “clean-up” of the atomic oxygen by the mercury of the gauge. 

A runwas then made with oxygen. The current in the primary of the trans- 
former was about 13 amperes and that in the magnet coils was about 3.5 
amperes. The re-entrant tube ~ in the discharge chamber permitted the 
insertion of a thermometer or thermocouple into the discharge for measuring 
the temperature near the first slit. The temperature remained fairly steady 
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Fig. 2. Diagrammatic shetch of pole pieces (cross-section). 


at about 340°C during the run. A.thermometer was employed since a ther- 
mocouple was found unsatisfactory due to a disturbing effect of the discharge 
upon it. 

An undeviated central line appeared after about two hours as a brown 
streak of lead dioxide on a pale yellow litharge background. After about 
nine hours faint lines on each side of the central line appeared. The run was 
continued for twenty hours,’ after which the target was removed and attached 
to a microscope slide. In order to facilitate measurements of the image, 
photomicrographs of the image (7X) were taken (Fig. 3). It was found that 
blue transmitted light gave the best photograph of the oxygen image. 
Cramer’s Super-contrast plates were employed. A density study of the re- 


7In a preliminary run without the magnet, with 20 amperes in the primary circuit, 
images were visible in the vacuum within fifteen minutes. In the run with the magnet de- 
scribed above it was necessary to use a smaller current for a longer time due to the low capacity 
of the transformer. It was at first a source of surprise that twenty hours of running should be 
necessary to bring out the side lines to an intensity which made photographing possible. This 
is probably due (1) to the fact that the magnet current reduces the intensity of the dis- 
charge considerably, (2) to the fact that the magnetic field broadens the deflected ray 
and thus reduces its intensity, and (3) to the inadequacy of the transformer for continuous 
duty with large primary currents. Beyond this, however, it is probable that some unknown 
condition was unfavorable in this particular run. 
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sulting photographic negative was then made with a Moll Recording Mic- 
rophotometer (Fig. 5a). A microphotometer study of the original image 
was unsuccessful, probably due to the slight difference in transmission coef- 
ficients of litharge and lead dioxide for the longer wave-lengths. 





Fig. 3. Photomicrograph of Fig. 4. Photomicrograph of 
oxygen image (7X). hydrogen image (7X). 


The method employed in this investigation involves a comparison under 
identical conditions of the relative deflections of rays of oxygen and hydro- 
gen atoms. The next step then was the substitution of hydrogen for oxygen 





Fig. 5a. Microphotometer record of negative of oxygen photomicrograph. 


in the discharge bulb. Molybdenum trioxide was used as a detector. The 
current in the primary of the transformer was 16 and that in the magnet 
coils was 3.5 amperes. The temperature of the discharge as recorded by the 
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thermometer near the first slit was 328°C. After fifty minutes the double- 
line splitting characteristic of the hydrogen atom appeared. The run was 
continued for three hours, after which the target was removed from the 
vacuum and photographed (Fig. 4). 

A ppearance of the images. Fig. 3 shows the photomicrograph of the oxygen 
image taken by transmitted blue light. From the image it can be seen that 
the inner edge of the deflected ray is relatively sharp compared with the 
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Fig. 5b. Smooth curve constructed from Fig. 5a. 


outer edge, where the density falls off very gradually. This condition is 
unfavorable to the direct determination of points of equal density in the 
deflected ray, a procedure necessary to the calculation of the results. For 
this reason the microphotometer record (Fig. 5a) was made from the nega- 
tive. Measurements were then made from Fig. 5b, which is a smooth curve 
constructed from Fig. 5a. The unsymmetrical appearance of the micro- 
photometer curve for the two side lines in Fig. 5 is due to the fact that the 














Fig. 6a. Microphotometer record of negative of hydrogen photomicrograph. 


value of 0H/ds is not constant over the entire region traversed by the ray. 
served by means of light scattered along the X axis. 

In the case of the hydrogen image also it was found advantageous to 
make the measurements from a microphotometer record of the negative. 
The density curve and the smooth curve constructed from it are given in 
Figs. 6a and 6b. 
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Since the value of 0H/ds is more constant on the side of the ray toward 
the slotted pole piece, all measurements were made from the corresponding 
side (left side) of the microphotometer curves. 


—— [Experimental 
——-— Theoretical 
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Fig. 6b. Smooth curve constructed from Fig. 6a. 


RESULTS AND CALCULATION OF MAGNETON VALUES 


In order to calculate the relative magnetic moment of the oxygen atom 
from Eq. (1) it is necessary to know the value of s,. for both hydrogen and 
oxygen. This value however cannot be measured directly from the micro- 
photographs of the images but must be calculated by means of an equation 
derived by Stern® from a consideration of the width of the undeflected ray 
and the Maxwellian distribution of velocities of the atoms in the ray. This 
equation is 


Sq = 35152/(S2—51) [In(s2/s1) +a?/s12{ 2— 454/351 +522/651} | (3) 


where s; and sz are points on the inner and outer edge, respectively, of the 
deflected streak, which have the same intensity; a is one-half of the width 
of the undeflected ray. Values of s; and sz for oxygen were obtained from 


TABLE I 

Dimensions of Apparatus 
nt DN A ng ck cc ccncvneicsseeascseweses 6.6 cm 
Distance from the second slit to the target................. 9.6cm 
ER ik cance decreas cderreadecens 6.0 cm 
Distance from the end of the field to the target, Lo.......... 2.3 cm 
reer Pere 0.0045 cm 
EE EL ee 0.005 cm 
RE IR RR OT 0.12 cm 
One-half the width of the undeflected line, a................ 0.0060 cm 


Measurements from Microphotometer Records 











T St S2 Sa M 
Substance °K m m m Bohr magnetons 
Hydrogen 601 100 214 383 
129 179 405 
Avg. 394 1 
Oxygen 613 Experimental 171 300 636 
186 286 649 
Avg. 642 1.67 
Theoretical 175 310 658 
190 290 666 
Avg. 662 1.71 





* Stern, Zeits. f. Physik 41, 563 (1927). 
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measurements upon Fig. 5b, and for hydrogen were obtained from Fig. 6b. 
These values are given in Table I. 

From these values of s; and sz corresponding values of s_ for oxygen and 
hydrogen were calculated by means of Eq. (3). These values were then 
substituted in Eq. (2) to obtain the average value of + 1.67 Bohr magnetons 
for the relative effective magnetic moment of the oxygen atom,—besides 
the value of zero. 

Discussion of errors. The maximum probable error of s, for oxygen is 
about three percent, and for hydrogen about four percent, and for each 
of the values of J about one percent. The value of 0H/Qs is sufficiently 
uniform throughout the region traversed by the ray so that any error intro- 
duced by the use of Eq. (1) instead of a more exact equation is negligible. 

Theoretical prediction of the magnetic moment. The ground term of the 
oxygen atom has been calculated by Hund?® from a consideration of atomic 
structure to be a *P triplet state. This prediction has been confirmed by 
spectroscopic data.'!® The atom has one stable state, the *P, state, and two 
metastable states, the *P; and *P, states. Since the term values for the 
triplet level are known, the relative probability of each state can be calcu- 
lated from the equation 


W= (2j+1)eF/*T = (2j+1)emMrel #7, (4) 
The values of the probability W calculated for each state appear in column 


4, Table II. The magneton values in column 5 have been calculated according 
to the usual method? for each of the P terms. 


TABLE II. Probabilities of the normal states. 











state term values!® Av W magneton values 
3P, 109833 cm7 0 5 +3, +3/2, 0 
sP, 109674 159 2.07 +3/2, 0 
5Po 109607 226 0.59 0 








By summing up the probabilities for the magneton values in the separate 
states the total probability for each of the values is obtained. These are 
given in Table III. They are used later in calculating the resultant effective 
magnetic moment of the oxygen atom. 


TABLE III. Probabilities of the magneton values. 











magneton values W 
+3 or —3 0.438 
+3/2 or —3/2 0.741 
0 1 








The theoretical intensity distribution curve for hydrogen is presented 
in Fig. 7, based upon the value of s given in Table I. The equation of Stern® 
sa/s+a 
I/Jo=3[e-“(y+1) Jeera (5) 
® Hund, “Linienspektren” (1928). 


1° Hopfield, Nature 112, 437 (1923); Phys. Rev. 21, 710 (1923); Astrophys. J. 59, 114 
(1924). 
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was employed, in which J is the intensity of the deflected ray at a distance 
s from the center of the undeflected ray; s, and a are the same quantities 
which appear in Eq. (3); and Jp is the average intensity of the undeflected 
ray. 

From the value of s, for the hydrogen atom (M=1, T=601°K) were 
calculated values of s, for the oxygen atom (M= +3/2 and +3, T=613°K). 
Similar intensity distribution curves were plotted for these values (curves 
B and C in Fig. 7). Curves B and C are intensity distribution curves re- 
spectively for M=3/2, and M=3 with the proper statistical weights given 
to the values of J/Jo, shown in Table III. Curve D which is a summation 
curve from curves B and C represents the intensity distribution which should 





J 
4 


OM 





rh 
010 4 am 


] / ry a 
\ 
005 7 \ 





























“™., 
~ 
att. —_— 
” it nn = 
0 Tt — 
200 400 600 800 sna 


Fig. 7. Intensity distribution curves calculated as follows 
A for M=1, s=394u. B for M=3/2, s=580,. 
C for M=3, s=1160u. D is a summation of B and C. 


obtain for the oxygen image. This curve appears also in Fig. 5b, for conven- 
ience in comparing it with the experimental intensity distribution curve. 
The theoretical curve is not expected to coincide with the experimental 
density distribution curve since the latter may be distorted in the direction 
of the intensity axis. However values of s, and hence of M calculated from 
corresponding values of s; and s, taken from the two curves should agree. 
The extent of the agreement can be seen from a comparison of Sq (exp.) 
with s, (theoret.) for comparable values of s; or sz, as given in Table I, or 
from a comparison of the two curves in Fig. 5b. The experimental values 
of the effective magnetic moment of the oxygen atom, zero and +1.67 agree 
within the limits of experimental error with the theoretical values zero and 
+1.71 Bohr magnetons. 

In conclusion the writers wish to thank Professors A. P. Carman and 
E. H. Williams of the Department of Physics for their kindness in allowing 
the use of the large Dubois magnet, and Mr. H. N. Swenson of the Depart- 
ment of Physics for his aid in making the microphotometer records. 
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EFFECT IN HYDROGEN-LIKE ATOMS 


By Boris PopoL_sky* AND VLADIMIR ROJANSKY 
UNIVERSITY OF CALIFORNIA, BERKELEY, AND WASHINGTON UNIVERSITY, St. Louis 


(Received October 11, 1929) 


ABSTRACT 


It has been shown by Schrédinger and O. Klein that the Kramers-Heisenberg 
formulas for dispersion and incoherent scattering can both be obtained by an ex- 
tension of Schrédinger’s method for treating dispersion, the terms in the Hamiltonian 
involving squares of the potentials being neglected. The resulting formulas agree with 
those of Dirac. They are, however, inconvenient for calculations, as they contain 
summations with respect to all energy levels combining with the pair of levels under 
consideration, and thus imply complicated integrations when a continuous spectrum is 
present. This paper treats the problem of a hydrogen-like atom, acted on by light the 
wave-length of which is large compared to the size of the atom and the frequency of 
which is not too near a resonance frequency, and develops a method which obviates the 
necessity of integrating over the continuous range. It is an extension of the method 
used by one of us in treating the dispersion by atomic hydrogen. It is applied in detail 
to the first two levels of H. Formulas are derived for the intensities of the Smekal- 
Raman lines with a shift corresponding to the first absorption line. 


T HAS been pointed out by several investigators that the Smekal-Raman 
effect is to be expected from the formula of Kramers and Heisenberg,’ 
which has also been derived entirely on the new mechanics by Dirac,? Schré- 
dinger® and O. Klein.‘ This formula, when applied to a system possessing 
both a discrete and a continuous spectrum, involves a summation over the 
discrete range and an integration over the continuous range, and, due to the 
complexity of the integration, is inconvenient for the actual calculation 
of the intensities of the expected lines. In the present paper a theory is 
developed for handling the case of a hydrogen-like atom which yields in- 
tensity formulas involving only summations, and thus permits the numerical 
evaluation of intensities with relative ease. The width of the spectral lines 
is neglected in the discussion, and the final formulas involve the resonance 
denominators characteristic of the formula of Kramers and Heisenberg. 
We start with the usual Hamiltonian function for a hydrogen-like atom 
in an electromagnetic field’ in which the effects of relativity and electron 
spin are neglected: 


H=1/2y-2(piteAi/c)?—ed, i=1,2,3 (1) 


* National Research Fellow. 

1H. A. Kramers and W. Heisenberg, Zeits. f. Physik 31, 681 (1925). 

2 P. A. M. Dirac, Proc. Roy. Soc. 114A, 710 (1927). 

+E. Schrédinger, Ann. d. Physik 81, 109 (1926). 

‘ O. Klein, Zeits. f. Physik, 41, 407 (1927). 

5 See for example H. Weyl, Gruppentheorie und Quantenmechanik, p. 88 (1928). 
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where yu is the mass of the electron, ¢@ the scalar and A=(A,, Az, As) the 
vector potential, all quantities throughout this paper being expressed in 
Gaussian units. If A is so chosen that ¢is merely the electrostatic potential 
due to the nucleus, e@ can be included in Ho, the Hamiltonian function for 
the unperturbed atom. Also, if divA=0, p; and A; commute, and (1) reduces 
to 


H=Hot+e/uc-LAipite’/Quc?-2A?, i=1,2,3. (2) 


If we assume the incident radiation to be a plane polarized light wave 
of frequency v propagated along the Y axis, with the electric vector along 
the Z axis, the components of the vector potential may be taken to be A, 
=A,=0, A.=—cF sin 2mv (t—y/c)/2mv, where c is the velocity of light, 
and F the amplitude of the electric vector. We now neglect the size of the 
atom in comparison with the wave-length of the incident light,*® neglect the 
last summation’ in (2), and obtain the wave equation 


[Ho+(h/2mi)d/dt|V =(eF/2nvy) sin 2xvt- pV, (3) 

the right hand side of which we regard as the perturbation term. Let a 
solution of the unperturbed problem be Wo(m) exp [27E(n)t/hi|, so that 

[Ho— E(n) ]¥o(n) =0, (4) 


where stands for the three quantum numbers. We assume that the solution 
of (3) is of the form 


V(n) = { Wo(m)+W¥1(n)} exp [—2riE(n)t/h]. (5) 

Substituting (5) into (3), using (4), and neglecting the term in F¥,(”), we 
obtain the equation to be satisfied by ¥;(m): 

[Ho— E(n)+(h/2mi)d/dt|V1(n) = (eF/2evp) sin 2evt- p(n). (6) 


It is seen from (6) that to be periodic, with the period 1/v, ¥,(m) must be 
of the form 


WV, (m) = uje?*** — uge27*"* —" (7) 
where the w’s are functions of coordinates only, satisfying the equation 


heF dV (n) 


bs 2 


[Ho— E(n) + bv ]u=— 
wv Oz 





(8) 


6 Definite limitations on the wave-length of the incident light will be imposed later in 
connection with Eq. (15). 

7 Neglecting this summation is contrary to Dirac’s treatment of this problem, for he con- 
siders that this term is necessary to obtain the correct dispersion formula. The fact that the 
present method also leads to the correct dispersion formula can be traced to the difference in 
the assumed physical conditions in the two methods. Here we look for that solution of the 
wave equation in which the terms due to the perturbation contain factors periodic in time with 
the period of the incident radiation. The solution so obtained is interpreted by a generaliza- 
tion of the method of Heisenberg, Zeits. f. Physik 33, 879 (1925), and E. Schrédinger, Ann. d. 
Physik 79, 734 (1926). 
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where ~.¥% has been eliminated through its equivalent (h/271)0V)/dz, and 
where, as in what follows, the upper sign refers to the subscript 1 and the 
lower to the subscript 2. Thus to an eigenfunction Vo(m) exp [27E(n)t/hi] of 
the unperturbed atom there corresponds an eigenfunction ¥(m) of the per- 
turbed atom of the form: 


V(n) = Vo(m) + uye2*”* — ge 27 } en 2FiE (nish. (9) 


The form of the matrix elements corresponding to the dipole moment 
er =P, on neglecting terms involving two of the w’s as factors, is now 


P(nn')se [ Hn)ev(n')ar 


. (10) 
= [pe2tiv (mnt Ji e2rilety(nn’)]t 4 Jog—2eily—v(nn’)}¢ 
where v(nn’) = |E(n) —E(n')|/h, and where the IJ’s are functions of the 
quantum numbers only. 

We now interpret 2|Jo|, 2|1,|, and 2|Z,| just as we would the coeffi- 
cients of the Fourier expansion of the dipole moment in the classical theory, 
except that this dipole moment is here associated with a transition from the 
state n to the state m’. Thus 


| 1, | 2. 83 [y+v(nn’) |4/c3 (11) 


is the energy radiated with the frequency »+v(nn’) per unit time, per unit 
solid angle, in the direction perpendicular to r. For the case n=n’, I, =], 
P(nn') =2I, cos 2xvt, and 2J, is the amplitude of the equivalent oscillating 
dipole in the state m. In the same way J) and J; give the energy radiated with 
frequency v(nn’) and v—v(nn’) respectively. 

In order to obtain Dirac’s result, it is necessary to expand the perturba- 
tion terms (5), as well as the right hand side of (6), into a series of the form 
>. ,c(s)Wo(s). As a consequence of the use of Wo’s for this expansion the inten- 
sity coefficients Io, J;, and J, will contain a series, each term of which con- 
nects the initial and the final state with some other state. This tempts one 
to interpret the process of transition from a state m to another state n’ as con- 
sisting not only of the direct jump n—n’, but also of the double jumps 
n—s—n’', where s stands for any state which combines with both n and n’. 

2. In using Eq. (9) it is inconvenient to expand the w’s in terms of the W's 
because these form a complete orthogonal set only if the continuous range of 
the functions, corresponding to the hyperbolic orbits, or to the imaginary 
values of the total quantum numbers, is included. Instead, we use a set of 
functions given by 


1/2 —|xk ! 1/2 
T(o,A,x;a@) = {|| ew FO | | | | P,'*\(cos 0) ar2ar),” 








2(a+| «| )! 
where (2a)1-(@—r—1) 91/2 
ar) = = dll < —argr i a 
oe ore (o+2)! | ee nfo = 





a 
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in which @ is any constant, ¢, A, and « are integers satisfying the conditions 
o@>d= |x|, and L234" (2ar) is the Laguerre associated function. The 
functions T satisfy the differential equation 


V°T —(a?—2ac/r)T=0, (14) 


and it can be shown® that for any constant real value of a the functions given 
by all allowed values of o, A, and « form an orthogonal and normalized set 
which is complete with respect to a function decreasing rapidly with in- 
creasing 7, such as r0Vo/dz. We set 


a=a(n)=+ | —8x%u[E(n) = hv |/h?} 2, (15) 


and restrict ourselves to the values of v which yield a real a. This restriction 
is physically equivalent to our earlier assumption that the size of the atom 
is negligible in comparison with the wave-length of the incident light, though 
the equivalence should be regarded as accidental. We also set 


n=4nrye°Z/h?, (16) 
and (8), upon the substitution of Ho and a slight simplification, becomes 
V2u(n, 1, m) —(a?—2n/r)u(n, 1, m) =(eF/ hv) dV o(n, 1, m)/dz, (17) 
where all the quantum numbers have been written out. If we now put 
rOVo(n, 1, m)/dz==,,r, a(n, l,m; o, dr, x)T(o,d,K; a), (18) 
u(n, 1, m) =~,,x, . b(n, 1, m; o,d, x) T(o, d, K; a), (19) 
we find, with the help of (14), that (17) is satisfied if we set 
b(n, l,m; o,d, x) =eF a(n, l,m; o,d, x)/2hv(n—ac). (20) 


Thus the use of the W's for the expansion of the u’s is not necessary. The 
essential thing is to obtain a solution of (3) in the form (9). Since the T’s 
are not the eigenfunctions, the states s in the process n—s—n’ will not be 
the eigenstates, and the interpretation of the process as a double jump 
loses its significance, being a spurious relation dependent on the use of the 
W's. 
The method here presented, with a slightly different notation, was applied 
‘to the treatment of dispersion in a previous paper.'® As an illustration of its 
application to the Smekal-Raman effect we shall consider here only transi- 
tions between the states of principal quantum numbers 1 and 2, when ob- 
served by means of light scattered along the X axis. 
3. We shall be concerned only with the following particular cases of 
(10): 
P(n,l,m; n’, 1’, m’) =Io(n,1,m; n’, Ul’, m’; z)e2riv(nnt 
+I1,(n, l, mM; n’, Fe m’ ; z)e2tilytr(nn’)]¢ 
+I2(n,1,m; n’,l’,m’; z)e2rib—vnn 1 
8 E. Schrédinger, Ann. d. Physik 80, 486 (1926). 
® Courant and Hilbert, Meth. d. Math. Physik I, Ch. 5, §5, p. 239 (1924). 


10 B. Podolsky, Proc. Nat. Acad. Sci., 14, 253, (1928) and F. Reiche, Zeits. f. Physik 53, 
168 (1929). 


(21) 
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and similar expressions with z replaced by y. We let n=1, m=/=0, and 
n'=2; on performing the calculations it then turns out that all the I’s 
corresponding to the permissible values of m’ and l’ vanish except the follow- 
ing: 


To(1, 0,0; 2, 1,0; 2) =e f ¥.(1,0, 0)z¥o(2, 1, 0)dr, 
(22) 
To(1,0, 0; 2,1, £15 9) =¢ ff Fol, 0,0) ya(2, 1, +1)ér, 


7,(1, 0,0; 2, 0,0; ame f ¥u(t,0, 0)zu,(2, 0, 0) dr 


~~ f Zi2(1, 0,0)2¥0(2,0,0)dr, (23) 


12(1, 0,0; 2, 0, 0; ame f m(1,0, 0)z¥o(2, 0, 0)dr 


—e f Wo(1, 0, 0)zu2(2,0,0)dr. (24) 


We may write Yo(n, 1, m), properly normalized, in the form 
Vo(n, l, m) - (n'/?/n) T(n, l, mM; n/n) ° (25) 


Since we need u(1, 0, 0), and u(2, 0, 0), we first obtain, by the use of (25): 


d¥o(1,0,0) 1 1/2 3 1/2 , n® 1/2 “a 

ae AG) HG) to HaG) erp 
and 

d¥o(2,0,0) 1\ 2) ¢/3\17 nbe-""\ /2(4r — nr?) 
ee AG) HG) Peemah(-) Sf om 


Comparing Eqs. (26), (27), (18), and (12), we obtain: 





2re-""(n8/3)/2= DS a(1, 0,0; ¢,1,0)M.*(2a(1)r), (28) 
o=2 
and 
(4r—nr?)(n5e-""/96) 2 = So a(2,0,0; 1,0) M3(2a(2)r). (29) 
o=2 


Since the M’s satisfy the orthogonality relation 


f M(2ar) M,(2ar)rdr = dee: , 
0 
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we find that 
a(1,0,0; £,1,0)=26nt/3)"" f r2e-”M,'(2a(1)r)dr, (31) 
0 
a(2,0,0; o,1,0)=(a'/96)""2 f (4r?—nr®)e—"!2M,1(2a(2)r)dr. (32) 
0 


The integrals in these equations are calculated with the help of the 
relation 


J(s,B; 0,4; a)= f r*e-8*M(2ar)dr 


0 





(0+)! ey d 
(s—d—1)! (2a)* \dx 
where x = —2a/(a+8). 

We introduce the notation B(n) = —hv/E(n), qi(n) =[1+6(n) |", and 
go(n) = [1 —B(n)]2; then, since E(n) = —27%wetZ?/h?n?, Eq. (15) becomes 
a(n) =nqi(n)/n and a(n) =nge(n)/n. Substituting the relations (31), (32), 
and (33) into Eqs. (19) and (20), and theninto Eqs. (23) and (24), and elimi- 
nating the a’s by means of the q’s, we finally obtain 

Mer  S (o+1)! qi(2)—— { [1—gu(2) ] [2—g(2) ] }o- 
3 2hvn? fae (0—2)!(2—ogi(2) { [1+¢1(2) ][2+91(2)] } +8 
q2(1) [1 —2ga(1) ]-*[1—go(1) Jo? 
1—ogo(1) [1+2g2(1) ]o*4[1+-¢2(1) ]o+? 


s—r—1 
) {xtt¥(14+x)e>-1} : (33) 





[4 es ogi(2) | 


1 





[1—ogi(2)+91°(2)] — [4-9eq2(1) (34) 


+4(2+02)q2°(1) — 40q25(1) \\ 


The expression for J, is similar, and can be obtained directly from the 
expression for J, by substituting g,(1) for g2(1), substituting g2(2) for g:(2), 
and changing the sign in front of the whole expression from plus to minus. 

‘These formulas involve only summations and, in conjunction with (11), 
enable one to calculate the intensity of the Smekal-Raman lines with a shift 
corresponding to the first absorption line, provided that the frequency of the 
incident light is not too near a resonance frequency of the atom (note the 
resonance denominators), and that the size of the atom is negligible in com- 
parison with the wave-length of the incident light. 

We wish to acknowledge our indebtedness to Mr. Harvey Hall and to 
Mr. G. G. Harvey for much assistance in checking this work. 
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ON LIFE AND CONCENTRATION OF METASTABLE 
ATOMS AND THE QUENCHING OF MER- 
CURY RESONANCE RADIATION 


By E. GAVIOLA 
DEPARTMENT OF TERRESTRIAL MAGNETISM, CARNEGIE INSTITUTION OF WASHINGTON 


(Received October 2, 1929). 


ABSTRACT 

A complete and general study of the dependence of life and concentration of met- 
astable mercury atoms, optically excited, on experimental conditions, on foreign gas- 
pressures, and on intensity of illumination is given. This study explains why the life 
can not be larger than 10~* second under laboratory conditions and why foreign gases 
act so differently in regard to the accumulation of excited atoms. The theoretical 
results are in good agreement with all experimental data available. An improved and 
more general formula for the quenching of mercury resonance radiation is obtained. 
The efficiencies of collisions of the second kind of metastable mercury atoms with Ne, 
A, and He molecules are calculated numerically. It is shown that the concentration 
of metastable atoms increases only with the square root of the intensity of the exciting 
light and not with the exciting light itself as generally assumed. 


I. LirFE oF METASTABLE ATOMS 


N A previous publication! I have mentioned that a provisional measure- 
ment of the transition probability of the “forbidden” line 2656 (2*P»—1'S,) 
indicates that the life of a metastable mercury atoms is, when undisturbed, 
of the order of magnitude of 10 seconds. On the other hand all the measure- 
ments of that life made under laboratory conditions give values for it that 
range from 10-* to 10-* second. The gap between the two values is consider- 
able. It is therefore of interest to analyze the factors that in the laboratory 
shorten the life so many times and to 























see how these factors work depending 2°F | t N, 
on experimental conditions. The me- ze, Pa 

tastable atoms that we shall consider ~ | Bi 

are the ones on the level 2*P». They are = s | | ~~ 
the only ones of importance for photo- A Ze, Yaz 
sensitized fluorescence and chemical 

reactions. Metastable atoms are pro- 71 | ata 6 

duced mainly and practically ex- | 

clusively in all the cases we are 1'S 

considering by collisions of the second Fig. 1. Diagram of transitions regard- 
kind of resonance atoms (2°P;) with ing metastable atoms. 


other molecules (Fig. 1). They are 

destroyed in many different ways which we shall consider one by one. 
Let us begin by considering a vessel with a quartz window containing 
saturated mercury vapor at room temperature and illuminated with the 


1 E. Gaviola, Phil. Mag. 6, 1181 (1928). 
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light of a water-cooled, magnetically-deflected, mercury arc. The intensity 
of illumination should not be too strong. The case of strong illumination 
will be considered further on. Under these conditions the few metastab!e 
atoms that will be produced will almost with certainty reach the wall 
owing to the long mean free path at 10-* mm pressure and they will be 
destroyed as a result of it. All experimental evidence favors the assumption 
that collisions with solid walls destroy all excited atoms including metastable 
ones. Now owing to the fact that the resonance line is reduced to one-half 
5 mm vapor, the mean distance of a metastable atom from the entrance wall 
when it is produced is 5 mm and since the mean thermal velocity is 1.7 K 10* 
cm/sec, the mean time that atoms will take for hitting the entrance wall 
will be of the order of a small multiple of 10-5 second. The presence of the 
other walls of the vessel will decrease that time. The mean life of a metastable 
atom “in vacuum” can therefore not be longer than say 1/2 X10~‘ second. 
This time could be increased only by lowering the temperature and at the 
same time increasing the dimensions of the vessel. The latter alone would not 
help owing to the fact that at room temperature most metastable atoms will, 
be produced only a few millimeters away from the entrance wall and they will 
have a good chance of hitting this wall independent of the distance of the 
others. Now let us assume that we introduce slowly into the vessel a foreign 
gas which does not or nearly does not destroy metastable atoms by collisions 
of the second kind. As such we may use H2O, Ne, A, or He, which have been 
freed from hydrogen. The life, 7o, of the metastable atoms will now increase 
in proportion to the foreign gas pressure owing to the diminishing of the 
diffusion rate towards the walls. At the same time the number of collisions 
of a metastable atom with foreign gas molecules will increase in proportion 
to the pressure. Now Foote? has called attention to the fact that at room 
temperature out of 6000 molecules one has a speed sufficient to perform a 
collision of the first kind bringing the metastable atom to the resonance level. 
A metastable atom will survive then on the average only 6000/£, collisions, 
if E, is the efficiency factor for collisions of the first kind. The life, 70, will 
then reach a maximum value when on the average a metastable atom suffers 
6000/F; collisions before hitting the walls. The maximum will be reached 
at a pressure of a few millimeters depending on the factor Z;. If we increase 
the foreign gas pressure, 7» will be still mainly determined by the time of 
6000/£E, collisions, but since this time decreases with the pressure, 70 will 
decrease in proportion to the inverse pressure. Let us put all this in better 
form: The probability of hitting the wall is 1/(c+2Zd), if c is the probability 
at pressure zero, d a constant inversely proportional to the diffusion coeff- 
cient and Z the number of gas kinetic collisions, which is proportional to the 
pressure. If collisions with the wall were the only death reason, the life, which 
we will call 7;, would be the inverse of the probability given above 


™m=c+dZ. (1) 


This is a straight line which cuts the ordinate axis at the point c and has a 
slope d that depends on the diffusion constant of the particular gas and especi- 


2 P. D. Foote, Phys. Rev. 30, 288 (1927). 
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ally on the position of the volume element that we are considering in regard 
to the walls; for a volume element near a wall,d will have asmaller value than 
for an element farther away from it. 

On the other hand the probability of suffering a collision of the first kind 
is ZE,/a. Here a is the average number of collisions required for one of them 
to have sufficient energy to perform a collision of the first kind. This number 
is given by Boltzmann’s formula a =e~“/*? where w is the difference of energy 
between the 2°P; and 2°P» states. The life, considering only collisions of the 
first kind, is then 


T2=a/ZE. (2) 


The resultant life vo is, if we take both processes into consideration, the in- 
verse of the sum of both probabilities 


1/ro= 1/rit+1/72. 


To simplify the formulas we will neglect c in (1) because it is small in compar- 
ison with dZ for all pressures above a tenth of a millimeter. We write then 


™1=dZ (3) 
and obtain with it 
to=1/(ZE,/at+1/dZ). (4) 


Figure 2 illustrates 71, 72 and 7 as functions of the pressure Z assuming arbit- 
rary values for d and E;. We see that 79 is mainly limited by wall collisions 
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Fig. 2. Dependence of life and concentration on foreign gas pressure. 


at pressures below one millimeter, while it is determined mostly by collisions 
of the first kind beyond 4 millimeters pressure. 7) can be increased only by 
lowering the temperature (increase of a) or by choosing a gas with small 
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efficiency E;. Differentiating 7) in regard to Z we find that 7) has a maximum 
for ZE,/a=1/dZ, that is, when the probabilities of hitting the wall and 
suffering a collision of the first kind are equal. Now, what are the numerical 
values of E, and d for the gases that come into consideration? Under the 
assumptions made so far, metastable atoms are destroyed at higher gas pres- 
sures only by collisions of the first kind. Under stationary conditions the 
number of transitions producing metastable atoms N,Z’E,’ (Fig. 1) must be 
equal to the number of transitions destroying those which are NoZE;/a. 
We have then 
N.ZE\/a = N,Z'Ey’ 


and since for higher pressures No =a JN, as we shall see presently, 
ZE,=2Z'E,’. (5) 


If we set arbitrarily Z = Z’, which amounts to saying that the collision section 
of the metastable atoms shall be set equal to that of resonance atoms in 
every case then we see that the efficiency of collisions of the first kind, of 
metastable atoms, is equal to the efficiency of collisions of the second kind, 
of resonance atoms. These last efficiencies have been calculated in a pre- 
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Fig. 3. Life, ro, as a function of the foreign gas pressure under ideal conditions. 


ceding paper*® to which we shall refer in the following, and we shall use the 
values given in Table II there. 

The value of d depends primarily, as we said before, on the position of the 
volume element considered in regard to the walls and on the diffusion con- 


3 E. Gaviola, Phys. Rev. 33, 309 (1929). 
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stant for metastable atoms in each particular gas. This last is not known. 
The dependence of d on position can be calculated exactly, but we shall not 
do it here to avoid lengthening the paper unnecessarily. We shall assume 
arbitrarily d=2X10-*P/Z sec/mm (P is the pressure in millimeters) in the 
following for all gases. This value is in accordance with previous calculations 
of similar cases‘ and represents more or less an average for the experimental 
conditions primarily considered here. Using these numerical values we ob- 
tain the curves of Fig. 3 fort. The ordinates give the life in seconds and the 
abscissae the pressure in millimeters. 

The influence of collisions of the second kind and of impurities. So far 
we have assumed that collisions of the second kind of metastable atoms 
with foreign gas molecules have an efficiency zero. Now if an impurity in 
the form of He: or Oz is present in the foreign gas, which will generally be the 
case if no special precautions are taken, this efficiency can not be zero. 
Moreover experimental evidence seems to show that some gases like CO 
have an efficiency of their own in destroying metastable atoms independent 
of impurities. Let us consider now how 7pis affected bycollisions of the second 
kind of metastable atoms. 7 will be limited not only by the probabilities 
of hitting the wall and collisions of the first kind but also by the probability 
of a collision of the second kind ZE, if Eo is the efficiency factor for this last 
type of collisions. Thus 


. i 1 Ei _ 
— =—+—+ZE,>=— + “(—+£0) 
dZ a 


TO dZ Q 
and if we make 
E\/at+E,y=E (6) 


then E is the total efficiency of a collision in destroying a metastable atom. 
T) now takes the form 


to=1/(ZE+1/dZ). (7) 


This expression is the same formula (4) of the previous case with the only 
difference that instead of E,/a we have now E, which is in general slightly 
larger. The curves of Figure 2 then will still represent the form in which 
To Varies with pressure, the ordinate will have to be reduced in general, and 
for higher pressures the reduction factor will be E;/aE. To give a numerical 
example, if we use N2 with so much impurity in it that E=2E,/a then the 
curve that will represent N, in this case is the one that corresponded to water 
in Figure 3. (This is because FE; water =2£, nitrogen.) 

The influence of formation of molecules, absorption of light and impurities 
of the vessel. It is known that in the case of a collision of a metastable mercury 
atom with a normal one there is a certain probability for the formation of 
a Hg: molecule. This probability will not depend on the pressure of the for- 
eign gas. The same is true for the probabilities of absorbing lines of the arc 
like 4046 and of colliding with some impurity (Hz) which develops in the ves- 


4 Meissner and Graffunder, Ann. d. Physik 84, 1041 (1927). 
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sel as a result of the illumination. To these we can add the probability of 
emission of 2656. We will denote by } the sum of the probabilities of all these 
processes which do not depend on the pressure. Then 


to=1/(ZE+1/dZ+5) (8) 


is the most general expression for the life of metastable atoms. 

The term 0 will be small in comparison with 1/dZ for small pressures and 
it will be negligible in comparison with ZE for large Z; only for intermediate 
pressures will it affect the curve 7» in which case the effect will be to flatten 
down the maxima of the curves of Figure 2. Altogether then we see that the 
form of the 7» curves is mainly determined by diffusion on one side and by 
collisions of the first and second kind on the other, and that the value of 
T) can not under laboratory conditions be increased practically beyond say 
10-* second at room temperature. 


II. RELATIVE NUMBER OF METASTABLE ATOMS 


The number of metastable atoms will be under all conditions given simply 
by the number of such atoms produced per second multiplied by their life 
To. Metastable atoms are produced in all cases considered here mainly and 
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Fig. 4. Relative concentration of metastable atoms as a function of the foreign gas pres- 
sure under ideal conditions. 
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practically exclusively by collisions of the second kind of resonance atoms 
with foreign gas molecules. The number of atoms brought to the 2°Po 
level per second is therefore N,ZE, [Figure 1 and Eq. (5)] and the number of 
metastable atoms 

No=NiZEj10. (9) 


The ratio No/N; has been plotted in Figure 2 as a function of the pressure. 
We see that this ratio will increase at first with the square of the pressure 
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owing to the fact that both ZE, and 79 increase with pressure and that it will 
rapidly approach the asymptotic value a. As we see, the form of the No/ MN; 
curve is very different from the one of the 7») curve and the maximum in the 
life will in general not correspond to a maximum in the concentration. 

Let us now calculate No under the ideal assumptions Ey=0 and b=0. 7 
is then given by formula (4) and 


No 1 1 “— 10 

Ni (- +i) ( ) 
In Figure 4 I have plotted this ratio for the gases we are considering as a 
function of the pressure using for a, d, and E, the same values as before. We 
see that all the curves tend more or less rapidly towards the asymptotic 
value @ given by Boltzmann’s formula. A comparison of Figures 3 and 4 
shows the curious result that for a given pressure the gas that allows the 
metastable atoms to live longest is the one which produces the smallest 
concentration of them. This shows again that life and amount by no means 
behave similarly.* 





CasE Ey) ¥ 0 


If we assume that E)~0 the life 79 is given by formula (7) and the amount 


No will be 
No E . + 
~=ZEn=(=+ =) 
N, E, dE,Z* 





and if we put 
E,/E=a! (11) 


No 1 ;T\ (12) 
N; (Sta , 
This formula is identical with formula (10) with the only difference that in- 
stead of a we have now a’ which is somewhat smaller than a. Formula (12) 
will be represented then by curves similar to the ones in Figure 4 that will 
converge asymptotically towards a constant value a’. This yields the inter- 
esting result that the Boltzmann distribution between 2*P» and 2°*P; will not 
hold now any more even at high pressures in spite of the rapid interchange 
between the two levels. 
The fact that No/N,; becomes constant for high pressure does not mean 
that Np itself will be constant. It was so in the ideal case considered before 


we have 





* Formula (10) and the curves of Figure 4 are not exact for any particular volume element 
v. They are exact for the average taken over the whole volume. This is due to the fact that 
the metastable atoms generated in a volume element v, and which have an average life ro, 
do not stay in v owing to diffusion; they distribute more or less all over the vessel, depending 
on the pressure. At higher pressures when the influence of diffusion disappears, formula (10) 
and curves of Figure 3 will be exact also for a particular volume element. An exact calculation 
of the dependence of No on pressure and position can be made as will be shown later. 
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when it was assumed that E=0. This however, leads to Ni= constant and 
N, will no longer be constant in the case Ey) =0. 


Case b + 0 


Let us consider the influence of formation of molecules, absorption of 
light, and of impurities of the vessel on the amount No. The life 70 is given 
now by formula (8). Therefore 


~=ZEsro=(—+ +) (13) 
nN, ii a’ dE,Z? ZE, 





The term 6/Z£E, will be small compared with 1/dE,Z? for small pressures and 
it will be negligible in comparison with 1/a’ for high pressures. The curves 
of Figure 4 will then represent (13) adequately for small and for large 
pressures (replacing a by a’) but for the intermediate range the presence of 
b/ZE, will flatten down the curves. 

To obtain No itself as a function of the pressure we have to multiply 
(13) by M;. The values of Ni=f(Z) can be obtained from the quenching 
curves of Stuart® or they can be calculated directly. We shall do the latter 
because with the help of Stuart’s curves we shall be able to obtain the numer- 
ical values of Ey and of } for some of the gases. 


III. THE QUENCHING OF RESONANCE RADIATION 


The quenching of mercury resonance radiation originally calculated by 
Foote in a very interesting paper” has been recently treated by Klumb and 
Pringsheim® in a more complete and very elegant way. The latter authors 
restricted themselves to the calculation of No and did not take into consider- 
ation the reabsorption factor g(1—f). To avoid repetition, I have used here 
Klumb and Pringsheim’s results in the calculation of Ni paying attention to 
the factor g(1—f), and without neglecting the term 7/7.° (in Klumb and 
Pringsheim’s symbols). I give here therefore the expression for N; and refer 
to the publication’ for its deduction. It is 


(1-4) M, -{1 (ZEo+b)r/(1—q/2) Ve 
2/N,° 1/a’+1/dE,Z?+6/ZE, 








=J (14) 


where N,°/(1—q/2) is the number of resonance atoms for Z=0, 7 their 
mean life at Z=0, and g the reabsorption factor. (1—g¢/2)Ni/M.°=J is 
proportional to the amount of resonance radiation observed and it rep- 
resents the value directly measured by Stuart. For pressures above a few 
millimeters diffusion can in general be neglected and therefore we can put 
1/dE,Z*?=0. Then (14) goes over into 


{ 1 Cte —q/2) ye 
1/o'+b/ZE; 


J= 





(15) 


5 H. Stuart, Zeits. f. Physik 32, 262 (1925). 
* H. Klumb and Peter Pringsheim, Zeits. f. Physik 52, 610 (1928). 
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Formula (15) is identical with formula (10) of Foote? if we translate his 
constants in the following way: p—Z, cb, b—-E,, d—E,/a, e—E» and fr- 
(1—q/2)/r. We can therefore rectify (15) in the same way in which Foote 
did with his formula. We can thus write it as follows: 

ar/(1—q/2) «a 1 


1/J-1 ZE, ZEo+b 








(16) 


If we plot therefore the inverse of the quantity at the left using Stuart’s 
values for J, we will obtain for pressures above a few millimeters a straight 
line, the slope of which will give us Eo, and its value for Z =0 will be equal to 
b. Table I gives the values used and the results obtained in that way. We 
have set g=1/4 as in the paper referred to in footnote 3. 

















TABLE I 
Gas Z ZE\ Ey b ZEo Eo 
Co 28. X10%p 28. X 10° p 1 -— = — 
H,O 6.4 is y ¥ 0.4 = + _- 
Ne 6 " 1 ws 0.2 1300 87 is. <10-* 
A 5 - 0.26 ie 0.05 570 10 2 ‘ 
He 14 ° 0.071 ” 0.005 330 3 0.2 . 


























The values of 6 and Ey for CO and H,0 cannot be obtained in this way 
due to the neglect of diffusion. We have seen that our calculation is correct 
only for pressures above a few millimeters; the quenching curves of CO and 
H.O are only known for a few millimeters pressure and therefore an estima- 
tion of 6 and of Eo cannot be made in this way for these gases. The values 
of b and E> obtained for Ne, A, and He are probably correct and give an 
idea of the experimental conditions of Stuart. We cannot ascribe to these 
values any specific meaning because 0 is the sum of the probabilities of four 
independent processes: emission of 2656, formation of molecules, absorption 
of light, and collisions with vessel impurities; and Ep embodies collisions of 
the second kind with foreign gas molecules and with impurities that may be 
present. If we arbitrarily interpret Eo as solely due to impurities (H2) we 
would have to assume that Stuart had one molecule of hydrogen in 10,600 
of Ng, in 66,500 of A, and in 1,860,000 of He. These values differ greatly from 
the ones obtained by Foote assuming E,=1 for all processes (see his Table 
II), and are more likely to correspond to experimental conditions. 


IV. ABSOLUTE NUMBER OF METASTABLE ATOMS 


Now that we have obtained the expression for N; and the numerical 
values of Ey and } we can return to our metastable atoms. Multiplying (14) 
and (13) we obtain 


No=N1°/(1—9/2)/ [1/a' +(ZEo+b)1/(1—9/2) +1/dE,Z?+b/ZE,| (17) 


as the final expression for the amount of metastable atoms as a function of the 
pressure. We know now all the constants entering in (17) and we can there- 
fore plot No as a function of the pressure Z. This has been done in Figures 
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5 and 6 for the gases we are considering. We would obtain the same curves 
for Ne, A, He, using formula (13) and for N; the values taken from Stuart’s 
curves. The ordinates of Figures 5 and 6 indicate the number of metastable 
atoms in a volume unit assuming that N,°=1; that is the unit of volume ab- 
sorbs 10’ light quanta of 2537 per second. As long as the illumination is not 
too strong the values of No for each pressure will be directly proportional 
to the absorbed light intensity. We shall consider presently the case of 
strong illumination. As said before we could not obtain the values of 6 and 
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Fig. 5. Number of metastable atoms as a function of the foreign gas pressure 
if 10’ quanta are absorbed per second. 
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Fig. 6. Number of metastable atoms as a function of the foreign gas pressure 
for N,°=1 (107 quanta absorbed per second). 


E, for CO and for H,O owing to the neglect of diffusion. The curves for these 
two gases in Figures 5 and 6 have been obtained using formula (17) and as- 
suming arbitrarily }=1300, Ey>=14X10-* for H.O and 6=1300, Ey=177 
X10-* for CO. These values have been chosen so that for No they give curves 
that are in agreement with the previous experimental results of the author.' 
They probably do not correspond to the experimental conditions of Stuart. 

Let us see how the values of Ey and b influence the shape of our curves No. 
A large Eo (impurities proportional to the pressure) will displace the maxi- 
mum of No towards smaller pressure without flattening it much and at the 
same time it will make Np» decrease rapidly for higher pressures; a large } 
(impurities independent of pressure) will displace the maximum towards 
higher pressures and it will flatten it considerably. This last has been experi- 
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mentally shown by Klumb and Pringsheim admitting a constant small 
amount of He together with the nitrogen® (Figure 4, p. 618 in Klumb and 
Pringsheim’s paper). The curves of Figures 5 and 6 are in general qualitative 
and sometimes even quantitative agreement with the experimental results 
of Donat,’ Loria,* Klumb and Pringsheim,® Pool,’ and Gaviola.! So for 
instance, the curves for the intensity of the photo-sensitized fluorescence of 
thallium as a function of the pressure of argon obtained by Donat are in 
very good agreement with the A curve of Figure 6; he found maxima for all 
the lines at 40-mm pressure and we see that this corresponds exactly to the 
maximum in the amount of metastable atoms. Since according to Loria at 
low pressure, Neis more efficient than Aand according to Donat the reverse 
is true at high pressures, we see that the apparent contradiction between 
their results is removed and conciliated by Figure 5. The fact that Klumb 
and Pringsheim obtain a maximum in the nitrogen curve already at 1 mm 
pressure, would indicate that they measured the absorption of 4047 at some 
distance from any wall so that their diffusion curve (7; in Figure 1) was very 
steep, which corresponds to a large d. The result obtained by Klumb and 
Pringsheim with water is almost certainly due to the pressure of free hydro- 
gen mixed with the water vapor. According to my experience water vapor 
contains always some free hydrogen even before it is illuminated with the 
arc. Finally the curve obtained by Pool for No(a) as a function of the nitrogen 
pressure is in very good agreement with our Ne» curve (Figure 6). The curve 
of Pool for 7o(b) as a function of the nitrogen pressure is also in good general 
agreement with our curve of Figure 3. It ought to be mentioned though that 
according to our theory the 7ocurve ought to increase linearly starting from 
zero, which the 7o(b) curve of Pool does not do. I am not able to explain this 
divergence. 

For a given illumination the maximal concentration of metastable atoms 
that can be obtained is then primarily limited by the ratio a=e-“!*? and by 
the amount of impurities contained in the admixed gas. a can be increased 
by lowering the temperature 7, but this would diminish rapidly the density 
of the mercury vapor in the vessel, so that no real gain would result. We can 
increase No only by increasing the intensity of the primary light, but this 
has also a limit as we shall see in the following paragraph. 


V. DEPENDENCE ON ILLUMINATION 


So far we have disregarded collisions between two excited atoms. As long 
as the illumination is not strong the number of excited atoms will be small 
in comparison with the number of normal atoms and collisions between two 
metastable atoms will occur rarely. But if we increase constantly the in- 
tensity of the exciting light, the number of metastable atoms will increase 
at first linearly with it and the probability of a collision of two metastable 
atoms will grow with the square of the illumination. Now Beutler and 

7K. Donat, Zeits. f. Physik 29, 345 (1924). 


8 S. Loria, Phys. Rev. 26, 573 (1925). 
® M. L.Pool, Phys. Rev. 33, 22 (1929). 
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Josephy’’ have shown that by a collision of two metastable atoms there 
is a great probability for one of them to take over the whole energy, the other 
becoming normal. Both metastable atoms cease then to be such as a result 
of the collision. This process will put an additional limit to the life and con- 
centration of metastable atoms. 

The life was given in the general case by formula (8). It will be diminished 
now by the additional probability of hitting another metastable atom, which 
will be proportional to their number No. We will call it cNo. Then 


to= (ZE+1/dZ+b+c¢No)-!=(a+cNo)— (18) 
if we put the part independent of the illumination 
(ZE+1/dZ+5b)-'=1/a. (19) 
On the other hand, the number N is given always by formula (9) 
No=ZE,\Nit0 


in which J, is proportional to the intensity of illumination 7. We can write 


then 
No=hl to (20) 


where h is a factor of proportionality. Combining (18) and (20) we obtain 


the quadratic expression 
cN,?+aNo—Al=0 (21) 
and its solution is 
2cNo=(a?+4chI)'/*—a. (22) 


As we see Np will increase only with the square root of the light intensity. 
We obtain in the same way for the life 


to= [(a?+4chI)'/?—a]/2chl. (23) 


To will then decrease at the same time as No increases. Now Klumb and 
Pringsheim® have observed that the absorption of the line 4047 does not in- 
crease in proportion to the illumination J if N2 at pressures above 0.05 mm 
is present. The absorption reaches a certain saturation at a value of 0.55. 
They explain this by pointing out that since the cold vapor absorbs only the 
core of the line of the arc, once the whole core is absorbed no further absorp- 
tion can take place. Their explanation is surely correct. It may be but one 
factor in the phenomenon since formula (22) shows that we can expect 
No to grow only with the square root of the illumination. The curves of 
Figure 6 in Klumb and Pringsheim’s paper could very well be reproduced 
using formula (22) assuming a suitable value for c. A repetition of Klumb 
and Pringsheim’s experiment using two absorption vessels of equal section 
and different length would show how much of the saturation in the absorp- 
tion is due to the different widths of the arc and absorption lines and how 
much to collisions between two metastable atoms. 

The considerations and formula developed here for the metastable atoms 
of mercury apply, mutatis mutandis, to any metastable molecules excited 
optically. 

10 H. Beutler and B. Josephy, Phil. Mag. 5, 222 (1928), Naturwiss. 15, 540 (1927). 
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ENERGY CHANGES BY A VARIATION FROM 
THE CRYSTALLOGRAPHIC GROUP 


By H. M. Evyjen 
CALIFORNIA INSTITUTE OF TECHNOLOGY, PASADENA 


(Received October 11, 1929). 


ABSTRACT 

The purpose of the paper may be briefly stated as follows: To investigate the 
possibility of energy minima by a systematic departure from the perfectly symmetrical 
crystalline lattice. A geometrically simple variation is considered. The surface 
(100-plane) of a heteropolar crystal of the face-centered type is contracted uniformly 
in patches, and the interior planes are likewise contracted in patches by amounts de- 
creasing linearly with the distance from the surface. In this way a “mosaic” structure 
of a type previously suggested in the literature! is obtained. A numerical calculation of 
the change in potential energy is carried out for rock-salt with certain simplifying as- 
sumptions. The change in energy is obtained as a power series in terms of the maxi- 
mum displacement of the ions. Due to the slow convergence of this alternating series 
an approximate remainder term is introduced to take care of powers higher than the 
third. The series contains two parameters, namely the linear extension of the patches 
and the depth of the cracks. The function, in general, is found to have two minima 
with respect to the displacement of the ions, one at zero displacement, and a second at 
a finite displacement. The second minimum, depending on the parameters, may be 
negative with respect to the former. The function also has minima with respect to 
the two parameters, thus proving the static stability of a mosaic structure of 
definite dimensions. These considerations with only minor changes, are also appli- 
cable to a mosaic structure in the interior of the crystal. 


: ; ‘HE potential energy of two ions may be written in the form: 
E=-—(e?/r)+A/re (1) 


where the first term represents the energy due to the inverse square Coulomb 
forces, and the second term phenomenologically stands for the energy due 
to the repulsive forces between ions, which are necessary for static equili- 
brium. A and a@ in general depend on the kind of ions considered.? In the 
following, however, we shall confine our attention to crystals of the rock- 
salt type and shall assume that these constants are the same for each of the 
three combinations of ions which occur. This greatly simplifies the calcu- 
lations, and is permissible in first approximation, in as much as we are not 
primarily interested in quantitative results. In the case of rock-salt then, 
the best value of a will be that corresponding to a Na-ion against a Cl-ion, 
namely 9.2, With this constant fixed and the value of the lattice constant 
known, the value of A may be computed in the usual manner from the 
minimum energy condition. With Bragg’s® value for the lattice constant of 
rock-salt, namely 2.81 X10-* cm, we obtain A = 25.8 X10-®. 


1 F, Zwicky, Proc. Nat. Acad. Sci. 15, No. 3, March, 1929. 

2 Lennard-Jones and Dent, Proc. Roy. Soc. Al21, 247 (1928). 
’ Bragg, X-Rays and Crystal Structure, 306, (1924). 
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This calculation is based on the assumption of an infinite space lattice, 
perfectly symmetrical, with a lattice constant everywhere the same. Need- 
less to say, this assumption has been amply justified in the case of a great 
number of physical properties. In the case of certain other properties,‘ 
there has been a great discrepancy between theory and facts, a discrepancy 
which has recently led to a questioning of the above mentioned assumption.° 
It is then a matter of considerable interest to investigate the energetic 
aspects of a systematic departure from the space lattice of perfect symmetry. 
As to the type of variation to be considered there are several possibilities. 
These might perhaps be classified as follows: (1) Variations utilizing the loss 
in potential energy by contraction of planes, lines, and dipoles, from the 
lattice constant of the solid to that of the plane, line, or dipole respectively. 
(2) Variations utilizing the loss in energy by expansion of planes of homo- 
genous polarity. (3) Variations utilizing the former two principles in com- 
bination with a rearrangement of the ions leading to losses in energy com- 
parable with the heat of sublimation. The most favorable variation probably 
would be of the third type, and probably would be geometrically somewhat 
complicated. If, however, we can establish the existence of energy minima 
outside of the crystallographic group by a geometrically simple variation of 
the first type, our main purpose will be fulfilled. Very reliable quantitative 
results can hardly be expected anyway, except perhaps as to order of magni- 
tude, since the existence of energy minima outside of the crystallographic 
group would in some measure invalidate our previous computation of the 
ionic constants a and A, constants which we shall have to use in the estab- 
lishment of the minima. 

The following variation, then, will be considered: Starting with an ideal 
rock-salt crystal, a boundary 100-plane will be contracted in square patches, 
the contraction of two neighboring ions being everywhere the same, namely 
d/(p—1)=6, except at the edges where the displacement will be d. is 
the number of ions in a line element of the patch. During this first step of 
the variation the interior lattice planes are kept rigid. The resulting change 
in energy will be due to two sources: (1) Change due to the ions in the con- 
tracting plane, and (2) change due to the ions of the interior planes. As a 
second step in the variation the interior planes are similarly contracted by 
amounts decreasing linearly with the distance from the surface. 

The energy of an isolated plane, contracted as outlined above, may 
be written: 

1 x x P P (-— 1)+ me? 
. 2 pa a 22 tee kid)? + (ma — kod)? |"/? 





=—~.9 M=—-@ r=l s=1 


(2) 





A 
+ 
[(la — ky5)? + (ma — kod)? |= 
where a is the lattice constant, and k; and ke are integers depending on / and 


r, mand s respectively in a manner such that 


* See f. ex. F. Zwicky, Phys. Zeits. 24, 131 (1923). 
5 See reference 1. 
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where the symbol R(m/p) stands for remainder in m/p. 

The change in energy now may be obtained by expanding the expression 
2 in powers of d/a=6,, and summing the pure number series occurring as 
coefficients in this expansion. The expansion converges for d/a<1. The 
coefficients of the three first terms in this expansion have been evaluated.® 
The higher orders in the expansion may, to good approximation, be said 
to arise entirely from the ions forming the boundaries of the patch. This 
enables us to replace the higher orders by an approximate remainder,® so 
that for reasonably great values of p(p> 10) the following functional relation 
is obtained between the average change in energy per ion, AE and 6,: 


padE/e?=1.9417—1.47536.+ {1.14p/(p—1) —0.6235} 8.2 
— {0.04—7.66/(p—1)?}6.8—2/(1+6.)+0.0583/(1+8,.)® (4) 


In addition to this there will be a decrease in energy due to polarization 
caused by the asymmetries of the mosaic configuration. This change in 
energy will be given by: 

Pp P 


pae=-< Saar} 


t=1 q=1 





where the first term in a Taylor expansion gives: ® 





“ 2ek © 5 ki(m*—2I*) —3kelm 
AXte— 2 et es Rempe 
so that the second order change in energy will be: 
padE/e* = —0.366(p/p—1)?-6.? (5) 


a being the average polarizability of the ions.’ 
The higher orders here may be neglected, so that the change in energy due 
to the ions of the contracting plane may be written: 


padE/e? = 1.9417 —1.47536.+ {1.14p/(p—1) —0.366(p/p—1)?—0.6235} 5.2 
— {0.04—7.66/(p—1)?}53 —2/(1+6.)+0.0583/(1+6,)® (6) 


The energy of a displaced ion in the surface due to the interior of the 
crystal may be written: 


*H. M. Evjen, Thesis, California Institute of Technology, 1929. 
7 Born and Heisenberg, Zeits. f. Physik. 23, 388, (1924). 
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where /; and he, depend on p in such a way that 
1 P 1 P (p? = 1) 
_— hy? =— he? = } P 
p 2; a, ~ P 12 (8) 


Again an expansion by the binominal theorem gives the change in energy. 
Only even powers of d occur in this expansion. The second order change has 
been evaluated.’ If account is taken of the Lennard-Jones contraction per- 
pendicular to the surface,’ it is: 


padE/e? =0.0228p(p+1)/(p—1)8.2 (9) 


The total change in energy by the first step in the variation now may be 
written: 


padE/e?=1.9417 —1.47536.+ {1.14p/(p—1) —0.366(p/p—1)? 
—0.6235+0.0228p(p+1)/(p—1) }5.2 
— {0.04—7.66/(p—1)?}5.8—2/(1+8.) 
+0.0583/(1+6,)°. 

From which it is seen that this first step in the variation, which in itself 

forms a highly artificial variation out of the crystallographic group, already 


gives rise to energy minima. For certain values of p the change in energy 
may be represented as depicted in Fig. 1. The function also will be seen to 


(10) 


BE 








Fig. 1. Change in energy by a variation from the crystallographic group. 


possess a minimum with respect to p, showing that the patches, if formed 
in this way, would tend to assume a definite size. 

The simultaneous minimum in 6, and p, however, falls outside of the 
region of convergence of the expansion, so that the value of p for which this 
occurs can only be roughly estimated. It is of the order of 10. 
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The smallness of the patches formed in this way is due to the restraining 
influence of the fixed space lattice represented by Eq. (9), which is distin- 
guished by the fact that it contains p to a higher power than the contri- 
bution of the plane, (Eq. 6). When we carry out the second step of our 
variation this restraining influence will be seen to be practically annihilated. 
That is, when the boundary plane is already contracted, the forces of inter- 
action between the planes will tend to contract the neighboring interior 
planes. It is not deemed necessary, therefore, actually to compute the energy 
of interaction between the planes by the second step of the variation, as an 
upper limit to the total change in energy of interaction is given by the expres- 
sion (9) divided by the number of planes contracted; and this will be seen to 
be a rather small contribution, since the term now has lost its distinction 
of containing p to a higher power than the contribution of the planes. In 
the place of the actual value of this positive contribution we shall therefore 
use the upper limit. 

Aside from this small change in energy due to the interaction between 
the planes, the change will be made up of a sum of expressions of the type 6. 
If this expression be denoted by f(6.), the sum in question will be: 


S= ¥ f(e-q/n). (11) 


In summing the last two terms of f(6,-q/n), the summation may be 
replaced by an integration, so that, writing =n, the total change in energy 
by the complete variation may be written: 


padE/e? = 1.9417 —0.73765.+0. 5026 2—0.01165,3 
—(2In(1+6.))/5.+(0.0073/6.) [1—1/[1+6.]®]+o(m,x) (12) 
where 


5. 0.136 
o(n, x) = 0.02282x6 .24+— | - .0.736+ (0. 121 +=) 6.+(0.02+0.26/x 


nN x 


. (1.915/x2?—0.022)76. 
+0.014/ —| -| 0.022 |= (13) 





n nN nN 


2 6. 2 1+6, 
——In (1+~) —-—ln— 
‘ n Me 1+6e/n 


n 0.0073 1 0.0073 1 
S| save al 
nm—1 6, (1+6,/n)8 nde [1+6.]§ 


This functional relation again may be represented graphically by the 
type of curves shown in Fig. 1. The function possesses minima with respect 
to x and as well as 6,, showing that the mosaic structure, if formed in this 
way, would have definite dimensions. Again the absolute minimum occurs 
outside of the region of convergence, and the values of x and m for which 
maximum stability obtains can be only roughly estimated. x seems to be of 
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the order 1, and m of the order 25, which agrees roughly with previous phe- 
nomenological estimates.® 

If the previous calculation had been carried out for planes very far 
removed from the boundaries of the crystal the final result would have been 
essentially the same. The only difference in the calculation would have been 
the introduction of factor 2 in the constant of expression 9, which again 
would mainly have the effect of decreasing somewhat the value of the ratio 
x at the absolute minimum, an effect which we should expect. 

Beyond establishing the possibility of energy minima outside of the 
crystallographic group, the previous calculation is of interest in that it 
suggests the possibility of a coexistence of different modifications of the 
same crystallographic group below certain temperatures. 


oe 


—- — 
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8 A. Smekal, Zeits. f. techn, Physik 7, 535 (1926); 8, 561(1927) and reference 1. 
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LETTERS TO THE EDITOR 


Prompt publication of brief reporis of important discoveries in physics may 
be secured by addressing them to this department. Closing dates for this depart- 
ment are, for the first issue of the month, the twenty-eighth of the preceding month; 
for the second issue, the thirteenth of the month. The Board of Editors does not 
hold itself responsible for the opinions expressed by the correspondents. 


The Nature of Cosmic Radiation 





In a communication in Naturwissenschaf- 
ten [16, 1044, (1928) ] Bothe and Kohlhérster 
have described absorption experiments with 
the new Geiger tube counter which indicate 
that cosmic radiation is of a corpuscular na- 
ture. In a more detailed description of the 
experiments given in the same _ periodical 
[17, 271, (1929)] they estimate that if this 
radiation is of high speed electrons, it must 
have an original energy of at least 10° electron 
volts. 

Although the absorption experiments seem 
definitely to indicate that the view of these 
authors is correct, more direct information 
on this point could be secured from experi- 
ments of the same kind with a magnetic field 
instead of absorption screens. With a mag- 
netic field of 7000 gauss and 24 centimeters in 
diameter, it should be possible to bend elec- 
trons of 10* volts velocity so that most of 
those passing through a tube counter above 


the field would fail to pass through a similar 
counter directly below it and the magnetic 
field. Hence, the softer portion of the radi- 
ation would be cut off from the lower counter 
and the number of coincidences reduced, 
since the radiation comes mainly in a vertical 
direction. 

Experiments have been started with such 
an arrangement using an air gap of 5 cm 
between the pole pieces of the electromagnet. 
Preliminary experiments indicate a definite 
decrease in the number of coincidences ob- 
served when the magnetic field is used, in 
confirmation of the views of Bothe and 
Kohlhérster. More data, however, are needed 
to make this point certain and a full account 
will appear later in the Bureau of Standards 
Journal of Research. 

L. F. Curtiss 

Bureau of Standards, 

October 30, 1929. 


The Kinetic Correction in Fluid Flow 


When the viscosity of a liquid is deter- 
mined from the rate of efflux, (£), through a 
capillary, a correction is made for the kinetic 
energy of the liquid. There is no existing 
method of doing this if the fluid deviates from 
the law of Poiseuille. Since a portion of the 
pressure, (P), is used in imparting energy of 
motion, this suggests that the correction be 
applied to the pressure. The work done per 
second on the liquid is 


W=EPg 
The kinetic energy of the same liquid is 
T = [2° /P (rdrdepv)v*/2 


v is the velocity component parallel to the 
axis of the capillary. The energy used in over- 
coming internal friction is T—-W. The 


pressure used in overcoming internal friction 
is therefore 


(1) P? =P —( xp/Eg) {Pv'rdr 


If the efflux is known as a function of the 
pressure, say E=f(P), where the rate of flow 
is small; then this equation can be corrected 
for kinetic energy by writing E =f(P*). 

The method requires that v be known as a 
function of r. For many fluids there exists re- 
lationships such as, E= F(R). The equation 


E=- af Pr*(av/ar)dr 


defines a value of dv/dr. This together with 
the condition that v=0, at r=R, gives a ve- 
locity function which can be substituted in 
Eq. (1). The value of v/dr, (or v), is however 
not unique; for to any function thus found, 
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as v, there can be added any number of other 
functions such as y, provided that /yrdr =0. 
The added functions will not affect the efflux, 
since E=2xfvrdr, but they will alter the 
kinetic energy term. Hypothetical y functions 
can be invented which satisfy all the known 
physical limitations to which »v is subject. 
Because of this, even the agreement be- 
tween theory and experiment shown by 
liquids obeying the law of Poisseuille, does 
not prove that the velocity function is 
uniquely determined. 


The kinetic energy correction can therefore 
assume values other than those calculated by 
the usual methods. It is possible that the 
anomalous turbidity shown by ammonium 
oleate solutions, where turbulent flow starts 
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at low values of the Reynold’s number, may 
be associated with a large kinetic energy 
correction which originates in the above man- 
ner. 

If the above limitation is ignored, the cal- 
culation for a viscous liquid yields the ac- 
cepted kinetic energy term. For plastic flow, 
as defined by Bingham, the kinetic energy 
correction is identical with an empirical term 
which is often employed. 

This paper will be published in one of the 
early numbers of “Rheology” (rhe-flow) a new 
quarterly dealing with phenomena connected 
with the flow of matter. 

S. BRADFORD STONE 

College of the City of New York 

October 30, 1929. 


The Nature of the Electrodeless Ring Discharge 


There seems to be a misunderstanding as to 
the nature of the electric forces which produce 
an electrodeless ring discharge. It is com- 
monly thought that the only force involved 
is the electromagnetic force due to the varia- 
tion of the magnetic field and hence it is 
supposed that the electrons in the discharge 
move in circles around the axis of the solenoid. 
J. S. Townsend and R. H. Donaldson (Phil. 
Mag. 5, 178, 1928) have shown that the 
electric force due to electric charges on the 
solenoid is generally very large compared to 
the electromagnetic force. During the course 
of some work on an electrodeless ring dis- 
charge, the writer has run across two lines of 
evidence which point to the same conclusion. 

To produce the electrodeless discharge, the 
writer used an ordinary Colpitts circuit whose 
inductance, a solenoid of 15 turns, surrounded 
a cylindrical discharge tube containing hy- 
drogen at 0.03 mm pressure. With the circuit 
oscillating at 50 meters and a tank current 
of about 7 amperes, a relatively intense dis- 
charge was produced. In order to increase the 
wave-length, additional inductance consisting 
of a coil of 7 turns was inserted in series with 
the exciting coil. However, these additional 
7 turns did not surround the discharge tube. 
Under these conditions, the circuit oscillated 
at 100 meters with over 10 amperes in the 
tank circuit but mo discharge was produced. 
As soon as the 7 additional turns were 
wrapped around the discharge tube as an ex- 


tension of the original coil, a very nice dis- 
charge was obtained, and all other conditions 
were apparently the same as before. This 
then, seems to indicate that the functioning 
of the electrodeless ring discharge depends 
primarily not on the current through the exciting 
coil but upon the potential difference between 
the ends of the coil. 

Another simple experiment, which would 
show whether the electric forces producing an 
electrodeless ring discharge were electro- 
magnetic or electrostatic in nature, was sug- 
gested to the writer by Dr. G. Breit. The 
exciting coil was increased in diameter until 
there was about 2 cm clearance between coil 
and tube. Under these conditions, a discharge 
could still be obtained in the tube. Now a 
grounded cylinder of copper gauze closely fit- 
ting the discharge tube was inserted between 
the exciting coil and tube. This shield should 
have no effect on electromagetic forces in 
the discharge tube, but would completely 
screen it from any electrostatic forces. It was 
found that as soon as the tube was surrounded 
by the shield, the discharge ceased. Thus we 
must conclude that the electrodeless ring dis- 
charge is an electrostatic and not an electro- 
magnetic phenomenon. 


Cuar es J. BRASEFIELD 


University of Michigan, 
October 30, 1929. 
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Raman Lines in the Spectrum of the Electric Discharge 


A suggestion has been made by H. S. Allen! 
and extended by Deodhar,? that many of the 
faint lines in the secondary spectrum of hydro- 
gen arise from Ramanscattering in theelectric 
discharge. It is the purpose of this note to 
call attention to a possible relationship be- 
tween the results of Allen and of Deodhar and 
other observations on the secondary spectrum 
and also to discuss several difficulties that 
might prevent the acceptance of the above 
suggestion. 

It may be recalled that Glitscher® plotted 
the intensity distribution in the secondary 
spectrum of hydrogen and pointed out that 
the frequency separations between the most 
intense maxima in the intensity curves and 
the Balmer lines, were approximately con- 
stant. This relationship was found to hold 
also for a group of secondary maxima and 
possibly for even lesser maxima. It has been 
pointed out since the observation of this fact 
that some theoretical significance may at some 
time be ascribed to it. It is interesting now to 
note that the work of Allen points out a simi- 
lar relationship between the Balmer lines and 
some of the weak secondary lines. Both Allen 
and Deodhar have pointed out some very re- 
markable agreement between the frequency 
separations of some of the lines in the sec- 
ondary spectrum of hydrogen and the Balmer 
lines. Furthermore they have found several 
groups of lines whose frequency separations 
from the Balmer lines are integral multiples of 
the lowest frequency separation. Although 
it is very tempting to call these lines Raman 
lines the frequency separations do not seem 
to be connected with anything that is known 
about the hydrogen molecule. There are other 
objections to the suggestion that the lines 
originate in Raman scattering and those ob- 
jections will now be discussed. 

The first of these objections arises when 
one compares the conditions in a discharge 
tube with those present in Rasetti’s* Raman 
scattering experiments in hydrogen. The 
pressure of hydrogen in Rasetti’s experiments 
was about 12 atmospheres and the exposure 
times varied from 10 to 40 hours. The pres- 
sure used was thus about 20,000 times as 
great as that in a discharge tube and the 
exposure times were about 1000 times as 
long as those usually employed in photograph- 


ing discharge-tube spectra. These factors 
alone should make it very difficult to see how 
Raman scattering that takes place in a dis- 
charge at low pressures could ever be ob- 
served. 

Granting however that there is some special 
and powerful scattering process in discharge 
tubes, one can mention another difficulty that 
arises when the intensity of these weak lines 
is considered. Since the frequency separations 
of these lines from the parent Balmer lines 
are integral multiples of a lowest frequency 
separation, it is reasonable to say that the 
energy levels that enter into the Raman 
effect combine with a common electronic 
level. One would then expect that if one of 
the Raman lines scattered by Ha, for example, 
is stronger than the corresponding one scat- 
tered by Hg, then all of the lines originating 
from Ha should be more intense than those 
originating from Hg. This, however, is not 
true as a glance at Allen’s and Deodhar’s 
notes will show. Furthermore, a group of lines 
should be scattered by the Lyman series and 
those scattered lines should be more intense 
than the ones scattered by the Balmer series. 
This follows from the fact that the frequencies 
of the Lyman series are nearer to the elec- 
tronic frequencies of the molecule than are 
those of the Balmer lines. From the theory 
of the intensity of Raman scattering the 
above conclusion follows immediately. 

Although the above arguments make it 
difficult to accept the notion that the weak 
secondary lines of Allen and of Deodhar are 
Raman lines, it seems that the relationship 
between them and the Balmer lines may 
prove of considerable interest. This is es- 
pecially true because on a fine scale, it is ex- 
actly the same as the one discussed by 
Glitscher a long time ago. 


JoserH KAPLAN 


University of California 
at Los Angeles, 


November 4, 1929. 


1 Allen, Nature, January 26, 1929. 

2 Deodhar, Zeits. f. Physik 57, 570 (1929). 

’ Glitscher, Sitzungsber. d. bayer. Akad., 
1916, p. 125. 

‘F. Rasetti, Phys. Rev. 34, 367 (1929). 
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To the Editor of the Physical Review 


Dear Sir: 

It is not until now that I have learned that 
by an oversight the letter 7 was used instead 
of J for the inner-quantum number in my re- 
cent paper on “The Arc Spectrum of Nickel”’ 
(Phys. Rev. 34, 821 (1929)). I write to place 


it on record that this was a slip, resu!ting from 
haste in preparation of the final manuscript, 
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just before my departure for a year’s leave, 
and by no means a deliberate departure from 
the recommendation of the informal commit- 
tee on nomenclature, of which I was a mem- 
ber. 
Very truly yours, 
Henry Norris RUSSELL 


November 1, 1929. 
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Infra-Red Analysis of Molecular Structure. F. I. G. RAwitns AND A. M. Taytor. Pp. 
176, figs. 35. Cambridge University Press, London. 1929. 

In recent years a great many papers both theoretical and experimental have appeared on 
the subject of infra-red spectra and in this book Rawlins and Taylor have performed the wel- 
come service of reviewing and examining critically the more important of them. They have 
divided their book into the four chapters of Gases, Liquids, Solids, and Experimental. Of 
these the chapters on Liquids and Solids are especially interesting and the authors have at- 
tempted to dispose of many of the fantastic theories which have sprung up from time to time, 
and have suggested new experiments which should yield fruitful results. 

Despite its excellence as a critical review of infra-red spectroscopy, the book lacks a cer- 
tain organic unity and does not seem to emphasize sufficiently the essential points in the 
theory of molecular spectra. The Mathematical Appendix in particular does not seem to be 
adequate but perhaps a fuller treatment would have been outside the scope of the book. 

Principally, however, the reviewer believes that it will prove valuable to specialists in 
the field of the infra-red, both because of its excellent bibliography and because of the -inter- 
correlation which is given between widely scattered pieces of experimental work. 

D. M. DENNISON 


Speech and Hearing. Harvey FLetcHer. Pp. 331, 151 figures. D. Van Nostrand Com- 
pany, Inc. New York, 1929. Price $5.50. 

This book is another of the series which has resulted from the “out-of hours courses” 
presented at the Bell Laboratories. It summarizes in part the results obtained during the past 
fifteen years in the comprehensive and fundamental survey and research program of the Bell 
Laboratories. In this work Dr. Fletcher has been the director. 

The book is divided into four parts. In the first part, entitled “Speech,” after brief de- 
scriptions of the organs of speech, of the formation of the English speech sounds, and of the 
methods of recording speech waves, a number of analyses of speech sounds are given. These 
give the power distribution as well as frequencies, or acoustic spectra. 

Part two, “Music and Voice” contains chapters on the “Physical Properties of Musical 
Sounds” and “Noise.” In these are given the wave forms and frequency spectra of some 
muscial instruments; also methods of measuring noise. 

Part three, “Hearing” opens with a description of the organs of hearing. A theory of hear- 
ing is given based on that originally given by Helmholtz and modified and extended by recent 
investigation. It may be called a resonance theory but one in which the effects of friction are 
fully recognized. Quoting, “The pitch or position of a tone on the muscial scale is then de- 
pendent upon the position of the maximum stimulation along the basilar membrane. The 
character of the sound depends upon the position and relative intensities of agitation of the 
various parts of the basilar membrane.” 

In the chapters on “Limits of Audition,” “Minimum Perceptible Differences in Sound” 
and “Masking Effects” the experimental results obtained are well summarized in graphs. 
The closing chapters of the third part treat “Binanal Effects” and “Methods of Testing the 
Acuity of Hearing.” In this chapter are described tests of value to the physician and in the 
public schools. 

In part four, “The Preception of Speech and Music” the loudness of sensation level of 
pure tones and complex sounds is first considered. In the chapter on “The Recognition of 
Pitch of Musical Tones” the pitch of a musical sound is defined as being that of a pure tone 
which has been adjusted to give the same sensation of pitch to the ear. By the use of a high 
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quality telephone system into which various filters could be introduced some interesting 
results on pitch recognition have been obtained. It was found that the pitch is determined by 
the common difference in the frequency of the harmonics, but only when these harmonics are 
multiples of the common difference. If this condition is not satisfied, no definite pitch can be 
recognized. These results have a bearing on the theory of hearing. 

The last chapters describe methods of measuring the recognition of speech sounds and 
the effects of changes in intensity with and without frequency distortion; and the effects of 
other distortions and of noise on the recognition of speech sounds. In a very interesting group 
of curves are shown the relations between articulation and energy and frequency. 

This book is a valuable contribution to our knowledge of “Speech and Hearing,” as it com- 
piles and records the recent work in that field. It should be of value to the physician and the 
psychologist as well as to the teacher of speech, the telephone engineer and the physicist. 

J. W. Bucuta 


Television. H. Horton SHELDON AND E, N. GrisEwoop. 129 figures. Pp. 194. D. Van 
Nostrand Company, Inc. New York, 1929. Price $2.75. 

The basic phenomena of television are described in a simple non-mathematical language. 
Then follow descriptions of the apparatus that is common to most systems such as the photo- 
electric cell, glow lamps and scanning and synchronizing devices. The particular arrangements 
used by the David, Bell, Jenkins and Alexanderson systems are described in separate chapters. 
One chapter is devoted to telephotography. The book contains a large number of photographs. 

J. W. Bucuta 





